A population study of integrate-and-fire-or-burst neurons.

Any realistic model of the neuronal pathway from the retina to the visual cortex (V1) must account for the bursting behavior of neurons in the lateral geniculate nucleus (LGN). A robust but minimal model, the integrate-and-fire-or-burst (IFB) model, has recently been proposed for individual LGN neurons. Based on this, we derive a dynamic population model and study a population of such LGN cells. This population model, the first simulation of its kind evolving in a two-dimensional phase space, is used to study the behavior of bursting populations in response to diverse stimulus conditions.

A population approach to cortical dynamics with an application to orientation tuning.

A typical functional region in cortex contains thousands of neurons, therefore direct neuronal simulation of the dynamics of such a region necessarily involves massive computation. A recent efficient alternative formulation is in terms of kinetic equations that describe the collective activity of the whole population of similar neurons. A previous paper has shown that these equations produce results that agree well with detailed direct simulations. Here we illustrate the power of this new technique by applying it to the investigation of the effect of recurrent connections upon the dynamics of orientation tuning in the visual cortex. Our equations express the kinetic counterpart of the hypercolumn model from which Somers et al (Somers D, Nelson S and Sur M 1995 J. Neurosci. 15 5448-65) computed steady-state cortical responses to static stimuli by direct simulation. We confirm their static results. Our method presents the opportunity to simulate the data-intensive dynamical experiments of Ringach et al (Ringach D, Hawken M and Shapley R 1997 Nature 387 281-4), in which 60 randomly oriented stimuli were presented each second for 15 min, to gather adequate statistics of responses to multiple presentations. Without readjustment of the previously defined parameters. our simulations yield substantial agreement with the experimental results. Our calculations suggest that differences in the experimental dynamical responses of cells in different cortical layers originate from differences in their recurrent connections with other cells. Thus our method of efficient simulation furnishes a variety of information that is not available from experiment alone.

The approach of a neuron population firing rate to a new equilibrium: an exact theoretical result.

The response of a noninteracting population of identical neurons to a step change in steady synaptic input can be analytically calculated exactly from the dynamical equation that describes the population's evolution in time. Here, for model integrate-and-fire neurons that undergo a fixed finite upward shift in voltage in response to each synaptic event, we compare the theoretical prediction with the result of a direct simulation of 90,000 model neurons. The degree of agreement supports the applicability of the population dynamics equation. The theoretical prediction is in the form of a series. Convergence is rapid, so that the full result is well approximated by a few terms.

On the simulation of large populations of neurons.

The dynamics of large populations of interacting neurons is investigated. Redundancy present in subpopulations of cortical networks is exploited through the introduction of a probabilistic description. A derivation of the kinetic equations for such subpopulations, under general transmembrane dynamics, is presented. The particular case of integrate-and-fire membrane dynamics is considered in detail. A variety of direct simulations of neuronal populations, under varying conditions and with as many as O(10(5)) neurons, is reported. Comparison is made with analogous kinetic equations under the same conditions. Excellent agreement, down to fine detail, is obtained. It is emphasized that no free parameters enter in the comparisons that are made.