Global dynamics of neural mass models

Neural mass models are used to simulate cortical dynamics and to explain the electrical and magnetic fields measured using electro- and magnetoencephalography. Simulations evince a complex phase-space structure for these kinds of models; including stationary points and limit cycles and the possibility for bifurcations and transitions among different modes of activity. This complexity allows neural mass models to describe the itinerant features of brain dynamics. However, expressive, nonlinear neural mass models are often difficult to fit to empirical data without additional simplifying assumptions: e.g., that the system can be modelled as linear perturbations around a fixed point. In this study we offer a mathematical analysis of neural mass models, specifically the canonical microcircuit model, providing analytical solutions describing dynamical itinerancy. We derive a perturbation analysis up to second order of the phase flow, together with adiabatic approximations. This allows us to describe amplitude modulations as gradient flows on a potential function of intrinsic connectivity. These results provide analytic proof-of-principle for the existence of semi-stable states of cortical dynamics at the scale of a cortical column. Crucially, this work allows for model inversion of neural mass models, not only around fixed points, but over regions of phase space that encompass transitions among semi or multi-stable states of oscillatory activity. In principle, this formulation of cortical dynamics may improve our understanding of the itinerancy that underwrites measures of cortical activity (through EEG or MEG). Crucially, these theoretical results speak to model inversion in the context of multiple semi-stable brain states, such as onset of seizure activity in epilepsy or beta bursts in Parkinsons disease.