No abstract

We revisit the CROS (“CRitical OScillations”) model which was recently proposed as an attempt to reproduce both scale-invariant neuronal avalanches and long-range temporal correlations. With excitatory and inhibitory stochastic neurons locally connected in a two-dimensional disordered network, the model exhibits a transition where alpha-band oscillations emerge. Precisely at the transition, the fluctuations of the network activity have nontrivial detrended fluctuation analysis (DFA) exponents, and avalanches (defined as supra-threshold activity) have power law distributions of size and duration. We show that, differently from previous results, the exponents governing the distributions of avalanche size and duration are not necessarily those of the mean-field directed percolation universality class (3/2 and 2, respectively). Instead, in a narrow region of parameter space, avalanche exponents obtained via a maximum-likelihood estimator vary continuously and follow a linear relation, in good agreement with results obtained from M/EEG data. In that region, moreover, the values of avalanche and DFA exponents display a spread with positive correlations, reproducing human MEG results.

Since the first measurements of neuronal avalanches [1], the critical brain hypothesis has gained traction [2]. However, if the brain is critical, what is the phase transition? For several decades it has been known that the cerebral cortex operates in a diversity of regimes [3], ranging from highly synchronous states (e.g. slow wave sleep [4], with higher spiking variability) to desynchronized states (e.g. alert waking [5], with lower spiking variability). Here, using independent signatures of criticality, we show that a phase transition occurs in an intermediate value of spiking variability. The critical exponents point to a universality class different from mean-field directed percolation (MF-DP). Importantly, as the cortex hovers around this critical point [6], it follows a linear relation between the avalanche exponents that encompasses previous experimental results from different setups [7,8] and is reproduced by a model. * AJF and NAPV contributed equally. †

Population models are a powerful mathematical tool to study the dynamics of neuronal networks and to simulate the brain at macroscopic scales. We present a mean-field model capable of quantitatively predicting the temporal dynamics of a network of complex spiking neuronal models, from Integrate-and-Fire to Hodgkin–Huxley, thus linking population models to neurons electrophysiology. This opens a perspective on generating biologically realistic mean-field models from electrophysiological recordings.