Universal critical dynamics in high resolution neuronal avalanche data

Critical and Ictal Phases in Simulated EEG Signals on a Small-World Network

Healthy brain function is marked by neuronal network dynamics at or near the critical phase, which separates regimes of instability and stasis. A failure to remain at this critical point can lead to neurological disorders such as epilepsy, which is associated with pathological synchronization of neuronal oscillations. Using full Hodgkin-Huxley (HH) simulations on a Small-World Network, we are able to generate synthetic electroencephalogram (EEG) signals with intervals corresponding to seizure (ictal) or non-seizure (interictal) states that can occur based on the hyperexcitability of the artificial neurons and the strength and topology of the synaptic connections between them. These interictal simulations can be further classified into scale-free critical phases and disjoint subcritical exponential phases. By changing the HH parameters, we can model seizures due to a variety of causes, including traumatic brain injury (TBI), congenital channelopathies, and idiopathic etiologies, as well as the effects of anticonvulsant drugs. The results of this work may be used to help identify parameters from actual patient EEG or electrocorticographic (ECoG) data associated with ictogenesis, as well as generating simulated data for training machine-learning seizure prediction algorithms.

Differential effects of propofol and ketamine on critical brain dynamics

Whether the brain operates at a critical "tipping" point is a long standing scientific question, with evidence from both cellular and systems-scale studies suggesting that the brain does sit in, or near, a critical regime. Neuroimaging studies of humans in altered states of consciousness have prompted the suggestion that maintenance of critical dynamics is necessary for the emergence of consciousness and complex cognition, and that reduced or disorganized consciousness may be associated with deviations from criticality. Unfortunately, many of the cellular-level studies reporting signs of criticality were performed in non-conscious systems (in vitro neuronal cultures) or unconscious animals (e.g. anaesthetized rats). Here we attempted to address this knowledge gap by exploring critical brain dynamics in invasive ECoG recordings from multiple sessions with a single macaque as the animal transitioned from consciousness to unconsciousness under different anaesthetics (ketamine and propofol). We use a previously-validated test of criticality: avalanche dynamics to assess the differences in brain dynamics between normal consciousness and both drug-states. Propofol and ketamine were selected due to their differential effects on consciousness (ketamine, but not propofol, is known to induce an unusual state known as "dissociative anaesthesia"). Our analyses indicate that propofol dramatically restricted the size and duration of avalanches, while ketamine allowed for more awake-like dynamics to persist. In addition, propofol, but not ketamine, triggered a large reduction in the complexity of brain dynamics. All states, however, showed some signs of persistent criticality when testing for exponent relations and universal shape-collapse. Further, maintenance of critical brain dynamics may be important for regulation and control of conscious awareness.

Hopf Bifurcation in Mean Field Explains Critical Avalanches in Excitation-Inhibition Balanced Neuronal Networks: A Mechanism for Multiscale Variability

Cortical neural circuits display highly irregular spiking in individual neurons but variably sized collective firing, oscillations and critical avalanches at the population level, all of which have functional importance for information processing. Theoretically, the balance of excitation and inhibition inputs is thought to account for spiking irregularity and critical avalanches may originate from an underlying phase transition. However, the theoretical reconciliation of these multilevel dynamic aspects in neural circuits remains an open question. Herein, we study excitation-inhibition (E-I) balanced neuronal network with biologically realistic synaptic kinetics. It can maintain irregular spiking dynamics with different levels of synchrony and critical avalanches emerge near the synchronous transition point. We propose a novel semi-analytical mean-field theory to derive the field equations governing the network macroscopic dynamics. It reveals that the E-I balanced state of the network manifesting irregular individual spiking is characterized by a macroscopic stable state, which can be either a fixed point or a periodic motion and the transition is predicted by a Hopf bifurcation in the macroscopic field. Furthermore, by analyzing public data, we find the coexistence of irregular spiking and critical avalanches in the spontaneous spiking activities of mouse cortical slice in vitro, indicating the universality of the observed phenomena. Our theory unveils the mechanism that permits complex neural activities in different spatiotemporal scales to coexist and elucidates a possible origin of the criticality of neural systems. It also provides a novel tool for analyzing the macroscopic dynamics of E-I balanced networks and its relationship to the microscopic counterparts, which can be useful for large-scale modeling and computation of cortical dynamics.

Scaling theory of armed-conflict avalanches

Armed conflict data display features consistent with scaling and universal dynamics in both social and physical properties like fatalities and geographic extent. We propose a randomly branching armed conflict model to relate the multiple properties to one another. The model incorporates a fractal lattice on which conflict spreads, uniform dynamics driving conflict growth, and regional virulence that modulates local conflict intensity. The quantitative constraints on scaling and universal dynamics we use to develop our minimal model serve more generally as a set of constraints for other models for armed conflict dynamics. We show how this approach akin to thermodynamics imparts mechanistic intuition and unifies multiple conflict properties, giving insight into causation, prediction, and intervention timing.

Network structure of cascading neural systems predicts stimulus propagation and recovery

OBJECTIVE

Many neural systems display spontaneous, spatiotemporal patterns of neural activity that are crucial for information processing. While these cascading patterns presumably arise from the underlying network of synaptic connections between neurons, the precise contribution of the network's local and global connectivity to these patterns and information processing remains largely unknown.

APPROACH

Here, we demonstrate how network structure supports information processing through network dynamics in empirical and simulated spiking neurons using mathematical tools from linear systems theory, network control theory, and information theory.

MAIN RESULTS

In particular, we show that activity, and the information that it contains, travels through cycles in real and simulated networks.

SIGNIFICANCE

Broadly, our results demonstrate how cascading neural networks could contribute to cognitive faculties that require lasting activation of neuronal patterns, such as working memory or attention.

Spike-phase coupling as an order parameter in a leaky integrate-and-fire model

It is known that the leaky integrate-and-fire neural model shows a transition from irregular to synchronous firing by increasing the coupling between the neurons. However, a quantitative characterization of this order-disorder transition, that is, the determination of the order of transition and also the critical exponents in the case of continuous transition, is not entirely known. In this work, we consider a network of N excitatory neurons with local connections, residing on a square lattice with periodic boundary conditions. The cooperation between neurons K plays the role of the control parameter that generates criticality when the critical cooperation strength K_{c} is adopted. We introduce the population-averaged voltage (PAV) as a representative value of the network's cooperative activity. Then, we show that the coupling between the timing of spikes and the phase of temporal fluctuations of PAV defined as m resorts to identify a Kuramoto order parameter. By increasing K, we find a continuous transition from irregular spiking to a phase-locked state at the critical point K_{c}. We deploy the finite-size scaling analysis to calculate the critical exponents of this transition. To explore the formal indicator of criticality, we study the neuronal avalanches profile at this critical point and find a scaling behavior with the exponents in a fair agreement with the experimental values both in vivo and in vitro.

Impact of Physical Obstacles on the Structural and Effective Connectivity of in silico Neuronal Circuits

Scaffolds and patterned substrates are among the most successful strategies to dictate the connectivity between neurons in culture. Here, we used numerical simulations to investigate the capacity of physical obstacles placed on a flat substrate to shape structural connectivity, and in turn collective dynamics and effective connectivity, in biologically-realistic neuronal networks. We considered μ-sized obstacles placed in mm-sized networks. Three main obstacle shapes were explored, namely crosses, circles and triangles of isosceles profile. They occupied either a small area fraction of the substrate or populated it entirely in a periodic manner. From the point of view of structure, all obstacles promoted short length-scale connections, shifted the in- and out-degree distributions toward lower values, and increased the modularity of the networks. The capacity of obstacles to shape distinct structural traits depended on their density and the ratio between axonal length and substrate diameter. For high densities, different features were triggered depending on obstacle shape, with crosses trapping axons in their vicinity and triangles funneling axons along the reverse direction of their tip. From the point of view of dynamics, obstacles reduced the capacity of networks to spontaneously activate, with triangles in turn strongly dictating the direction of activity propagation. Effective connectivity networks, inferred using transfer entropy, exhibited distinct modular traits, indicating that the presence of obstacles facilitated the formation of local effective microcircuits. Our study illustrates the potential of physical constraints to shape structural blueprints and remodel collective activity, and may guide investigations aimed at mimicking organizational traits of biological neuronal circuits.

Between Perfectly Critical and Fully Irregular: A Reverberating Model Captures and Predicts Cortical Spike Propagation

Knowledge about the collective dynamics of cortical spiking is very informative about the underlying coding principles. However, even most basic properties are not known with certainty, because their assessment is hampered by spatial subsampling, i.e., the limitation that only a tiny fraction of all neurons can be recorded simultaneously with millisecond precision. Building on a novel, subsampling-invariant estimator, we fit and carefully validate a minimal model for cortical spike propagation. The model interpolates between two prominent states: asynchronous and critical. We find neither of them in cortical spike recordings across various species, but instead identify a narrow "reverberating" regime. This approach enables us to predict yet unknown properties from very short recordings and for every circuit individually, including responses to minimal perturbations, intrinsic network timescales, and the strength of external input compared to recurrent activation "thereby informing about the underlying coding principles for each circuit, area, state and task.

Avalanche properties in striplike ferromagnetic systems

We present numerical findings on the behavior of the athermal nonequilibrium random-field Ising model of spins at the thin striplike L_{1}×L_{2}×L_{3} cubic lattices with L_{1}<L_{2}<L_{3}. Changing of system sizes highly influences the evolution and shape of avalanches. The smallest avalanches [classified as three-dimension- (3D) like] are unaffected by the system boundaries, the larger are sandwiched between the top and bottom system faces so are 2D-like, while the largest are extended over the system lateral cross section and propagate along the length L_{3} like in 1D systems. Such a structure of avalanches causes double power-law distributions of their size, duration, and energy with larger effective critical exponent corresponding to 3D-like and smaller to 2D-like avalanches. The distributions scale with thickness L_{1} and are collapsible following the proposed scaling predictions which, together with the distributions' shape, might be important for analysis of the Barkhausen noise experimental data for striplike samples. Finally, the impact of system size on external field that triggers the largest avalanche for a given disorder is presented and discussed.

Stability of motor cortex network states during learning-associated neural reorganizations

A substantial reorganization of neural activity and neuron-to-movement relationship in motor cortical circuits accompanies the emergence of reproducible movement patterns during motor learning. Little is known about how this tempest of neural activity restructuring impacts the stability of network states in recurrent cortical circuits. To investigate this issue, we reanalyzed data in which we recorded for 14 days via population calcium imaging the activity of the same neural populations of pyramidal neurons in layer 2/3 and layer 5 of forelimb motor and premotor cortex in mice during the daily learning of a lever-press task. We found that motor cortex network states remained stable with respect to the critical network state during the extensive reorganization of both neural population activity and its relation to lever movement throughout learning. Specifically, layer 2/3 cortical circuits unceasingly displayed robust evidence for operating at the critical network state, a regime that maximizes information capacity and transmission and provides a balance between network robustness and flexibility. In contrast, layer 5 circuits operated away from the critical network state for all 14 days of recording and learning. In conclusion, this result indicates that the wide-ranging malleability of synapses, neurons, and neural connectivity during learning operates within the constraint of a stable and layer-specific network state regarding dynamic criticality, and suggests that different cortical layers operate under distinct constraints because of their specialized goals.NEW & NOTEWORTHY The neural activity reorganizes throughout motor learning, but how this reorganization impacts the stability of network states is unclear. We used two-photon calcium imaging to investigate how the network states in layer 2/3 and layer 5 of forelimb motor and premotor cortex are modulated by motor learning. We show that motor cortex network states are layer-specific and constant regarding criticality during neural activity reorganization, and suggests that layer-specific constraints could be motivated by different functions.

Why Brain Criticality Is Clinically Relevant: A Scoping Review

The past 25 years have seen a strong increase in the number of publications related to criticality in different areas of neuroscience. The potential of criticality to explain various brain properties, including optimal information processing, has made it an increasingly exciting area of investigation for neuroscientists. Recent reviews on this topic, sometimes termed brain criticality, make brief mention of clinical applications of these findings to several neurological disorders such as epilepsy, neurodegenerative disease, and neonatal hypoxia. Other clinicallyrelevant domains - including anesthesia, sleep medicine, developmental-behavioral pediatrics, and psychiatry - are seldom discussed in review papers of brain criticality. Thorough assessments of these application areas and their relevance for clinicians have also yet to be published. In this scoping review, studies of brain criticality involving human data of all ages are evaluated for their current and future clinical relevance. To make the results of these studies understandable to a more clinical audience, a review of the key concepts behind criticality (e.g., phase transitions, long-range temporal correlation, self-organized criticality, power laws, branching processes) precedes the discussion of human clinical studies. Open questions and forthcoming areas of investigation are also considered.

Time-dependent branching processes: a model of oscillating neuronal avalanches

Recently, neuronal avalanches have been observed to display oscillations, a phenomenon regarded as the co-existence of a scale-free behaviour (the avalanches close to criticality) and scale-dependent dynamics (the oscillations). Ordinary continuous-time branching processes with constant extinction and branching rates are commonly used as models of neuronal activity, yet they lack any such time-dependence. In the present work, we extend a basic branching process by allowing the extinction rate to oscillate in time as a new model to describe cortical dynamics. By means of a perturbative field theory, we derive relevant observables in closed form. We support our findings by quantitative comparison to numerics and qualitative comparison to available experimental results.

Firing rate homeostasis counteracts changes in stability of recurrent neural networks caused by synapse loss in Alzheimer's disease

The impairment of cognitive function in Alzheimer's disease is clearly correlated to synapse loss. However, the mechanisms underlying this correlation are only poorly understood. Here, we investigate how the loss of excitatory synapses in sparsely connected random networks of spiking excitatory and inhibitory neurons alters their dynamical characteristics. Beyond the effects on the activity statistics, we find that the loss of excitatory synapses on excitatory neurons reduces the network's sensitivity to small perturbations. This decrease in sensitivity can be considered as an indication of a reduction of computational capacity. A full recovery of the network's dynamical characteristics and sensitivity can be achieved by firing rate homeostasis, here implemented by an up-scaling of the remaining excitatory-excitatory synapses. Mean-field analysis reveals that the stability of the linearised network dynamics is, in good approximation, uniquely determined by the firing rate, and thereby explains why firing rate homeostasis preserves not only the firing rate but also the network's sensitivity to small perturbations.

Temporal Signatures of Criticality in Human Cortical Excitability as Probed by Early Somatosensory Responses

Brain responses vary considerably from moment to moment, even to identical sensory stimuli. This has been attributed to changes in instantaneous neuronal states determining the system's excitability. Yet the spatiotemporal organization of these dynamics remains poorly understood. Here we test whether variability in stimulus-evoked activity can be interpreted within the framework of criticality, which postulates dynamics of neural systems to be tuned toward the phase transition between stability and instability as is reflected in scale-free fluctuations in spontaneous neural activity. Using a novel noninvasive approach in 33 male human participants, we tracked instantaneous cortical excitability by inferring the magnitude of excitatory postsynaptic currents from the N20 component of the somatosensory evoked potential. Fluctuations of cortical excitability demonstrated long-range temporal dependencies decaying according to a power law across trials, a hallmark of systems at critical states. As these dynamics covaried with changes in prestimulus oscillatory activity in the alpha band (8-13 Hz), we establish a mechanistic link between ongoing and evoked activity through cortical excitability and argue that the co-emergence of common temporal power laws may indeed originate from neural networks poised close to a critical state. In contrast, no signatures of criticality were found in subcortical or peripheral nerve activity. Thus, criticality may represent a parsimonious organizing principle of variability in stimulus-related brain processes on a cortical level, possibly reflecting a delicate equilibrium between robustness and flexibility of neural responses to external stimuli.SIGNIFICANCE STATEMENT Variability of neural responses in primary sensory areas is puzzling, as it is detrimental to the exact mapping between stimulus features and neural activity. However, such variability can be beneficial for information processing in neural networks if it is of a specific nature, namely, if dynamics are poised at a so-called critical state characterized by a scale-free spatiotemporal structure. Here, we demonstrate the existence of a link between signatures of criticality in ongoing and evoked activity through cortical excitability, which fills the long-standing gap between two major directions of research on neural variability: the impact of instantaneous brain states on stimulus processing on the one hand and the scale-free organization of spatiotemporal network dynamics of spontaneous activity on the other.

Self-regulated critical brain dynamics originate from high frequency-band activity in the MEG

Brain function requires the flexible coordination of billions of neurons across multiple scales. This could be achieved by scale-free, critical dynamics balanced at the edge of order and disorder. Criticality has been demonstrated in several, often reduced neurophysiological model systems. In the intact human brain criticality has yet been only verified for the resting state. A more direct link between the concept of criticality and oscillatory brain physiology, which is strongly related to cognition, is yet missing. In the present study we therefore carried out a frequency-specific analysis of criticality in the MEG, recorded while subjects were in a defined cognitive state through mindfulness meditation. In a two-step approach we assessed whether the macroscopic neural avalanche dynamics is scale-free by evaluating the goodness of a power-law fits of cascade size and duration distributions of MEG deflections in different frequency bands. In a second step we determined the closeness of the power-law exponents to a critical value of -1.5. Power-law fitting was evaluated by permutation testing, fitting of alternative distributions, and cascade shape analysis. Criticality was verified by defined relationships of exponents of cascade size and duration distributions. Behavioral relevance of criticality was tested by correlation of indices of criticality with individual scores of the Mindful Attention Awareness Scale. We found that relevant scale-free near-critical dynamics originated only from broad-band high-frequency (> 100 Hz) MEG activity, which has been associated with action potential firing, and therefore links criticality on the macroscopic level of MEG to critical spike avalanches on a microscopic level. Whereas a scale-free dynamics was found under mindfulness meditation and rest, avalanche dynamics shifted towards a critical point during meditation by reduction of neural noise. Together with our finding that during mindfulness meditation avalanches show differences in topography relative to rest, our results show that self-regulated attention as required during meditation can serve as a control parameter of criticality in scale-free brain dynamics.

Sensitivity of collective outcomes identifies pivotal components

A social system is susceptible to perturbation when its collective properties depend sensitively on a few pivotal components. Using the information geometry of minimal models from statistical physics, we develop an approach to identify pivotal components to which coarse-grained, or aggregate, properties are sensitive. As an example, we introduce our approach on a reduced toy model with a median voter who always votes in the majority. The sensitivity of majority-minority divisions to changing voter behaviour pinpoints the unique role of the median. More generally, the sensitivity identifies pivotal components that precisely determine collective outcomes generated by a complex network of interactions. Using perturbations to target pivotal components in the models, we analyse datasets from political voting, finance and Twitter. Across these systems, we find remarkable variety, from systems dominated by a median-like component to those whose components behave more equally. In the context of political institutions such as courts or legislatures, our methodology can help describe how changes in voters map to new collective voting outcomes. For economic indices, differing system response reflects varying fiscal conditions across time. Thus, our information-geometric approach provides a principled, quantitative framework that may help assess the robustness of collective outcomes to targeted perturbation and compare social institutions, or even biological networks, with one another and across time.

Atomic Scale Dynamics Drive Brain-like Avalanches in Percolating Nanostructured Networks

Self-assembled networks of nanoparticles and nanowires have recently emerged as promising systems for brain-like computation. Here, we focus on percolating networks of nanoparticles which exhibit brain-like dynamics. We use a combination of experiments and simulations to show that the brain-like network dynamics emerge from atomic-scale switching dynamics inside tunnel gaps that are distributed throughout the network. The atomic-scale dynamics emulate leaky integrate and fire (LIF) mechanisms in biological neurons, leading to the generation of critical avalanches of signals. These avalanches are quantitatively the same as those observed in cortical tissue and are signatures of the correlations that are required for computation. We show that the avalanches are associated with dynamical restructuring of the networks which self-tune to balanced states consistent with self-organized criticality. Our simulations allow visualization of the network states and detailed mechanisms of signal propagation.

Conservation laws by virtue of scale symmetries in neural systems

In contrast to the symmetries of translation in space, rotation in space, and translation in time, the known laws of physics are not universally invariant under transformation of scale. However, a special case exists in which the action is scale invariant if it satisfies the following two constraints: 1) it must depend upon a scale-free Lagrangian, and 2) the Lagrangian must change under scale in the same way as the inverse time, [Formula: see text]. Our contribution lies in the derivation of a generalised Lagrangian, in the form of a power series expansion, that satisfies these constraints. This generalised Lagrangian furnishes a normal form for dynamic causal models-state space models based upon differential equations-that can be used to distinguish scale symmetry from scale freeness in empirical data. We establish face validity with an analysis of simulated data, in which we show how scale symmetry can be identified and how the associated conserved quantities can be estimated in neuronal time series.

Long-range temporal correlations in scale-free neuromorphic networks

Biological neuronal networks are the computing engines of the mammalian brain. These networks exhibit structural characteristics such as hierarchical architectures, small-world attributes, and scale-free topologies, providing the basis for the emergence of rich temporal characteristics such as scale-free dynamics and long-range temporal correlations. Devices that have both the topological and the temporal features of a neuronal network would be a significant step toward constructing a neuromorphic system that can emulate the computational ability and energy efficiency of the human brain. Here we use numerical simulations to show that percolating networks of nanoparticles exhibit structural properties that are reminiscent of biological neuronal networks, and then show experimentally that stimulation of percolating networks by an external voltage stimulus produces temporal dynamics that are self-similar, follow power-law scaling, and exhibit long-range temporal correlations. These results are expected to have important implications for the development of neuromorphic devices, especially for those based on the concept of reservoir computing.

Biological neuronal networks exhibit well-defined properties such as hierarchical structures and scale-free topologies, as well as a high degree of local clustering and short path lengths between nodes. These structural properties are intimately connected to the observed long-range temporal correlations in the network dynamics. Fabrication of artificial networks with similar structural properties would facilitate brain-like (“neuromorphic”) computing. Here we show experimentally that percolating networks of nanoparticles exhibit similar long-range temporal correlations to those of biological neuronal networks and use simulations to demonstrate that the dynamics arise from an underlying scale-free network architecture. We discuss similarities between the biological and percolating systems and highlight the potential for the percolating networks to be used in neuromorphic computing applications.

Extreme value theory of evolving phenomena in complex dynamical systems: Firing cascades in a model of a neural network

We extend the scope of the dynamical theory of extreme values to include phenomena that do not happen instantaneously but evolve over a finite, albeit unknown at the onset, time interval. We consider complex dynamical systems composed of many individual subsystems linked by a network of interactions. As a specific example of the general theory, a model of a neural network, previously introduced by other authors to describe the electrical activity of the cerebral cortex, is analyzed in detail. On the basis of this analysis, we propose a novel definition of a neuronal cascade, a physiological phenomenon of primary importance. We derive extreme value laws for the statistics of these cascades, both from the point of view of exceedances (that satisfy critical scaling theory in a certain regime) and of block maxima.

The Critically Tuned Cortex

The criticality hypothesis predicts that cortex operates near a critical point for optimum information processing. In this issue of Neuron, Ma et al. (2019) find evidence consistent with a mechanism that tunes cortex to criticality, even in the face of a strong perturbation over several days.

Stability of neuronal avalanches and long-range temporal correlations during the first year of life in human infants

During infancy, the human brain rapidly expands in size and complexity as neural networks mature and new information is incorporated at an accelerating pace. Recently, it was shown that single-electrode EEG in preterms at birth exhibits scale-invariant intermittent bursts. Yet, it is currently not known whether the normal infant brain, in particular, the cortex, maintains a distinct dynamical state during development that is characterized by scale-invariant spatial as well as temporal aspects. Here we employ dense-array EEG recordings acquired from the same infants at 6 and 12 months of age to characterize brain activity during an auditory odd-ball task. We show that suprathreshold events organize as spatiotemporal clusters whose size and duration are power-law distributed, the hallmark of neuronal avalanches. Time series of local suprathreshold EEG events display significant long-range temporal correlations (LRTCs). No differences were found between 6 and 12 months, demonstrating stability of avalanche dynamics and LRTCs during the first year after birth. These findings demonstrate that the infant brain is characterized by distinct spatiotemporal dynamical aspects that are in line with expectations of a critical cortical state. We suggest that critical state dynamics, which theory and experiments have shown to be beneficial for numerous aspects of information processing, are maintained by the infant brain to process an increasingly complex environment during development.

Description of spreading dynamics by microscopic network models and macroscopic branching processes can differ due to coalescence

Spreading processes are conventionally monitored on a macroscopic level by counting the number of incidences over time. The spreading process can then be modeled either on the microscopic level, assuming an underlying interaction network, or directly on the macroscopic level, assuming that microscopic contributions are negligible. The macroscopic characteristics of both descriptions are commonly assumed to be identical. In this work we show that these characteristics of microscopic and macroscopic descriptions can be different due to coalescence, i.e., a node being activated at the same time by multiple sources. In particular, we consider a (microscopic) branching network (probabilistic cellular automaton) with annealed connectivity disorder, record the macroscopic activity, and then approximate this activity by a (macroscopic) branching process. In this framework we analytically calculate the effect of coalescence on the collective dynamics. We show that coalescence leads to a universal nonlinear scaling function for the conditional expectation value of successive network activity. This allows us to quantify the difference between the microscopic model parameter and established estimates of the macroscopic branching parameter. To overcome this difference, we propose a nonlinear estimator that correctly infers the microscopic model parameter for all system sizes.

Emergence of power laws in noncritical neuronal systems

Experimental and computational studies provide compelling evidence that neuronal systems are characterized by power-law distributions of neuronal avalanche sizes. This fact is interpreted as an indication that these systems are operating near criticality, and, in turn, typical properties of critical dynamical processes, such as optimal information transmission and stability, are attributed to neuronal systems. The purpose of this Rapid Communication is to show that the presence of power-law distributions for the size of neuronal avalanches is not a sufficient condition for the system to operate near criticality. Specifically, we consider a simplistic model of neuronal dynamics on networks and show that the degree distribution of the underlying neuronal network may trigger power-law distributions for neuronal avalanches even when the system is not in its critical regime. To certify and explain our findings we develop an analytical approach based on percolation theory and branching processes techniques.

Through neuronal avalanches to consciousness: Conjectures and elaborations

The concept of the "conduction avalanche" as a mechanism of signal propagation in the central nervous system was introduced in 1895 by Santiago Ramón y Cajal. Through such a mechanism, stimuli received by sensory receptor cells would be augmented as they reached the cerebral cortex through the corresponding anatomical pathways, leading to conscious perception. Cajal applied the concept to the visual, auditory, olfactory and somatosensory systems, and in 1896 extended it to the cerebellar cortex. Beginning in 1899 and up to 2003, prominent neuroscientists, including Lewellys F. Barker, C. Judson Herrick, Francis X. Dercum, Rafael Lorente de Nó, Cornelius U. Ariëns Kappers, Hartwig Kuhlenbeck, Gordon M. Shepherd, and Rodolfo Llinás, have referred in their writings to the principle of conduction avalanches, crediting Cajal. In 2003, John Beggs and Dietmar Plenz introduced the concept of the "neuronal avalanche," modelled after the power law of physics, as a property of neocortical networks and a new mode of spontaneous activity, distinct from the oscillatory, synchronized and wave states previously conceived to underpin the integrative function of the cerebral cortex. The topic has become an issue of intensive research over the past 15 years. In this paper, we discuss the basic tenets of the Cajalian principle, followed by an exposition of ideas throughout the 20th century, and an overview, in a modern perspective, of the neuronal avalanche as a mechanism in the current study of the neural bases of consciousness.

Critical synchronization dynamics of the Kuramoto model on connectome and small world graphs

The hypothesis, that cortical dynamics operates near criticality also suggests, that it exhibits universal critical exponents which marks the Kuramoto equation, a fundamental model for synchronization, as a prime candidate for an underlying universal model. Here, we determined the synchronization behavior of this model by solving it numerically on a large, weighted human connectome network, containing 836733 nodes, in an assumed homeostatic state. Since this graph has a topological dimension d < 4, a real synchronization phase transition is not possible in the thermodynamic limit, still we could locate a transition between partially synchronized and desynchronized states. At this crossover point we observe power-law–tailed synchronization durations, with τ_t ≃ 1.2(1), away from experimental values for the brain. For comparison, on a large two-dimensional lattice, having additional random, long-range links, we obtain a mean-field value: τ_t ≃ 1.6(1). However, below the transition of the connectome we found global coupling control-parameter dependent exponents 1 < τ_t ≤ 2, overlapping with the range of human brain experiments. We also studied the effects of random flipping of a small portion of link weights, mimicking a network with inhibitory interactions, and found similar results. The control-parameter dependent exponent suggests extended dynamical criticality below the transition point.

Spike-Timing-Dependent Plasticity With Axonal Delay Tunes Networks of Izhikevich Neurons to the Edge of Synchronization Transition With Scale-Free Avalanches

Critical brain hypothesis has been intensively studied both in experimental and theoretical neuroscience over the past two decades. However, some important questions still remain: (i) What is the critical point the brain operates at? (ii) What is the regulatory mechanism that brings about and maintains such a critical state? (iii) The critical state is characterized by scale-invariant behavior which is seemingly at odds with definitive brain oscillations? In this work we consider a biologically motivated model of Izhikevich neuronal network with chemical synapses interacting via spike-timing-dependent plasticity (STDP) as well as axonal time delay. Under generic and physiologically relevant conditions we show that the system is organized and maintained around a synchronization transition point as opposed to an activity transition point associated with an absorbing state phase transition. However, such a state exhibits experimentally relevant signs of critical dynamics including scale-free avalanches with finite-size scaling as well as critical branching ratios. While the system displays stochastic oscillations with highly correlated fluctuations, it also displays dominant frequency modes seen as sharp peaks in the power spectrum. The role of STDP as well as time delay is crucial in achieving and maintaining such critical dynamics, while the role of inhibition is not as crucial. In this way we provide possible answers to all three questions posed above. We also show that one can achieve supercritical or subcritical dynamics if one changes the average time delay associated with axonal conduction.

Efficient message passing for cascade size distributions

How big is the risk that a few initial failures of networked nodes amplify to large cascades that endanger the functioning of the system? Common answers refer to the average final cascade size. Two analytic approaches allow its computation: (a) (heterogeneous) mean field approximation and (b) belief propagation. The former applies to (infinitely) large locally tree-like networks, while the latter is exact on finite trees. Yet, cascade sizes can have broad and multi-modal distributions that are not well represented by their average. Full distribution information is essential to identify likely events and to estimate the tail risk, i.e. the probability of extreme events. We therefore present an efficient message passing algorithm that calculates the cascade size distribution in finite networks. It is exact on finite trees and for a large class of cascade processes. An approximate version applies to any network structure and performs well on locally tree-like networks, as we show with several examples.

The scale-invariant, temporal profile of neuronal avalanches in relation to cortical γ-oscillations

Activity cascades are found in many complex systems. In the cortex, they arise in the form of neuronal avalanches that capture ongoing and evoked neuronal activities at many spatial and temporal scales. The scale-invariant nature of avalanches suggests that the brain is in a critical state, yet predictions from critical theory on the temporal unfolding of avalanches have yet to be confirmed in vivo. Here we show in awake nonhuman primates that the temporal profile of avalanches follows a symmetrical, inverted parabola spanning up to hundreds of milliseconds. This parabola constrains how avalanches initiate locally, extend spatially and shrink as they evolve in time. Importantly, parabolas of different durations can be collapsed with a scaling exponent close to 2 supporting critical generational models of neuronal avalanches. Spontaneously emerging, transient γ-oscillations coexist with and modulate these avalanche parabolas thereby providing a temporal segmentation to inherently scale-invariant, critical dynamics. Our results identify avalanches and oscillations as dual principles in the temporal organization of brain activity.

Avalanches in Compressed Sandstone: Crackling Noise under Confinement

The acoustic emission, AE, from avalanches of local cracks and microstructural changes of sandstone under confined compression have been reported. These avalanches soften the underlying minerals and play a key role as indicators for the prediction of geo-engineering disasters, such as mining collapses, rock outbursts caused by high ground stress, and man-made quakes by fracking. Compressed sandstone is a model material for the investigation of avalanches. The avalanche energies, amplitudes, and waiting times show the probability distributions that allow us to distinguish between three compression stages; namely, (I) pre-failure, (II) correlated failure, and (III) post-failure. The energy of stage I and stage II is power-law distributed and scale invariant, while post-failure experiments show power laws with high exponential damping (friction). The scaling behavior is close to the predictions of a mean-field (MF) model (stage II) and a force-integrated mean-field model (stage I). Confinement shifts the value of the energy exponent closer to the MF prediction. Omori’s law and waiting time distributions are independent of stress during the compression; their scaling exponents are very similar to those found in seismological studies.

Avalanches and criticality in self-organized nanoscale networks

Current efforts to achieve neuromorphic computation are focused on highly organized architectures, such as integrated circuits and regular arrays of memristors, which lack the complex interconnectivity of the brain and so are unable to exhibit brain-like dynamics. New architectures are required, both to emulate the complexity of the brain and to achieve critical dynamics and consequent maximal computational performance. We show here that electrical signals from self-organized networks of nanoparticles exhibit brain-like spatiotemporal correlations and criticality when fabricated at a percolating phase transition. Specifically, the sizes and durations of avalanches of switching events are power law distributed, and the power law exponents satisfy rigorous criteria for criticality. These signals are therefore qualitatively and quantitatively similar to those measured in the cortex. Our self-organized networks provide a low-cost platform for computational approaches that rely on spatiotemporal correlations, such as reservoir computing, and are an important step toward creating neuromorphic device architectures.

Critical Dynamics Mediate Learning of New Distributed Memory Representations in Neuronal Networks

We explore the possible role of network dynamics near a critical point in the storage of new information in silico and in vivo, and show that learning and memory may rely on neuronal network features mediated by the vicinity of criticality. Using a mean-field, attractor-based model, we show that new information can be consolidated into attractors through state-based learning in a dynamical regime associated with maximal susceptibility at the critical point. Then, we predict that the subsequent consolidation process results in a shift from critical to sub-critical dynamics to fully encapsulate the new information. We go on to corroborate these findings using analysis of rodent hippocampal CA1 activity during contextual fear memory (CFM) consolidation. We show that the dynamical state of the CA1 network is inherently poised near criticality, but the network also undergoes a shift towards sub-critical dynamics due to successful consolidation of the CFM. Based on these findings, we propose that dynamical features associated with criticality may be universally necessary for storing new memories.

Coarse Graining, Fixed Points, and Scaling in a Large Population of Neurons

We develop a phenomenological coarse-graining procedure for activity in a large network of neurons, and apply this to recordings from a population of 1000+ cells in the hippocampus. Distributions of coarse-grained variables seem to approach a fixed non-Gaussian form, and we see evidence of scaling in both static and dynamic quantities. These results suggest that the collective behavior of the network is described by a nontrivial fixed point.

Critical disorder and critical magnetic field of the nonequilibrium athermal random-field Ising model in thin systems

In the present study of the nonequilibrium athermal random-field Ising model we focus on the behavior of the critical disorder R_{c}(l) and the critical magnetic field H_{c}(l) under different boundary conditions when the system thickness l varies. We propose expressions for R_{c}(l) and H_{c}(l) as well as for the effective critical disorder R_{c}^{eff}(l,L) and effective critical magnetic field H_{c}^{eff}(l,L) playing the role of the effective critical parameters for the L×L×l lattices of finite lateral size L. We support these expressions by the scaling collapses of the magnetization and susceptibility curves obtained in extensive simulations. The collapses are achieved with the two-dimensional (2D) exponents for l below some characteristic value, providing thus a numerical evidence that the thin systems exhibit a 2D-like criticality which should be relevant for the experimental analyses of thin ferromagnetic samples.

Cortical Circuit Dynamics Are Homeostatically Tuned to Criticality In Vivo

Homeostatic mechanisms stabilize neuronal activity in vivo, but whether this process gives rise to balanced network dynamics is unknown. Here, we continuously monitored the statistics of network spiking in visual cortical circuits in freely behaving rats for 9 days. Under control conditions in light and dark, networks were robustly organized around criticality, a regime that maximizes information capacity and transmission. When input was perturbed by visual deprivation, network criticality was severely disrupted and subsequently restored to criticality over 48 h. Unexpectedly, the recovery of excitatory dynamics preceded homeostatic plasticity of firing rates by >30 h. We utilized model investigations to manipulate firing rate homeostasis in a cell-type-specific manner at the onset of visual deprivation. Our results suggest that criticality in excitatory networks is established by inhibitory plasticity and architecture. These data establish that criticality is consistent with a homeostatic set point for visual cortical dynamics and suggest a key role for homeostatic regulation of inhibition.

Identifying epidemic threshold by temporal profile of outbreaks on networks

Identifying epidemic threshold is of great significance in preventing and controlling disease spreading on real-world networks. Previous studies have proposed different theoretical and numerical approaches to determine the epidemic threshold for the susceptible-infected-recovered (SIR) model, but the numerical study of the critical points on networks by utilizing temporal characteristics of epidemic outbreaks is still lacking. Here, we study the temporal profile of epidemic outbreaks, i.e., the average avalanche shapes of a fixed duration. At the critical point, the rescaled average terminating and nonterminating avalanche shapes for different durations collapse onto two universal curves, respectively, while the average number of subsequent events essentially remains constant. We propose two numerical measures to determine the epidemic threshold by analyzing the convergence of the rescaled average nonterminating avalanche shapes for varying durations and the stability of the average number of subsequent events, respectively. Extensive numerical simulations demonstrate that our methods can accurately identify the numerical threshold for the SIR dynamics on synthetic and empirical networks. Compared with traditional numerical measures, our methods are more efficient due to the constriction of observation duration and thus are more applicable to large-scale networks. This work helps one to understand the temporal profile of disease propagation and would promote further studies on the phase transition of epidemic dynamics.

25 years of criticality in neuroscience - established results, open controversies, novel concepts

Twenty-five years ago, Dunkelmann and Radons (1994) showed that neural networks can self-organize to a critical state. In models, the critical state offers a number of computational advantages. Thus this hypothesis, and in particular the experimental work by Beggs and Plenz (2003), has triggered an avalanche of research, with thousands of studies referring to it. Nonetheless, experimental results are still contradictory. How is it possible, that a hypothesis has attracted active research for decades, but nonetheless remains controversial? We discuss the experimental and conceptual controversy, and then present a parsimonious solution that (i) unifies the contradictory experimental results, (ii) avoids disadvantages of a critical state, and (iii) enables rapid, adaptive tuning of network properties to task requirements.

Time-series thresholding and the definition of avalanche size

Avalanches whose sizes and durations are distributed as power laws appear in many contexts, from physics to geophysics and biology. Here we show that there is a hidden peril in thresholding continuous times series-from either empirical or synthetic data-for the identification of avalanches. In particular, we consider two possible alternative definitions of avalanche size used, e.g., in the empirical determination of avalanche exponents in the analysis of neural-activity data. By performing analytical and computational studies of an Ornstein-Uhlenbeck process (taken as a guiding example) we show that (1) if relatively large threshold values are employed to determine the beginning and ending of avalanches and (2) if-as sometimes done in the literature-avalanche sizes are defined as the total area (above zero) of the avalanche, then true asymptotic scaling behavior is not seen, instead the observations are dominated by transient effects. This problem-that we have detected in some recent works-leads to misinterpretations of the resulting scaling regimes.

Stability of neuronal avalanches and long-range temporal correlations during the first year of life in human infant

During infancy, the human brain rapidly expands in size and complexity as neural networks mature and new information is incorporated at an accelerating pace. Recently, it was shown that single electrode EEG in preterms at birth exhibits scale-invariant intermittent bursts. Yet, it is currently not known whether the normal infant brain, in particular, the cortex maintains a distinct dynamical state during development that is characterized by scale-invariant spatial as well as temporal aspects. Here we employ dense-array EEG recordings acquired from the same infants at 6 and 12 months of age to characterize brain activity during an auditory oddball task. We show that suprathreshold events organize as spatiotemporal clusters whose size and duration are power-law distributed, the hallmark of neuronal avalanches. Time series of local suprathreshold EEG events display significant long-range temporal correlations (LRTCs). No differences were found between 6 and 12 months, demonstrating stability of avalanche dynamics and LRTCs during the first year after birth. These findings demonstrate that the infant brain is characterized by distinct spatiotemporal dynamical aspects that are in line with expectations of a critical cortical state. We suggest that critical state dynamics, which theory and experiments have shown to be beneficial for numerous aspects of information processing, are maintained by the infant brain to process an increasingly complex environment during development.

Heterogeneity and Delayed Activation as Hallmarks of Self-Organization and Criticality in Excitable Tissue

Self-organized critical dynamics is assumed to be an attractive mode of functioning for several real-life systems and entails an emergent activity in which the extent of observables follows a power-law distribution. The hallmarks of criticality have recently been observed in a plethora of biological systems, including beta cell populations within pancreatic islets of Langerhans. In the present study, we systematically explored the mechanisms that drive the critical and supercritical behavior in networks of coupled beta cells under different circumstances by means of experimental and computational approaches. Experimentally, we employed high-speed functional multicellular calcium imaging of fluorescently labeled acute mouse pancreas tissue slices to record calcium signals in a large number of beta cells simultaneously, and with a high spatiotemporal resolution. Our experimental results revealed that the cellular responses to stimulation with glucose are biphasic and glucose-dependent. Under physiological as well as under supraphysiological levels of stimulation, an initial activation phase was followed by a supercritical plateau phase with a high number of global intercellular calcium waves. However, the activation phase displayed fingerprints of critical behavior under lower stimulation levels, with a progressive recruitment of cells and a power-law distribution of calcium wave sizes. On the other hand, the activation phase provoked by pathophysiologically high glucose concentrations, differed considerably and was more rapid, less continuous, and supercritical. To gain a deeper insight into the experimentally observed complex dynamical patterns, we built up a phenomenological model of coupled excitable cells and explored empirically the model’s necessities that ensured a good overlap between computational and experimental results. It turned out that such a good agreement between experimental and computational findings was attained when both heterogeneous and stimulus-dependent time lags, variability in excitability levels, as well as a heterogeneous cell-cell coupling were included into the model. Most importantly, since our phenomenological approach involved only a few parameters, it naturally lends itself not only for determining key mechanisms of self-organized criticality at the tissue level, but also points out various features for comprehensive and realistic modeling of different excitable systems in nature.

Critical excitation-inhibition balance in dense neural networks

The "edge of chaos" phase transition in artificial neural networks is of renewed interest in light of recent evidence for criticality in brain dynamics. Statistical mechanics traditionally studied this transition with connectivity k as the control parameter and an exactly balanced excitation-inhibition ratio. While critical connectivity has been found to be low in these model systems, typically around k=2, which is unrealistic for natural neural systems, a recent study utilizing the excitation-inhibition ratio as the control parameter found a new, nearly degree independent, critical point when connectivity is large. However, the new phase transition is accompanied by an unnaturally high level of activity in the network. Here we study random neural networks with the additional properties of (i) a high clustering coefficient and (ii) neurons that are solely either excitatory or inhibitory, a prominent property of natural neurons. As a result, we observe an additional critical point for networks with large connectivity, regardless of degree distribution, which exhibits low activity levels that compare well with neuronal brain networks.

Simple model of complex dynamics of activity patterns in developing networks of neuronal cultures

Living neuronal networks in dissociated neuronal cultures are widely known for their ability to generate highly robust spatiotemporal activity patterns in various experimental conditions. Such patterns are often treated as neuronal avalanches that satisfy the power scaling law and thereby exemplify self-organized criticality in living systems. A crucial question is how these patterns can be explained and modeled in a way that is biologically meaningful, mathematically tractable and yet broad enough to account for neuronal heterogeneity and complexity. Here we derive and analyse a simple network model that may constitute a response to this question. Our derivations are based on few basic phenomenological observations concerning the input-output behavior of an isolated neuron. A distinctive feature of the model is that at the simplest level of description it comprises of only two variables, the network activity variable and an exogenous variable corresponding to energy needed to sustain the activity, and few parameters such as network connectivity and efficacy of signal transmission. The efficacy of signal transmission is modulated by the phenomenological energy variable. Strikingly, this simple model is already capable of explaining emergence of network spikes and bursts in developing neuronal cultures. The model behavior and predictions are consistent with published experimental evidence on cultured neurons. At the larger, cellular automata scale, introduction of the energy-dependent regulatory mechanism results in the overall model behavior that can be characterized as balancing on the edge of the network percolation transition. Network activity in this state shows population bursts satisfying the scaling avalanche conditions. This network state is self-sustainable and represents energetic balance between global network-wide processes and spontaneous activity of individual elements.

Criticality between Cortical States

Since the first measurements of neuronal avalanches, the critical brain hypothesis has gained traction. 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, ranging from highly synchronous states (with higher spiking variability) to desynchronized states (with lower spiking variability). Here, using both new and publicly available data, we test independent signatures of criticality and show that a phase transition occurs in an intermediate value of spiking variability, in both anesthetized and freely moving animals. The critical exponents point to a universality class different from mean-field directed percolation. Importantly, as the cortex hovers around this critical point, the avalanche exponents follow a linear relation that encompasses previous experimental results from different setups and is reproduced by a model.

Single-Cell Membrane Potential Fluctuations Evince Network Scale-Freeness and Quasicriticality

What information single neurons receive about general neural circuit activity is a fundamental question for neuroscience. Somatic membrane potential (V _m) fluctuations are driven by the convergence of synaptic inputs from a diverse cross-section of upstream neurons. Furthermore, neural activity is often scale-free, implying that some measurements should be the same, whether taken at large or small scales. Together, convergence and scale-freeness support the hypothesis that single V _m recordings carry useful information about high-dimensional cortical activity. Conveniently, the theory of "critical branching networks" (one purported explanation for scale-freeness) provides testable predictions about scale-free measurements that are readily applied to V _m fluctuations. To investigate, we obtained whole-cell current-clamp recordings of pyramidal neurons in visual cortex of turtles with unknown genders. We isolated fluctuations in V _m below the firing threshold and analyzed them by adapting the definition of "neuronal avalanches" (i.e., spurts of population spiking). The V _m fluctuations which we analyzed were scale-free and consistent with critical branching. These findings recapitulated results from large-scale cortical population data obtained separately in complementary experiments using microelectrode arrays described previously (Shew et al., 2015). Simultaneously recorded single-unit local field potential did not provide a good match, demonstrating the specific utility of V _m Modeling shows that estimation of dynamical network properties from neuronal inputs is most accurate when networks are structured as critical branching networks. In conclusion, these findings extend evidence of critical phenomena while also establishing subthreshold pyramidal neuron V _m fluctuations as an informative gauge of high-dimensional cortical population activity.SIGNIFICANCE STATEMENT The relationship between membrane potential (V _m) dynamics of single neurons and population dynamics is indispensable to understanding cortical circuits. Just as important to the biophysics of computation are emergent properties such as scale-freeness, where critical branching networks offer insight. This report makes progress on both fronts by comparing statistics from single-neuron whole-cell recordings with population statistics obtained with microelectrode arrays. Not only are fluctuations of somatic V _m scale-free, they match fluctuations of population activity. Thus, our results demonstrate appropriation of the brain's own subsampling method (convergence of synaptic inputs) while extending the range of fundamental evidence for critical phenomena in neural systems from the previously observed mesoscale (fMRI, LFP, population spiking) to the microscale, namely, V _m fluctuations.

Modeling neuronal avalanches and long-range temporal correlations at the emergence of collective oscillations: Continuously varying exponents mimic M/EEG results

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.

Pattern phase diagram of spiking neurons on spatial networks

We study an abstracted model of neuronal activity via numerical simulation and report spatiotemporal pattern formation and criticallike dynamics. A population of pulse coupled, discretized, relaxation oscillators is simulated over networks with varying edge density and spatial embeddedness. For intermediate edge density and sufficiently strong spatial embeddedness, we observe a spatiotemporal pattern in the field of oscillator phases, visually resembling the surface of a frothing liquid. Increasing the edge density results in a distribution of neuronal avalanche sizes which follows a power law with exponent one (Zipf's law). Further increasing edge density yields metastability between pattern formation and synchronization, before transitioning entirely into synchrony.

Robustness of Griffiths effects in homeostatic connectome models

I provide numerical evidence for the robustness of the Griffiths phase (GP) reported previously in dynamical threshold model simulations on a large human brain network with N=836733 connected nodes. The model, with equalized network sensitivity, is extended in two ways: introduction of refractory states or by randomized time-dependent thresholds. The nonuniversal power-law dynamics in an extended control parameter region survives these modifications for a short refractory state and weak disorder. In case of temporal disorder the GP shrinks and for stronger heterogeneity disappears, leaving behind a mean-field type of critical transition. Activity avalanche size distributions below the critical point decay faster than in the original model, but the addition of inhibitory interactions sets it back to the range of experimental values.