Critical behaviour of the stochastic Wilson-Cowan model

Inhibitory neurons and the asymmetric shape of neuronal avalanches

In the last twenty years neuronal avalanches have been deeply investigated, both experimentally and numerically, also framing the results in the context of the avalanche scaling theory. In particular the avalanche shape has recently received a wide attention, also because the existence of a universal shape is an indication of the brain acting at a critical point. Within this scope, the detection of the shape asymmetry and the understanding of the mechanisms leading to it can provide useful insights into brain activity. Experimental data evidence, either symmetric or leftward asymmetry in the shape, results are not confirmed by numerical studies. Here we analyze the role of inhibition, connectivity range, and short term plasticity in determining the avalanche shape in an integrate and fire model. Results indicate that, not only the physiological fraction of inhibitory neurons is crucial to observe leftward asymmetry, but also the different synaptic recovery rates between excitatory and inhibitory neurons, confirming the importance of a dynamic balance between excitation and inhibition in brain activity.

Threshold crossing time theory for quasicycles with application to brain rhythms

The amplitude of a two-dimensional Ornstein-Uhlenbeck colored noise process evolves according to the one-dimensional Rayleigh process. This is a general model for the random amplitude fluctuations of a quasicycle, i.e., of a noise-induced oscillation around an equilibrium with complex eigenvalues in physical and biological systems. We consider the probability density of time intervals during which the amplitude is either below or above a fixed threshold. The statistics of such first return times (FRTs) are of particular interest in neuroscience to characterize brain rhythm power excursions known as bursts, as well as avalanches and other branching processes. In contrast with the density of first passage times computed using Fokker-Planck theory between a start point and a different endpoint, the density of FRTs is non-normalizable. A recently proposed technique reframes the problem using an expansion of the Fokker-Planck eigenfunctions along with a correction to the normalization. Analytical expressions for the FRT density for the Rayleigh process are shown to be in good agreement with those computed from numerical realizations over a wide range of parameters, both for trajectories above and below threshold. Special care is required to evaluate the theory above threshold due to the crowded roots of the Tricomi confluent hypergeometric function. The results provide insight into the statistics of threshold crossing times in quasicycles generally, and in the stochastic Wilson-Cowan neural equations in particular. Surprisingly, FRTs are governed by a single meta-parameter Δ given by the ratio of the noise strength and the linear stability coefficient. We find the universal property that the mean FRT is invariant to the ratio of threshold to sqrt[Δ]. The FRT density further exhibits exponential behavior over medium to long timescales, and mixtures of exponentials at shorter FRTs, thereby establishing the absence of strict power-law scaling in these threshold-crossing statistics.

A General, Noise-Driven Mechanism for the 1/f-Like Behavior of Neural Field Spectra

Consistent observations across recording modalities, experiments, and neural systems find neural field spectra with 1/f-like scaling, eliciting many alternative theories to explain this universal phenomenon. We show that a general dynamical system with stochastic drive and minimal assumptions generates 1/f-like spectra consistent with the range of values observed in vivo without requiring a specific biological mechanism or collective critical behavior.

Unconventional criticality, scaling breakdown, and diverse universality classes in the Wilson-Cowan model of neural dynamics

The Wilson-Cowan model constitutes a paradigmatic approach to understanding the collective dynamics of networks of excitatory and inhibitory units. It has been profusely used in the literature to analyze the possible phases of neural networks at a mean-field level, e.g., assuming large fully connected networks. Moreover, its stochastic counterpart allows one to study fluctuation-induced phenomena, such as avalanches. Here we revisit the stochastic Wilson-Cowan model paying special attention to the possible phase transitions between quiescent and active phases. We unveil eight possible types of such transitions, including continuous ones with scaling behavior belonging to known universality classes-such as directed percolation and tricritical directed percolation-as well as six distinct ones. In particular, we show that under some special circumstances, at a so-called "Hopf tricritical directed percolation" transition, rather unconventional behavior is observed, including the emergence of scaling breakdown. Other transitions are discontinuous and show different types of anomalies in scaling and/or exhibit mixed features of continuous and discontinuous transitions. These results broaden our knowledge of the possible types of critical behavior in networks of excitatory and inhibitory units and are, thus, of relevance to understanding avalanche dynamics in actual neuronal recordings. From a more general perspective, these results help extend the theory of nonequilibrium phase transitions into quiescent or absorbing states.

Less is different: Why sparse networks with inhibition differ from complete graphs

In neuronal systems, inhibition contributes to stabilizing dynamics and regulating pattern formation. Through developing mean-field theories of neuronal models, using complete graph networks, inhibition is commonly viewed as one "control parameter" of the system, promoting an absorbing phase transition. Here, we show that, for low connectivity sparse networks, inhibition weight is not a control parameter of the absorbing transition. We present analytical and simulation results using generic stochastic integrate-and-fire neurons that, under specific restrictions, become other simpler stochastic neuron models common in literature, which allows us to show that our results are valid for those models as well. We also give a simple explanation about why the inhibition role depends on topology, even when the topology has a dimensionality greater than the critical one. The absorbing transition independence of the inhibitory weight may be an important feature of a sparse network, as it will allow the network to maintain a near-critical regime, self-tuning average excitation, but at the same time have the freedom to adjust inhibitory weights for computation, learning, and memory, exploiting the benefits of criticality.

Fluctuation-dissipation relations in the imbalanced Wilson-Cowan model

The relation between spontaneous and stimulated brain activity is a fundamental question in neuroscience which has received wide attention in experimental studies. Recently, it has been suggested that the evoked response to external stimuli can be predicted from temporal correlations of spontaneous activity. Previous theoretical results, confirmed by the comparison with magnetoencephalography data for human brains, were obtained for the Wilson-Cowan model in the condition of balance of excitation and inhibition, a signature of a healthy brain. Here we extend previous studies to imbalanced conditions by examining a region of parameter space around the balanced fixed point. Analytical results are compared to numerical simulations of Wilson-Cowan networks. We evidence that in imbalanced conditions the functional form of the time correlation and response functions can show several behaviors, exhibiting also an oscillating regime caused by the emergence of complex eigenvalues. The analytical predictions are fully in agreement with numerical simulations, validating the role of cross-correlations in the response function. Furthermore, we identify the leading role of inhibitory neurons in controlling the overall activity of the system, tuning the level of excitability and imbalance.

Bistability and criticality in the stochastic Wilson-Cowan model

We study a stochastic version of the Wilson-Cowan model of neural dynamics, where the response function of neurons grows faster than linearly above the threshold. The model shows a region of parameters where two attractive fixed points of the dynamics exist simultaneously. One fixed point is characterized by lower activity and scale-free critical behavior, while the second fixed point corresponds to a higher (supercritical) persistent activity, with small fluctuations around a mean value. When the number of neurons is not too large, the system can switch between these two different states with a probability depending on the parameters of the network. Along with alternation of states, the model displays a bimodal distribution of the avalanches of activity, with a power-law behavior corresponding to the critical state, and a bump of very large avalanches due to the high-activity supercritical state. The bistability is due to the presence of a first-order (discontinuous) transition in the phase diagram, and the observed critical behavior is connected with the line where the low-activity state becomes unstable (spinodal line).

Power spectrum and critical exponents in the 2D stochastic Wilson-Cowan model

The power spectrum of brain activity is composed by peaks at characteristic frequencies superimposed to a background that decays as a power law of the frequency, [Formula: see text], with an exponent [Formula: see text] close to 1 (pink noise). This exponent is predicted to be connected with the exponent [Formula: see text] related to the scaling of the average size with the duration of avalanches of activity. "Mean field" models of neural dynamics predict exponents [Formula: see text] and [Formula: see text] equal or near 2 at criticality (brown noise), including the simple branching model and the fully-connected stochastic Wilson-Cowan model. We here show that a 2D version of the stochastic Wilson-Cowan model, where neuron connections decay exponentially with the distance, is characterized by exponents [Formula: see text] and [Formula: see text] markedly different from those of mean field, respectively around 1 and 1.3. The exponents [Formula: see text] and [Formula: see text] of avalanche size and duration distributions, equal to 1.5 and 2 in mean field, decrease respectively to [Formula: see text] and [Formula: see text]. This seems to suggest the possibility of a different universality class for the model in finite dimension.

Self-organized criticality in a mesoscopic model of excitatory-inhibitory neuronal populations by short-term and long-term synaptic plasticity

Dynamics of an interconnected population of excitatory and inhibitory spiking neurons wandering around a Bogdanov-Takens (BT) bifurcation point can generate the observed scale-free avalanches at the population level and the highly variable spike patterns of individual neurons. These characteristics match experimental findings for spontaneous intrinsic activity in the brain. In this paper, we address the mechanisms causing the system to get and remain near this BT point. We propose an effective stochastic neural field model which captures the dynamics of the mean-field model. We show how the network tunes itself through local long-term synaptic plasticity by STDP and short-term synaptic depression to be close to this bifurcation point. The mesoscopic model that we derive matches the directed percolation model at the absorbing state phase transition.

Addressing skepticism of the critical brain hypothesis

The hypothesis that living neural networks operate near a critical phase transition point has received substantial discussion. This “criticality hypothesis” is potentially important because experiments and theory show that optimal information processing and health are associated with operating near the critical point. Despite the promise of this idea, there have been several objections to it. While earlier objections have been addressed already, the more recent critiques of Touboul and Destexhe have not yet been fully met. The purpose of this paper is to describe their objections and offer responses. Their first objection is that the well-known Brunel model for cortical networks does not display a peak in mutual information near its phase transition, in apparent contradiction to the criticality hypothesis. In response I show that it does have such a peak near the phase transition point, provided it is not strongly driven by random inputs. Their second objection is that even simple models like a coin flip can satisfy multiple criteria of criticality. This suggests that the emergent criticality claimed to exist in cortical networks is just the consequence of a random walk put through a threshold. In response I show that while such processes can produce many signatures criticality, these signatures (1) do not emerge from collective interactions, (2) do not support information processing, and (3) do not have long-range temporal correlations. Because experiments show these three features are consistently present in living neural networks, such random walk models are inadequate. Nevertheless, I conclude that these objections have been valuable for refining research questions and should always be welcomed as a part of the scientific process.

Scaling of avalanche shape and activity power spectrum in neuronal networks

Many systems in nature exhibit avalanche dynamics with scale-free features. A general scaling theory has been proposed for critical avalanche profiles in crackling noise, predicting the collapse onto a universal avalanche shape, as well as the scaling behavior of the activity power spectrum as Brown noise. Recently, much attention has been given to the profile of neuronal avalanches, measured in neuronal systems in vitro and in vivo. Although a universal profile was evidenced, confirming the validity of the general scaling theory, the parallel study of the power spectrum scaling under the same conditions was not performed. The puzzling observation is that in the majority of healthy neuronal systems the power spectrum exhibits a behavior close to 1/f, rather than Brown, noise. Here we perform a numerical study of the scaling behavior of the avalanche shape and the power spectrum for a model of integrate and fire neurons with a short-term plasticity parameter able to tune the system to criticality. We confirm that, at criticality, the average avalanche size and the avalanche profile fulfill the general avalanche scaling theory. However, the power spectrum consistently exhibits Brown noise behavior, for both fully excitatory networks and systems with 30% inhibitory networks. Conversely, a behavior closer to 1/f noise is observed in systems slightly off criticality. Results suggest that the power spectrum is a good indicator to determine how close neuronal activity is to criticality.

Disentangling the critical signatures of neural activity

The critical brain hypothesis has emerged as an attractive framework to understand neuronal activity, but it is still widely debated. In this work, we analyze data from a multi-electrodes array in the rat's cortex and we find that power-law neuronal avalanches satisfying the crackling-noise relation coexist with spatial correlations that display typical features of critical systems. In order to shed a light on the underlying mechanisms at the origin of these signatures of criticality, we introduce a paradigmatic framework with a common stochastic modulation and pairwise linear interactions inferred from our data. We show that in such models power-law avalanches that satisfy the crackling-noise relation emerge as a consequence of the extrinsic modulation, whereas scale-free correlations are solely determined by internal interactions. Moreover, this disentangling is fully captured by the mutual information in the system. Finally, we show that analogous power-law avalanches are found in more realistic models of neural activity as well, suggesting that extrinsic modulation might be a broad mechanism for their generation.