Skip to main content


We're creating a new version of this page. See preview

  • Poster presentation
  • Open Access

Biologically plausible statistics from a Markov model of spiking cortical networks

  • 1Email author,
  • 2,
  • 3,
  • 2 and
  • 1
BMC Neuroscience20089 (Suppl 1) :P38

  • Published:


  • Critical Exponent
  • Neuron Network
  • Spike Train
  • Cortical Slice
  • Stochastic Simulation Algorithm


Understanding the spatio-temporal correlations observed in nervous tissue is a major challenge in computational neuroscience. We approach this challenge by modeling the global network state as a single Markov process [1, 2]. Each neuron is modeled as a two (active or inactive) – or three-state (active, inactive, or refractory) random variable, with each neuron's spike probability a function of its input current and internal threshold, updated at time steps dependent upon the current network state. Using the stochastic simulation algorithm [3], we simulate the network model with excitatory and inhibitory neurons and with random connectivity tuned so that each spike triggers on average one new spike, i.e. in the critical regime.

In order to study the model analytically, we recast the model in the language of statistical field theory [2]. The result is a moment-generating functional which permits the calculation of correlation functions and critical exponents. Analytical calculation of the first moment produces a form of the Wilson-Cowan equations [4]. The analytical and simulation results can be compared to those obtained from cortical slice recordings [5].


Figure 1 shows that model simulations produce spike trains, as part of highly variable, correlated, network activity resembling the "neuronal avalanches" observed in cortical slice tissue [5]. The simulations reveal a power-law distribution of avalanche sizes with a critical exponent of -2.15 when a 900 neuron network is simulated (data set represented by dots in Figure 2) and -2.168 when a 1600 neuron network is simulated (data set represented by circles in Figure 2). These compare well with the critical exponents calculated from unsorted spike data taken from multiunit recordings. Figure 3 shows one such experiment, where the critical exponent was calculated to be -2.14.
Figure 1
Figure 1

Raster Plot.

Figure 2
Figure 2

Model size distribution.

Figure 3
Figure 3

Experimental size distribution.



This work was supported in part by the Falk Foundation and the Linn family. Also, we would like to thank Dr. John Beggs, Dr. Alan Litke, and Aonan Tang for providing the experimental data on critical exponents and Figure 3.

Authors’ Affiliations

Department of Pediatrics, University of Chicago, Chicago, IL 60637, USA
Department of Mathematics, University of Chicago, Chicago, IL 60637, USA
The Salk Institute for Biological Studies, La Jolla, CA 92037, USA


  1. Cowan JD: Stochastic neurodynamics. Proc of the 1990 conf on adv in neural inf process sys 3. 1990, 62-69.Google Scholar
  2. Buice MA, Cowan JD: Field-theoretic approach to fluctuation effects in neural networks. Phys Rev E. 2007, 75 (5): 051919-10.1103/PhysRevE.75.051919.View ArticleGoogle Scholar
  3. Gillespie DT: Exact stochastic simulation of coupled chemical reactions. J Phys Chem. 1977, 81 (25): 2340-2361. 10.1021/j100540a008.View ArticleGoogle Scholar
  4. Wilson HR, Cowan JD: Excitatory and inhibitory interactions in localized populations of model neurons. Biophys J. 1972, 12 (1): 1-24.PubMed CentralView ArticlePubMedGoogle Scholar
  5. Beggs JM, Plenz D: Neuronal avalanches in neocortical circuits. J Neurosci. 2003, 23 (35): 11167-11177.PubMedGoogle Scholar


© Benayoun et al; licensee BioMed Central Ltd. 2008

This article is published under license to BioMed Central Ltd.