- Poster presentation
Turing instabilities in a mean field model of electrocortical activity
BMC Neurosciencevolume 13, Article number: P128 (2012)
The mean field model formulated by Liley et al.  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 .
In this study, we fix the model parameters to values that have previously been shown to lead to physiologically interesting γ-range activity . 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.
We have implemented a time-stepper in the open-source PETSc environment . 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.
Liley DTJ, Cadusch PJ, Dafilis MP: A spatially continuous mean field theory of electrocortical activity. Network: comput Neural syst. 2002, 13: 67-113.
Coombes S: Large-scale neural dynamics: Simple and complex. NeuroImage. 2010, 52: 731-739. 10.1016/j.neuroimage.2010.01.045.
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.
Balay S, et al: PETSc users manual. Argonne National Laboratory. [http://www.ncs.anl.gov/petsc]
The authors would like to thank Ingo Bojak and David Liley for useful discussions.