Email updates

Keep up to date with the latest news and content from BMC Neuroscience and BioMed Central.

This article is part of the supplement: Twenty First Annual Computational Neuroscience Meeting: CNS*2012

Open Access Poster presentation

Turing instabilities in a mean field model of electrocortical activity

Lennaert van Veen* and Kevin Green

Author Affiliations

Faculty of Science, University of Ontario Institute of Technology, Oshawa, Ontario L1H 7K4, Canada

For all author emails, please log on.

BMC Neuroscience 2012, 13(Suppl 1):P128  doi:10.1186/1471-2202-13-S1-P128

The electronic version of this article is the complete one and can be found online at:

Published:16 July 2012

© 2012 van Veen and Green; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Poster presentation

The mean field model formulated by Liley et al. [1] describes the membrane potential and synaptic interaction of excitatory and inhibitory neuron populations in a two-dimensional slab of cortical tissue. When complemented with periodic boundary conditions, it yields a set of fourteen coupled partial differential equations that encompass both highly nonlinear local interaction and long-range interaction through cable equations. We consider these equations as an autonomous dynamical system and aim to parse it using bifurcation analysis. We focus in particular on Turing-type bifurcations from spatially homogeneous equilibria to time-periodic patterns. Such bifurcations have been considered in simplified models but, to the best of our knowledge, not in a physiologically realistic setting [2].

In this study, we fix the model parameters to values that have previously been shown to lead to physiologically interesting γ-range activity [3]. We locate the primary instability when varying the inhibitory-to-inhibitory connection density, known to strongly affect the power spectrum, and the system size.

Although for small system size, spatially homogeneous periodic solutions appear, for system sizes of a few square centimeters we observe centimeter-scale, oscillating patterns, which appear in subcritical Hopf bifurcations and have a strong γ-band component in their power spectrum. An example is shown in Fig. 1, for a domain of 12.8 by 12.8cm. The black curve denotes stable (solid) and unstable (dashed) equilibria, and the blue curve denotes spatially homogeneous periodic solutions. The latter can be stable to spatially homogeneous perturbations. The primary instability corresponds to a spatial wave number two in both directions. The inlay shows a snapshot of the excitatory membrane potential for the corresponding, mildly unstable, periodic solution near the Hopf bifurcation.

thumbnailFigure 1. Partial bifurcation diagram of the mean-field model. A periodic solution (blue) bifurcates from an equilibrium (black), both spatially homogeneous. The primary instability, however, is a Turing bifurcation to the pattern in the left inlay. The second inlay shows the power spectrum of the excitatory membrane potential at one point in the domain, with a peak at 40Hz as well as resonance peaks. The periodic solution was computed with 1mm resolution.

We have implemented a time-stepper in the open-source PETSc environment [4]. It uses a finite difference approximation for the Laplacian and first order implicit time stepping, both for the field equations and for the tangent linear model. The distribution of the grid over CPUs is automated, allowing for efficient parallel computation. This code can be used for the computation of equilibria, traveling waves, periodic orbits, manifolds and other building blocks of dynamical systems analysis, leading to a greater insight in the model's behaviour than could previously be extracted from analysis of spatially homogeneous solutions alone.


The authors would like to thank Ingo Bojak and David Liley for useful discussions.


  1. Liley DTJ, Cadusch PJ, Dafilis MP: A spatially continuous mean field theory of electrocortical activity.

    Network: comput Neural syst 2002, 13:67-113. OpenURL

  2. Coombes S: Large-scale neural dynamics: Simple and complex.

    NeuroImage 2010, 52:731-739. PubMed Abstract | Publisher Full Text OpenURL

  3. Bojak I, Liley DTJ: Self-organized 40Hz synchronization in a physiological theory of EEG.

    Neurocomput 2007, 70:2085-2090. Publisher Full Text OpenURL

  4. Balay S, et al.: PETSc users manual. [] webcite

    Argonne National Laboratory OpenURL