Volume 10 Supplement 1

Eighteenth Annual Computational Neuroscience Meeting: CNS*2009

Open Access

Analysis of the power spectra, autocorrelation function and EEG time-series signal of a network of leaky integrate-and-fire neurons with conductance-based synapses

  • Andre DH Peterson1, 2Email author,
  • Hamish Meffin4,
  • Anthony N Burkitt1, 2,
  • Iven MY Mareels1,
  • David B Grayden1, 2,
  • Levin Kuhlmann1 and
  • Mark J Cook2, 3
BMC Neuroscience200910(Suppl 1):P167

DOI: 10.1186/1471-2202-10-S1-P167

Published: 13 July 2009


Focal epilepsy is characterized by the spread of seizure activity from pathological cortical tissue (focus) to other parts of the surrounding cortex and is typically diagnosed via the EEG [1]. The research described below will form the basis of a mathematical description of a mesoscopic network of cortical columns where the network dynamics will be examined as seizure-like behaviour spreads from a focal (pathological) column to other columns. In particular, the emphasis will be on how the local dynamics, network topology and physiological regulatory (control) mechanisms influence the overall global dynamics of the seizure spread. This study examines the dynamics of a network of neurons that approximate a single cortical column. Both the time series and power spectrum of the network are calculated and used to approximate the EEG signal of a cortical column.


The power spectrum is calculated from the autocorrelation function of a network of leaky integrate-and-fire neurons with conductance-based synapses that receive Poisson distributed synaptic input [2]. This is then generalized to a mean-field network approximation that includes both excitatory and inhibitory neurons [3]. This requires the calculation of the first passage time density, which is found numerically by solving a nonlinear Volterra integral of the first kind using Fourier transform methods. This also yields the EEG time-series resulting from the spikes generated by the network.


The analytical results of the power spectra, autocorrelation function, first-passage time density and EEG time series are compared with network simulation results. Results were obtained using parameter values that represent typical cortical in vivo neurons [4].


This work is the first stage necessary for constructing a physiologically plausible mathematical model of a mesoscopic network of cortical columns. The results presented here will be used as a mathematical approximation of a single cortical column and be generalized into a nonlinear network of columns. Future research will be directed toward incorporating an epileptic focus into a network of columns to investigate seizure propagation dynamics as described in the introduction.



This work was funded by the Australian Research Council (ARC Linkage Project #LP0560684).

Authors’ Affiliations

Department of Electrical & Electronic Engineering, The University of Melbourne
The Bionic Ear Institute
Department of Clinical Neurosciences, St. Vincent's Hospital
NICTA VRL, c/- Dept of Electrical & Electronic Engineering, University of Melbourne


  1. Milton J, Jung P: Epilepsy as a Dynamic Disease. 2003, SpringerView ArticleGoogle Scholar
  2. Burkitt AN: A review of the integrate-and-fire neuron model: II. Inhomogeneous synaptic input and network properties. Biol Cybern. 2006, 95: 97-112. 10.1007/s00422-006-0082-8.PubMedView ArticleGoogle Scholar
  3. Brunel N: Dynamics of sparsely connected networks of excitatory and inhibitory spiking neurons. J Comput Neurosci. 2000, 8: 183-208. 10.1023/A:1008925309027.PubMedView ArticleGoogle Scholar
  4. Meffin H, Burkitt AN, Grayden DB: An analytical model for the 'large, fluctuating synaptic conductance state' typical of neocortical neurons in vivo. J Comput Neurosci. 2004, 16: 159-175. 10.1023/B:JCNS.0000014108.03012.81.PubMedView ArticleGoogle Scholar


© Peterson et al; licensee BioMed Central Ltd. 2009

This article is published under license to BioMed Central Ltd.