A typical neuron, e.g., a pyramidal cell in the vertebrate cortex, receives input spikes from thousands of other neurons, which in turn receive input from their presynaptic neurons and so forth; see Fig. 8.4. While it is not impossible to incorporate millions of integrate-and-fire neurons into a huge network model, it is often reasonable to focus the modeling efforts on a specific subset of neurons, e.g., a column in the visual cortex, and describe input from other parts of the brain as a stochastic background activity.
Let us consider a nonlinear integrate-and-fire neuron with index that is part of a large network. Its input consists of (i) an external input ; (ii) input spikes from other neurons of the network; and (iii) stochastic spike arrival due to the background activity in other parts of the brain. The membrane potential evolves according to
where is the Dirac function and is the coupling strength from a presynaptic neurons in the network to neuron . Input from background neurons is weighted by the factor . While the firing times are generated by the threshold crossings of presynaptic integrate-and-fire neurons, the firing times of a background neuron are generated by a Poisson process with mean rate .
To simplify the following discussions we adopt three simplifications. First, we focus on a leaky integrate-and-fire neuron and shift the voltage so that the resting potential is at zero. Hence we can set . Second, we concentrate on a single neuron receiving stochastic input from background neurons. Hence we can drop the sum over which represents input from the network and also drop the index of our specific neuron. We therefore arrive at
In Stein’s model, each input spike generates a postsynaptic potential with , i.e., the potential jumps upon spike arrival by an amount and decays exponentially thereafter. Integration of Eq. (8.20) yields
for where is the last firing time of the neuron. It is straightforward to generalize the model so as to include a synaptic time constant and work with arbitrary postsynaptic potentials that are generated by stochastic spike arrival; cf. Fig. 8.5A.
In order to calculate the fluctuations of the membrane potential caused by stochastic spike arrival, we assume that the firing threshold is relatively high and the input weak so that the neuron does not reach its firing threshold. Hence, we can safely neglect both threshold and reset. The leaky integrate-and-fire model of Stein (Eq. (8.21)) is then equivalent to a model of a passive membrane driven by stochastic spike arrival.
We assume that each input spike evokes a postsynaptic potential of the same amplitude and shape, independent of . The input statistics is assumed to be Poisson, i.e., firing times are independent. Thus, the total input spike train (summed across all synapses)
that arrives at neuron is a random process with expectation
cf. Eq. (7.46).
Obviously, the integration over the -function in the last term on the right-hand side is possible and would lead back to the more compact representation . The advantage of having the spike train appear explicitly is that we can exploit the definition of the random process , in particular, its mean and variance.
In Fig. 8.5 we have simulated a neuron which receives input from background neurons with rate Hz. The total spike arrival rate is therefore 1 kHz. Each spike evokes an EPSP with ms. The evaluation of Eqs. (8.26) and (8.2.1) for constant input yields and .
Let us suppose that an integrate-and-fire neuron defined by Eq. (8.20) with ms receives input from 100 excitatory neurons () and 100 inhibitory neurons (). Each background neuron fires at a rate of = 10 Hz. Thus, in each millisecond, the neuron receives on average one excitatory and one inhibitory input spike. Each spike leads to a jump of the membrane potential of . The trajectory of the membrane potential is therefore similar to that of a random walk subject to a return force caused by the leak term that drives the membrane potential always back to zero; cf. Fig. 8.6A.
If, in addition, a constant stimulus is applied so that the mean membrane potential (in the absence of the background spikes) is just below threshold, then the presence of random background spikes may drive toward the firing threshold. Whenever , the membrane potential is reset to .
Since firing is driven by the fluctuations of the membrane potential, the interspike intervals vary considerably; cf. Fig. 8.6. Balanced excitatory and inhibitory spike input could thus contribute to the large variability of interspike intervals in cortical neurons; see Section 8.3.
With the above set of parameters, the mean of the stochastic background input vanishes since . Using the same arguments as in the previous example, we can convince ourselves that the stochastic arrival of background spikes generates fluctuations of the voltage with variance
cf. Section 8.4 for a different derivation.
Let us now increase all rates by a factor of and multiply at the same time the synaptic efficacies by a factor . Then both mean and variance of the stochastic background input are the same as before, but the size of the jumps is decreased; cf. Fig. 8.6B. In the limit of the jump process turns into a diffusion process and we arrive at the stochastic model of Eq. (8.7). In other words, the balanced action of the excitatory and inhibitory spike trains, and respectively, arriving at the synapses with Poisson input rate yields in the limit a white noise input
The above transition is called the diffusion limit and will be systematically discussed in Section 8.4. Intuitively, the limit process implies that in each short time interval a large number of excitatory and inhibitory input spikes arrive, each one causing the membrane potential to jump by a tiny amount upward or downward.
© Cambridge University Press. This book is in copyright. No reproduction of any part of it may take place without the written permission of Cambridge University Press.