Multiple oscillatory states in models of collective neuronal dynamics.

Abstract

In our previous studies, we showed that the both realistic and analytical computational models of neural dynamics can display multiple sustained states (attractors) for the same values of model parameters. Some of these states can represent normal activity while other, of oscillatory nature, may represent epileptic types of activity. We also showed that a simplified, analytical model can mimic this type of behavior and can be used instead of the realistic model for large scale simulations. The primary objective of the present work is to further explore the phenomenon of multiple stable states, co-existing in the same operational model, or phase space, in systems consisting of large number of interconnected basic units. As a second goal, we aim to specify the optimal method for state control of the system based on inducing state transitions using appropriate external stimulus. We use here interconnected model units that represent the behavior of neuronal populations as an effective dynamic system. The model unit is an analytical model (S. Kalitzin et al., Epilepsy Behav.22 (2011) S102-S109) and does not correspond directly to realistic neuronal processes (excitatory-inhibitory synaptic interactions, action potential generation). For certain parameter choices however it displays bistable dynamics imitating the behavior of realistic neural mass models. To analyze the collective behavior of the system we applied phase synchronization analysis (PSA), principal component analysis (PCA) and stability analysis using Lyapunov exponent (LE) estimation. We obtained a large variety of stable states with different dynamic characteristics, oscillatory modes and phase relations between the units. These states can be initiated by appropriate initial conditions; transitions between them can be induced stochastically by fluctuating variables (noise) or by specific inputs. We propose a method for optimal reactive control, allowing forced transitions from one state (attractor) into another.