Volume 13 Supplement 1

Twenty First Annual Computational Neuroscience Meeting: CNS*2012

Open Access

Turing instabilities in a mean field model of electrocortical activity

BMC Neuroscience201213(Suppl 1):P128

DOI: 10.1186/1471-2202-13-S1-P128

Published: 16 July 2012

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.
Figure 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.

Authors’ Affiliations

Faculty of Science, University of Ontario Institute of Technology


  1. Liley DTJ, Cadusch PJ, Dafilis MP: A spatially continuous mean field theory of electrocortical activity. Network: comput Neural syst. 2002, 13: 67-113.View ArticleGoogle Scholar
  2. Coombes S: Large-scale neural dynamics: Simple and complex. NeuroImage. 2010, 52: 731-739. 10.1016/j.neuroimage.2010.01.045.View ArticlePubMedGoogle Scholar
  3. Bojak I, Liley DTJ: Self-organized 40Hz synchronization in a physiological theory of EEG. Neurocomput. 2007, 70: 2085-2090. 10.1016/j.neucom.2006.10.087.View ArticleGoogle Scholar
  4. Balay S, et al: PETSc users manual. Argonne National Laboratory. [http://www.ncs.anl.gov/petsc]


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

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.