Chaos and generalised multistability in a mesoscopic model of the electroencephalogram

We present evidence for chaos and generalised multistability in a mesoscopic model of the electroencephalogram (EEG). Two limit cycle attractors and one chaotic attractor were found to coexist in a two-dimensional plane of the ten-dimensional volume of initial conditions. The chaotic attractor was found to have a moderate value of the largest Lyapunov exponent (3.4 s⁻¹ base e) with an associated Kaplan-Yorke (Lyapunov) dimension of 2.086. There are two different limit cycles appearing in conjunction with this particular chaotic attractor: one multiperiodic low amplitude limit cycle whose largest spectral peak is within the alpha band (8-13 Hz) of the EEG; and another multiperiodic large-amplitude limit cycle which may correspond to epilepsy. The cause of the coexistence of these structures is explained with a one-parameter bifurcation analysis. Each attractor has a basin of differing complexity: the large-amplitude limit cycle has a basin relatively uncomplicated in its structure while the small-amplitude limit cycle and chaotic attractor each have much more finely structured basins of attraction, but none of the basin boundaries appear to be fractal. The basins of attraction for the chaotic and small-amplitude limit cycle dynamics apparently reside within each other. We briefly discuss the implications of these findings in the context of theoretical attempts to understand the dynamics of brain function and behaviour.