Power-law statistics and universal scaling in the absence of criticality

Sub-Second Fluctuation between Top-Down and Bottom-Up Modes Distinguishes Diverse Human Brain States

Information continuously flows between regions of the human brain, exhibiting distinct patterns that dynamically shift across states of consciousness, cognitive modes, and neuropsychiatric conditions. In this study, we introduce Relative Phase Analysis (RPA), a method that leverages phase-lead/lag relationships to reveal the real-time dynamics of dominant directional patterns and their rapid transitions. We demonstrate that the human brain switches on a sub-second timescale between a top-down mode-where anterior regions drive posterior activity-and a bottom-up mode, characterized by reverse directionality. These dynamics are most pronounced during full consciousness and gradually become less distinct as awareness diminishes. Furthermore, we find from simultaneous EEG-fMRI recordings that the top-down mode is expressed when higher-order cognitive networks are more active while the bottom-up mode is expressed when sensory systems are more active. Moreover, comparisons of an attention deficit hyperactivity disorder (ADHD) inattentive cohort with typically developing individuals reveal distinct imbalances in these transition dynamics, highlighting the potential of RPA as a diagnostic biomarker. Complementing our empirical findings, a coupled-oscillator model of the structural brain network recapitulates these emergent patterns, suggesting that such directional modes and transitions may arise naturally from inter-regional neural interactions. Altogether, this study provides a framework for understanding whole-brain dynamics in real-time and identifies sub-second fluctuations in top-down versus bottom-up directionality as a fundamental mechanism underlying human information processing.

Focal Infrared Neural Stimulation Propagates Dynamical Transformations in Auditory Cortex

Significance

Infrared neural stimulation (INS) has emerged as a potent neuromodulation technology, offering safe and focal stimulation with superior spatial recruitment profiles compared to conventional electrical methods. However, the neural dynamics induced by INS stimulation remain poorly understood. Elucidating these dynamics will help develop new INS stimulation paradigms and advance its clinical application.

Aim

In this study, we assessed the local network dynamics of INS entrainment in the auditory thalamocortical circuit using the chronically implanted rat model; our approach focused on measuring INS energy-based local field potential (LFP) recruitment induced by focal thalamocortical stimulation. We further characterized linear and nonlinear oscillatory LFP activity in response to single-pulse and periodic INS and performed spectral decomposition to uncover specific LFP band entrainment to INS. Finally, we examined spike-field transformations across the thalamocortical synapse using spike-LFP coherence coupling.

Results

We found that INS significantly increases LFP amplitude as a log-linear function of INS energy per pulse, primarily entraining to LFP β and γ bands with synchrony extending to 200 Hz in some cases. A subset of neurons demonstrated nonlinear, chaotic oscillations linked to information transfer across cortical circuits. Finally, we utilized spike-field coherences to correlate spike coupling to LFP frequency band activity and suggest an energy-dependent model of network activation resulting from INS stimulation.

Conclusions

We show that INS reliably drives robust network activity and can potently modulate cortical field potentials across a wide range of frequencies in a stimulus parameter-dependent manner. Based on these results, we propose design principles for developing full coverage, all-optical thalamocortical auditory neuroprostheses.

Excitation-Inhibition Balance, Neural Criticality, and Activities in Neuronal Circuits

Neural activities in local circuits exhibit complex and multilevel dynamic features. Individual neurons spike irregularly, which is believed to originate from receiving balanced amounts of excitatory and inhibitory inputs, known as the excitation-inhibition balance. The spatial-temporal cascades of clustered neuronal spikes occur in variable sizes and durations, manifested as neural avalanches with scale-free features. These may be explained by the neural criticality hypothesis, which posits that neural systems operate around the transition between distinct dynamic states. Here, we summarize the experimental evidence for and the underlying theory of excitation-inhibition balance and neural criticality. Furthermore, we review recent studies of excitatory-inhibitory networks with synaptic kinetics as a simple solution to reconcile these two apparently distinct theories in a single circuit model. This provides a more unified understanding of multilevel neural activities in local circuits, from spontaneous to stimulus-response dynamics.

The vitals for steady nucleation maps of spontaneous spiking coherence in autonomous two-dimensional neuronal networks

Thin pancake-like neuronal networks cultured on top of a planar microelectrode array have been extensively tried out in neuroengineering, as a substrate for the mobile robot's control unit, i.e., as a cyborg's brain. Most of these attempts failed due to intricate self-organizing dynamics in the neuronal systems. In particular, the networks may exhibit an emergent spatial map of steady nucleation sites ("n-sites") of spontaneous population spikes. Being unpredictable and independent of the surface electrode locations, the n-sites drastically change local ability of the network to generate spikes. Here, using a spiking neuronal network model with generative spatially-embedded connectome, we systematically show in simulations that the number, location, and relative activity of spontaneously formed n-sites ("the vitals") crucially depend on the samplings of three distributions: (1) the network distribution of neuronal excitability, (2) the distribution of connections between neurons of the network, and (3) the distribution of maximal amplitudes of a single synaptic current pulse. Moreover, blocking the dynamics of a small fraction (about 4%) of non-pacemaker neurons having the highest excitability was enough to completely suppress the occurrence of population spikes and their n-sites. This key result is explained theoretically. Remarkably, the n-sites occur taking into account only short-term synaptic plasticity, i.e., without a Hebbian-type plasticity. As the spiking network model used in this study is strictly deterministic, all simulation results can be accurately reproduced. The model, which has already demonstrated a very high richness-to-complexity ratio, can also be directly extended into the three-dimensional case, e.g., for targeting peculiarities of spiking dynamics in cerebral (or brain) organoids. We recommend the model as an excellent illustrative tool for teaching network-level computational neuroscience, complementing a few benchmark models.

Failure in a population: Tauopathy disrupts homeostatic set-points in emergent dynamics despite stability in the constituent neurons

Homeostatic regulation of neuronal activity is essential for robust computation; set-points, such as firing rate, are actively stabilized to compensate for perturbations. The disruption of brain function central to neurodegenerative disease likely arises from impairments of computationally essential set-points. Here, we systematically investigated the effects of tau-mediated neurodegeneration on all known set-points in neuronal activity. We continuously tracked hippocampal neuronal activity across the lifetime of a mouse model of tauopathy. We were unable to detect effects of disease in measures of single-neuron firing activity. By contrast, as tauopathy progressed, there was disruption of network-level neuronal activity, quantified by measuring neuronal pairwise interactions and criticality, a homeostatically controlled, ideal computational regime. Deviations in criticality correlated with symptoms, predicted underlying anatomical pathology, occurred in a sleep-wake-dependent manner, and could be used to reliably classify an animal's genotype. This work illustrates how neurodegeneration may disrupt the computational capacity of neurobiological systems.

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.

The recovery of parabolic avalanches in spatially subsampled neuronal networks at criticality

Scaling relationships are key in characterizing complex systems at criticality. In the brain, they are evident in neuronal avalanches-scale-invariant cascades of neuronal activity quantified by power laws. Avalanches manifest at the cellular level as cascades of neuronal groups that fire action potentials simultaneously. Such spatiotemporal synchronization is vital to theories on brain function yet avalanche synchronization is often underestimated when only a fraction of neurons is observed. Here, we investigate biases from fractional sampling within a balanced network of excitatory and inhibitory neurons with all-to-all connectivity and critical branching process dynamics. We focus on how mean avalanche size scales with avalanche duration. For parabolic avalanches, this scaling is quadratic, quantified by the scaling exponent, χ = 2, reflecting rapid spatial expansion of simultaneous neuronal firing over short durations. However, in networks sampled fractionally, χ is significantly lower. We demonstrate that applying temporal coarse-graining and increasing a minimum threshold for coincident firing restores χ = 2, even when as few as 0.1% of neurons are sampled. This correction crucially depends on the network being critical and fails for near sub- and supercritical branching dynamics. Using cellular 2-photon imaging, our approach robustly identifies χ = 2 over a wide parameter regime in ongoing neuronal activity from frontal cortex of awake mice. In contrast, the common 'crackling noise' approach fails to determine χ under similar sampling conditions at criticality. Our findings overcome scaling bias from fractional sampling and demonstrate rapid, spatiotemporal synchronization of neuronal assemblies consistent with scale-invariant, parabolic avalanches at criticality.

Maximum-entropy-based metrics for quantifying critical dynamics in spiking neuron data

An important working hypothesis to investigate brain activity is whether it operates in a critical regime. Recently, maximum-entropy phenomenological models have emerged as an alternative way of identifying critical behavior in neuronal data sets. In the present paper, we investigate the signatures of criticality from a firing rate-based maximum-entropy approach on data sets generated by computational models, and we compare them to experimental results. We found that the maximum entropy approach consistently identifies critical behavior around the phase transition in models and rules out criticality in models without phase transition. The maximum-entropy-model results are compatible with results for cortical data from urethane-anesthetized rats data, providing further support for criticality in the brain.

Nondifferentiable activity in the brain

Spike raster plots of numerous neurons show vertical stripes, indicating that neurons exhibit synchronous activity in the brain. We seek to determine whether these coherent dynamics are caused by smooth brainwave activity or by something else. By analyzing biological data, we find that their cross-correlograms exhibit not only slow undulation but also a cusp at the origin, in addition to possible signs of monosynaptic connectivity. Here we show that undulation emerges if neurons are subject to smooth brainwave oscillations while a cusp results from nondifferentiable fluctuations. While modern analysis methods have achieved good connectivity estimation by adapting the models to slow undulation, they still make false inferences due to the cusp. We devise a new analysis method that may solve both problems. We also demonstrate that oscillations and nondifferentiable fluctuations may emerge in simulations of large-scale neural networks.

The recovery of parabolic avalanches in spatially subsampled neuronal networks at criticality

Scaling relationships are key in characterizing complex systems at criticality. In the brain, they are evident in neuronal avalanches-scale-invariant cascades of neuronal activity quantified by power laws. Avalanches manifest at the cellular level as cascades of neuronal groups that fire action potentials simultaneously. Such spatiotemporal synchronization is vital to theories on brain function yet avalanche synchronization is often underestimated when only a fraction of neurons is observed. Here, we investigate biases from fractional sampling within a balanced network of excitatory and inhibitory neurons with all-to-all connectivity and critical branching process dynamics. We focus on how mean avalanche size scales with avalanche duration. For parabolic avalanches, this scaling is quadratic, quantified by the scaling exponent, χ = 2 , reflecting rapid spatial expansion of simultaneous neuronal firing over short durations. However, in networks sampled fractionally, χ is significantly lower. We demonstrate that applying temporal coarse-graining and increasing a minimum threshold for coincident firing restores χ = 2 , even when as few as 0.1% of neurons are sampled. This correction crucially depends on the network being critical and fails for near sub- and supercritical branching dynamics. Using cellular 2-photon imaging, our approach robustly identifies χ = 2 over a wide parameter regime in ongoing neuronal activity from frontal cortex of awake mice. In contrast, the common 'crackling noise' approach fails to determine χ under similar sampling conditions at criticality. Our findings overcome scaling bias from fractional sampling and demonstrate rapid, spatiotemporal synchronization of neuronal assemblies consistent with scale-invariant, parabolic avalanches at criticality.

Neural criticality from effective latent variables

Observations of power laws in neural activity data have raised the intriguing notion that brains may operate in a critical state. One example of this critical state is 'avalanche criticality', which has been observed in various systems, including cultured neurons, zebrafish, rodent cortex, and human EEG. More recently, power laws were also observed in neural populations in the mouse under an activity coarse-graining procedure, and they were explained as a consequence of the neural activity being coupled to multiple latent dynamical variables. An intriguing possibility is that avalanche criticality emerges due to a similar mechanism. Here, we determine the conditions under which latent dynamical variables give rise to avalanche criticality. We find that populations coupled to multiple latent variables produce critical behavior across a broader parameter range than those coupled to a single, quasi-static latent variable, but in both cases, avalanche criticality is observed without fine-tuning of model parameters. We identify two regimes of avalanches, both critical but differing in the amount of information carried about the latent variable. Our results suggest that avalanche criticality arises in neural systems in which activity is effectively modeled as a population driven by a few dynamical variables and these variables can be inferred from the population activity.

Sleep restores an optimal computational regime in cortical networks

Sleep is assumed to subserve homeostatic processes in the brain; however, the set point around which sleep tunes circuit computations is unknown. Slow-wave activity (SWA) is commonly used to reflect the homeostatic aspect of sleep; although it can indicate sleep pressure, it does not explain why animals need sleep. This study aimed to assess whether criticality may be the computational set point of sleep. By recording cortical neuron activity continuously for 10-14 d in freely behaving rats, we show that normal waking experience progressively disrupts criticality and that sleep functions to restore critical dynamics. Criticality is perturbed in a context-dependent manner, and waking experience is causal in driving these effects. The degree of deviation from criticality predicts future sleep/wake behavior more accurately than SWA, behavioral history or other neural measures. Our results demonstrate that perturbation and recovery of criticality is a network homeostatic mechanism consistent with the core, restorative function of sleep.

Topological data analysis of the firings of a network of stochastic spiking neurons

Topological data analysis is becoming more and more popular in recent years. It has found various applications in many different fields, for its convenience in analyzing and understanding the structure and dynamic of complex systems. We used topological data analysis to analyze the firings of a network of stochastic spiking neurons, which can be in a sub-critical, critical, or super-critical state depending on the value of the control parameter. We calculated several topological features regarding Betti curves and then analyzed the behaviors of these features, using them as inputs for machine learning to discriminate the three states of the network.

[Connectome in stroke patients]

Stroke is the main cause of disability among neurological diseases. There are questions of the accuracy of topical diagnosis and rehabilitation prognosis in clinical practice. Answers to these questions may be given by an approach to the study of the nervous system as a dynamic network consisting of a set of brain regions with anatomical and functional connections between them. Active study of the connectome in neurological patients in recent years became possible due to the availability of noninvasive neuroimaging methods. This review covers types of connectome and most accessible methods of obtaining research data for their construction in a neurological hospital. The review also describes resting-state networks that reflect basic brain activity in the absence of tasks. Resting-state connectivity can be used for the diagnosis of patients with severe neurological deficits. Also, changes in resting-state connectivity may indicate recovery after a stroke. The connectome analysis uses graph theory, representing the nervous system as a set of nodes and connections between them, and providing a mathematical framework allowing to study it, and methods of algebraic topology that expand the possibilities of analyzing neuroimaging data beyond graph theory. Attention is paid to the concept of self-organized criticality, which describes the brain as a system located near the critical point, where the transmission of information is most optimized. Also presented are data from studies of self-organized criticality in relation to the dynamics of recovery of patients after stroke.

Changes in functional connectivity preserve scale-free neuronal and behavioral dynamics

Does the brain optimize itself for storage and transmission of information, and if so, how? The critical brain hypothesis is based in statistical physics and posits that the brain self-tunes its dynamics to a critical point or regime to maximize the repertoire of neuronal responses. Yet, the robustness of this regime, especially with respect to changes in the functional connectivity, remains an unsolved fundamental challenge. Here, we show that both scale-free neuronal dynamics and self-similar features of behavioral dynamics persist following significant changes in functional connectivity. Specifically, we find that the psychedelic compound ibogaine that is associated with an altered state of consciousness fundamentally alters the functional connectivity in the retrosplenial cortex of mice. Yet, the scale-free statistics of movement and of neuronal avalanches among behaviorally related neurons remain largely unaltered. This indicates that the propagation of information within biological neural networks is robust to changes in functional organization of subpopulations of neurons, opening up a new perspective on how the adaptive nature of functional networks may lead to optimality of information transmission in the brain.

Critical Growth of Cerebral Tissue in Organoids: Theory and Experiments

We develop a Fokker-Planck theory of tissue growth with three types of cells (symmetrically dividing, asymmetrically dividing, and nondividing) as main agents to study the growth dynamics of human cerebral organoids. Fitting the theory to lineage tracing data obtained in next generation sequencing experiments, we show that the growth of cerebral organoids is a critical process. We derive analytical expressions describing the time evolution of clonal lineage sizes and show how power-law distributions arise in the limit of long times due to the vanishing of a characteristic growth scale. We discuss that the independence of critical growth on initial conditions could be biologically advantageous.

Fluctuating landscapes and heavy tails in animal behavior

Animal behavior is shaped by a myriad of mechanisms acting on a wide range of scales. This immense variability hampers quantitative reasoning and renders the identification of universal principles elusive. Through data analysis and theory, we here show that slow non-ergodic drives generally give rise to heavy-tailed statistics in behaving animals. We leverage high-resolution recordings of C. elegans locomotion to extract a self-consistent reduced order model for an inferred reaction coordinate, bridging from sub-second chaotic dynamics to long-lived stochastic transitions among metastable states. The slow mode dynamics exhibits heavy-tailed first passage time distributions and correlation functions, and we show that such heavy tails can be explained by dynamics on a time-dependent potential landscape. Inspired by these results, we introduce a generic model in which we separate faster mixing modes that evolve on a quasi-stationary potential, from slower non-ergodic modes that drive the potential landscape, and reflect slowly varying internal states. We show that, even for simple potential landscapes, heavy tails emerge when barrier heights fluctuate slowly and strongly enough. In particular, the distribution of first passage times and the correlation function can asymptote to a power law, with related exponents that depend on the strength and nature of the fluctuations. We support our theoretical findings through direct numerical simulations.

A pattern of cognitive resource disruptions in childhood psychopathology

The Hurst exponent (H) isolated in fractal analyses of neuroimaging time series is implicated broadly in cognition. Within this literature, H is associated with multiple mental disorders, suggesting that H is transdimensionally associated with psychopathology. Here, we unify these results and demonstrate a pattern of decreased H with increased general psychopathology and attention-deficit/hyperactivity factor scores during a working memory task in 1,839 children. This pattern predicts current and future cognitive performance in children and some psychopathology in 703 adults. This pattern also defines psychological and functional axes associating psychopathology with an imbalance in resource allocation between fronto-parietal and sensorimotor regions, driven by reduced resource allocation to fronto-parietal regions. This suggests the hypothesis that impaired working memory function in psychopathology follows from a reduced cognitive resource pool and a reduction in resources allocated to the task at hand.

Tauopathy severely disrupts homeostatic set-points in emergent neural dynamics but not in the activity of individual neurons

The homeostatic regulation of neuronal activity is essential for robust computation; key set-points, such as firing rate, are actively stabilized to compensate for perturbations. From this perspective, the disruption of brain function central to neurodegenerative disease should reflect impairments of computationally essential set-points. Despite connecting neurodegeneration to functional outcomes, the impact of disease on set-points in neuronal activity is unknown. Here we present a comprehensive, theory-driven investigation of the effects of tau-mediated neurodegeneration on homeostatic set-points in neuronal activity. In a mouse model of tauopathy, we examine 27,000 hours of hippocampal recordings during free behavior throughout disease progression. Contrary to our initial hypothesis that tauopathy would impact set-points in spike rate and variance, we found that cell-level set-points are resilient to even the latest stages of disease. Instead, we find that tauopathy disrupts neuronal activity at the network-level, which we quantify using both pairwise measures of neuron interactions as well as measurement of the network's nearness to criticality, an ideal computational regime that is known to be a homeostatic set-point. We find that shifts in network criticality 1) track with symptoms, 2) predict underlying anatomical and molecular pathology, 3) occur in a sleep/wake dependent manner, and 4) can be used to reliably classify an animal's genotype. Our data suggest that the critical set-point is intact, but that homeostatic machinery is progressively incapable of stabilizing hippocampal networks, particularly during waking. This work illustrates how neurodegenerative processes can impact the computational capacity of neurobiological systems, and suggest an important connection between molecular pathology, circuit function, and animal behavior.

Cluster scaling and critical points: A cautionary tale

Many systems in nature are conjectured to exist at a critical point, including the brain and earthquake faults. The primary reason for this conjecture is that the distribution of clusters (avalanches of firing neurons in the brain or regions of slip in earthquake faults) can be described by a power law. Because there are other mechanisms such as 1/f noise that can produce power laws, other criteria that the cluster critical exponents must satisfy can be used to conclude whether or not the observed power-law behavior indicates an underlying critical point rather than an alternate mechanism. We show how a possible misinterpretation of the cluster scaling data can lead one to incorrectly conclude that the measured critical exponents do not satisfy these criteria. Examples of the possible misinterpretation of the data for one-dimensional random site percolation and the one-dimensional Ising model are presented. We stress that the interpretation of a power-law cluster distribution indicating the presence of a critical point is subtle and its misinterpretation might lead to the abandonment of a promising area of research.

Critical dynamics arise during structured information presentation within embodied in vitro neuronal networks

Understanding how brains process information is an incredibly difficult task. Amongst the metrics characterising information processing in the brain, observations of dynamic near-critical states have generated significant interest. However, theoretical and experimental limitations associated with human and animal models have precluded a definite answer about when and why neural criticality arises with links from attention, to cognition, and even to consciousness. To explore this topic, we used an in vitro neural network of cortical neurons that was trained to play a simplified game of 'Pong' to demonstrate Synthetic Biological Intelligence (SBI). We demonstrate that critical dynamics emerge when neural networks receive task-related structured sensory input, reorganizing the system to a near-critical state. Additionally, better task performance correlated with proximity to critical dynamics. However, criticality alone is insufficient for a neuronal network to demonstrate learning in the absence of additional information regarding the consequences of previous actions. These findings offer compelling support that neural criticality arises as a base feature of incoming structured information processing without the need for higher order cognition.

Cell-autonomous diversification in bacteria arises from calcium dynamics self-organizing at a critical point

How dynamic bacterial calcium is regulated, with kinetics faster than typical mechanisms of cellular adaptation, is unknown. We discover bacterial calcium fluctuations are temporal-fractals resulting from a property known as self-organized criticality (SOC). SOC processes are poised at a phase transition separating ordered and chaotic dynamical regimes and are observed in many natural and anthropogenic systems. SOC in bacterial calcium emerges due to calcium channel coupling mediated via membrane voltage. Environmental or genetic perturbations modify calcium dynamics and the critical exponent suggesting a continuum of critical attractors. Moving along this continuum alters the collective information capacity of bacterial populations. We find that the stochastic transition from motile to sessile lifestyle is partially mediated by SOC-governed calcium fluctuations through the regulation of c-di-GMP. In summary, bacteria co-opt the physics of phase transitions to maintain dynamic calcium equilibrium, and this enables cell-autonomous population diversification during surface colonization by leveraging the stochasticity inherent at a boundary between phases.

Neuronal avalanches: Sandpiles of self-organized criticality or critical dynamics of brain waves?

Analytical expressions for scaling of brain wave spectra derived from the general nonlinear wave Hamiltonian form show excellent agreement with experimental "neuronal avalanche" data. The theory of the weakly evanescent nonlinear brain wave dynamics [Phys. Rev. Research 2, 023061 (2020); J. Cognitive Neurosci. 32, 2178 (2020)] reveals the underlying collective processes hidden behind the phenomenological statistical description of the neuronal avalanches and connects together the whole range of brain activity states, from oscillatory wave-like modes, to neuronal avalanches, to incoherent spiking, showing that the neuronal avalanches are just the manifestation of the different nonlinear side of wave processes abundant in cortical tissue. In a more broad way these results show that a system of wave modes interacting through all possible combinations of the third order nonlinear terms described by a general wave Hamiltonian necessarily produces anharmonic wave modes with temporal and spatial scaling properties that follow scale free power laws. To the best of our knowledge this has never been reported in the physical literature and may be applicable to many physical systems that involve wave processes and not just to neuronal avalanches.

Parabolic avalanche scaling in the synchronization of cortical cell assemblies

Neurons in the cerebral cortex fire coincident action potentials during ongoing activity and in response to sensory inputs. These synchronized cell assemblies are fundamental to cortex function, yet basic dynamical aspects of their size and duration are largely unknown. Using 2-photon imaging of neurons in the superficial cortex of awake mice, we show that synchronized cell assemblies organize as scale-invariant avalanches that quadratically grow with duration. The quadratic avalanche scaling was only found for correlated neurons, required temporal coarse-graining to compensate for spatial subsampling of the imaged cortex, and suggested cortical dynamics to be critical as demonstrated in simulations of balanced E/I-networks. The corresponding time course of an inverted parabola with exponent of χ = 2 described cortical avalanches of coincident firing for up to 5 s duration over an area of 1 mm². These parabolic avalanches maximized temporal complexity in the ongoing activity of prefrontal and somatosensory cortex and in visual responses of primary visual cortex. Our results identify a scale-invariant temporal order in the synchronization of highly diverse cortical cell assemblies in the form of parabolic avalanches.

Microscale Neuronal Activity Collectively Drives Chaotic and Inflexible Dynamics at the Macroscale in Seizures

Neuronal activity propagates through the network during seizures, engaging brain dynamics at multiple scales. Such propagating events can be described through the avalanches framework, which can relate spatiotemporal activity at the microscale with global network properties. Interestingly, propagating avalanches in healthy networks are indicative of critical dynamics, where the network is organized to a phase transition, which optimizes certain computational properties. Some have hypothesized that the pathologic brain dynamics of epileptic seizures are an emergent property of microscale neuronal networks collectively driving the brain away from criticality. Demonstrating this would provide a unifying mechanism linking microscale spatiotemporal activity with emergent brain dysfunction during seizures. Here, we investigated the effect of drug-induced seizures on critical avalanche dynamics, using in vivo whole-brain two-photon imaging of GCaMP6s larval zebrafish (males and females) at single neuron resolution. We demonstrate that single neuron activity across the whole brain exhibits a loss of critical statistics during seizures, suggesting that microscale activity collectively drives macroscale dynamics away from criticality. We also construct spiking network models at the scale of the larval zebrafish brain, to demonstrate that only densely connected networks can drive brain-wide seizure dynamics away from criticality. Importantly, such dense networks also disrupt the optimal computational capacities of critical networks, leading to chaotic dynamics, impaired network response properties and sticky states, thus helping to explain functional impairments during seizures. This study bridges the gap between microscale neuronal activity and emergent macroscale dynamics and cognitive dysfunction during seizures.SIGNIFICANCE STATEMENT Epileptic seizures are debilitating and impair normal brain function. It is unclear how the coordinated behavior of neurons collectively impairs brain function during seizures. To investigate this we perform fluorescence microscopy in larval zebrafish, which allows for the recording of whole-brain activity at single-neuron resolution. Using techniques from physics, we show that neuronal activity during seizures drives the brain away from criticality, a regime that enables both high and low activity states, into an inflexible regime that drives high activity states. Importantly, this change is caused by more connections in the network, which we show disrupts the ability of the brain to respond appropriately to its environment. Therefore, we identify key neuronal network mechanisms driving seizures and concurrent cognitive dysfunction.

Time-resolved correlation of distributed brain activity tracks E-I balance and accounts for diverse scale-free phenomena

Much of systems neuroscience posits the functional importance of brain activity patterns that lack natural scales of sizes, durations, or frequencies. The field has developed prominent, and sometimes competing, explanations for the nature of this scale-free activity. Here, we reconcile these explanations across species and modalities. First, we link estimates of excitation-inhibition (E-I) balance with time-resolved correlation of distributed brain activity. Second, we develop an unbiased method for sampling time series constrained by this time-resolved correlation. Third, we use this method to show that estimates of E-I balance account for diverse scale-free phenomena without need to attribute additional function or importance to these phenomena. Collectively, our results simplify existing explanations of scale-free brain activity and provide stringent tests on future theories that seek to transcend these explanations.

The fractal brain: scale-invariance in structure and dynamics

The past 40 years have witnessed extensive research on fractal structure and scale-free dynamics in the brain. Although considerable progress has been made, a comprehensive picture has yet to emerge, and needs further linking to a mechanistic account of brain function. Here, we review these concepts, connecting observations across different levels of organization, from both a structural and functional perspective. We argue that, paradoxically, the level of cortical circuits is the least understood from a structural point of view and perhaps the best studied from a dynamical one. We further link observations about scale-freeness and fractality with evidence that the environment provides constraints that may explain the usefulness of fractal structure and scale-free dynamics in the brain. Moreover, we discuss evidence that behavior exhibits scale-free properties, likely emerging from similarly organized brain dynamics, enabling an organism to thrive in an environment that shares the same organizational principles. Finally, we review the sparse evidence for and try to speculate on the functional consequences of fractality and scale-freeness for brain computation. These properties may endow the brain with computational capabilities that transcend current models of neural computation and could hold the key to unraveling how the brain constructs percepts and generates behavior.

Critical brain wave dynamics of neuronal avalanches

Analytical expressions for scaling of brain wave spectra derived from the general non-linear wave Hamiltonian form show excellent agreement with experimental "neuronal avalanche" data. The theory of the weakly evanescent non-linear brain wave dynamics reveals the underlying collective processes hidden behind the phenomenological statistical description of the neuronal avalanches and connects together the whole range of brain activity states, from oscillatory wave-like modes, to neuronal avalanches, to incoherent spiking, showing that the neuronal avalanches are just the manifestation of the different non-linear side of wave processes abundant in cortical tissue. In a more broad way these results show that a system of wave modes interacting through all possible combinations of the third order non-linear terms described by a general wave Hamiltonian necessarily produces anharmonic wave modes with temporal and spatial scaling properties that follow scale free power laws. To the best of our knowledge this has never been reported in the physical literature and may be applicable to many physical systems that involve wave processes and not just to neuronal avalanches.

Finite-size correlation behavior near a critical point: A simple metric for monitoring the state of a neural network

In this article, a correlation metric κ_{c} is proposed for the inference of the dynamical state of neuronal networks. κ_{C} is computed from the scaling of the correlation length with the size of the observation region, which shows qualitatively different behavior near and away from the critical point of a continuous phase transition. The implementation is first studied on a neuronal network model, where the results of this new metric coincide with those obtained from neuronal avalanche analysis, thus well characterizing the critical state of the network. The approach is further tested with brain optogenetic recordings in behaving mice from a publicly available database. Potential applications and limitations for its use with currently available optical imaging techniques are discussed.

Theoretical foundations of studying criticality in the brain

Criticality is hypothesized as a physical mechanism underlying efficient transitions between cortical states and remarkable information-processing capacities in the brain. While considerable evidence generally supports this hypothesis, nonnegligible controversies persist regarding the ubiquity of criticality in neural dynamics and its role in information processing. Validity issues frequently arise during identifying potential brain criticality from empirical data. Moreover, the functional benefits implied by brain criticality are frequently misconceived or unduly generalized. These problems stem from the nontriviality and immaturity of the physical theories that analytically derive brain criticality and the statistic techniques that estimate brain criticality from empirical data. To help solve these problems, we present a systematic review and reformulate the foundations of studying brain criticality, that is, ordinary criticality (OC), quasi-criticality (qC), self-organized criticality (SOC), and self-organized quasi-criticality (SOqC), using the terminology of neuroscience. We offer accessible explanations of the physical theories and statistical techniques of brain criticality, providing step-by-step derivations to characterize neural dynamics as a physical system with avalanches. We summarize error-prone details and existing limitations in brain criticality analysis and suggest possible solutions. Moreover, we present a forward-looking perspective on how optimizing the foundations of studying brain criticality can deepen our understanding of various neuroscience questions.

Scale-free avalanche dynamics possibly generated by randomly jumping among many stable states

A large amount of research has used the scale-free statistics of neuronal avalanches as a signature of the criticality of neural systems, which bears criticisms. For instance, the work of Touboul and Destexhe demonstrated that non-critical systems could also display such scale-free dynamics, which passed their rigorous statistical analyses. In this paper, we show that a fully connected stochastic neural network may also generate scale-free dynamics simply by jumping among many stable states.

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.

Turing's cascade instability supports the coordination of the mind, brain, and behavior

Turing inspired a computer metaphor of the mind and brain that has been handy and has spawned decades of empirical investigation, but he did much more and offered behavioral and cognitive sciences another metaphor-that of the cascade. The time has come to confront Turing's cascading instability, which suggests a geometrical framework driven by power laws and can be studied using multifractal formalism and multiscale probability density function analysis. Here, we review a rapidly growing body of scientific investigations revealing signatures of cascade instability and their consequences for a perceiving, acting, and thinking organism. We review work related to executive functioning (planning to act), postural control (bodily poise for turning plans into action), and effortful perception (action to gather information in a single modality and action to blend multimodal information). We also review findings on neuronal avalanches in the brain, specifically about neural participation in body-wide cascades. Turing's cascade instability blends the mind, brain, and behavior across space and time scales and provides an alternative to the dominant computer metaphor.

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.

Self-organized criticality as a framework for consciousness: A review study

Objective

No current model of consciousness is univocally accepted on either theoretical or empirical grounds, and the need for a solid unifying framework is evident. Special attention has been given to the premise that self-organized criticality (SOC) is a fundamental property of neural system. SOC provides a competitive model to describe the physical mechanisms underlying spontaneous brain activity, and thus, critical dynamics were proposed as general gauges of information processing representing a strong candidate for a surrogate measure of consciousness. As SOC could be a neurodynamical framework, which may be able to bring together existing theories and experimental evidence, the purpose of this work was to provide a comprehensive overview of progress of research on SOC in association with consciousness.

Methods

A comprehensive search of publications on consciousness and SOC published between 1998 and 2021 was conducted. The Web of Science database was searched, and annual number of publications and citations, type of articles, and applied methods were determined.

Results

A total of 71 publications were identified. The annual number of citations steadily increased over the years. Original articles comprised 50.7% and reviews/theoretical articles 43.6%. Sixteen studies reported on human data and in seven studies data were recorded in animals. Computational models were utilized in n = 12 studies. EcoG data were assessed in n = 4 articles, fMRI in n = 4 studies, and EEG/MEG in n = 10 studies. Notably, different analytical tools were applied in the EEG/MEG studies to assess a surrogate measure of criticality such as the detrended fluctuation analysis, the pair correlation function, parameters from the neuronal avalanche analysis and the spectral exponent.

Conclusion

Recent studies pointed out agreements of critical dynamics with the current most influencing theories in the field of consciousness research, the global workspace theory and the integrated information theory. Thus, the framework of SOC as a neurodynamical parameter for consciousness seems promising. However, identified experimental work was small in numbers, and a heterogeneity of applied analytical tools as a surrogate measure of criticality was observable, which limits the generalizability of findings.

Memory in quantum dot blinking

The photoluminescence intermittency (blinking) of quantum dots is interesting because it is an easily measured quantum process whose transition statistics cannot be explained by Fermi's golden rule. Commonly, the transition statistics are power-law distributed, implying that quantum dots possess at least trivial memories. By investigating the temporal correlations in the blinking data, we demonstrate with high statistical confidence that there is nontrivial memory between the on and off brightness duration data of blinking quantum dots. We define nontrivial memory to be statistical complexity greater than one. We show that this memory cannot be discovered using the transition distribution. We show by simulation that this memory does not arise from standard data manipulations. Finally, we conclude that at least three physical mechanisms can explain the measured nontrivial memory: (1) storage of state information in the chemical structure of a quantum dot; (2) the existence of more than two intensity levels in a quantum dot; and (3) the overlap in the intensity distributions of the quantum dot states, which arises from fundamental photon statistics.

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.

May the 4C's be with you: an overview of complexity-inspired frameworks for analysing resting-state neuroimaging data

Competing and complementary models of resting-state brain dynamics contribute to our phenomenological and mechanistic understanding of whole-brain coordination and communication, and provide potential evidence for differential brain functioning associated with normal and pathological behaviour. These neuroscientific theories stem from the perspectives of physics, engineering, mathematics and psychology and create a complicated landscape of domain-specific terminology and meaning, which, when used outside of that domain, may lead to incorrect assumptions and conclusions within the neuroscience community. Here, we review and clarify the key concepts of connectivity, computation, criticality and coherence-the 4C's-and outline a potential role for metastability as a common denominator across these propositions. We analyse and synthesize whole-brain neuroimaging research, examined through functional magnetic imaging, to demonstrate that complexity science offers a principled and integrated approach to describe, and potentially understand, macroscale spontaneous brain functioning.

Robustness and Flexibility of Neural Function through Dynamical Criticality

In theoretical biology, robustness refers to the ability of a biological system to function properly even under perturbation of basic parameters (e.g., temperature or pH), which in mathematical models is reflected in not needing to fine-tune basic parameter constants; flexibility refers to the ability of a system to switch functions or behaviors easily and effortlessly. While there are extensive explorations of the concept of robustness and what it requires mathematically, understanding flexibility has proven more elusive, as well as also elucidating the apparent opposition between what is required mathematically for models to implement either. In this paper we address a number of arguments in theoretical neuroscience showing that both robustness and flexibility can be attained by systems that poise themselves at the onset of a large number of dynamical bifurcations, or dynamical criticality, and how such poising can have a profound influence on integration of information processing and function. Finally, we examine critical map lattices, which are coupled map lattices where the coupling is dynamically critical in the sense of having purely imaginary eigenvalues. We show that these map lattices provide an explicit connection between dynamical criticality in the sense we have used and "edge of chaos" criticality.

The Analysis of Mammalian Hearing Systems Supports the Hypothesis That Criticality Favors Neuronal Information Representation but Not Computation

In the neighborhood of critical states, distinct materials exhibit the same physical behavior, expressed by common simple laws among measurable observables, hence rendering a more detailed analysis of the individual systems obsolete. It is a widespread view that critical states are fundamental to neuroscience and directly favor computation. We argue here that from an evolutionary point of view, critical points seem indeed to be a natural phenomenon. Using mammalian hearing as our example, we show, however, explicitly that criticality does not describe the proper computational process and thus is only indirectly related to the computation in neural systems.

Scale-invariant representation of machine learning

The success of machine learning has resulted from its structured representation of data. Similar data have close internal representations as compressed codes for classification or emerged labels for clustering. We observe that the frequency of internal codes or labels follows power laws in both supervised and unsupervised learning models. This scale-invariant distribution implies that machine learning largely compresses frequent typical data, and simultaneously, differentiates many atypical data as outliers. In this study, we derive the process by which these power laws can naturally arise in machine learning. In terms of information theory, the scale-invariant representation corresponds to a maximally uncertain data grouping among possible representations that guarantee a given learning accuracy.

Neuronal Population Transitions Across a Quiescent-to-Active Frontier and Bifurcation

The mechanistic understanding of why neuronal population activity hovers on criticality remains unresolved despite the availability of experimental results. Without a coherent mathematical framework, the presence of power-law scaling is not straightforward to reconcile with findings implying epileptiform activity. Although multiple pictures have been proposed to relate the power-law scaling of avalanche statistics to phase transitions, the existence of a phase boundary in parameter space is until now an assumption. Herein, a framework based on differential inclusions, which departs from approaches constructed from differential equations, is shown to offer an adequate consolidation of evidences apparently connected to criticality and those linked to hyperexcitability. Through this framework, the phase boundary is elucidated in a parameter space spanned by variables representing levels of excitation and inhibition in a neuronal network. The interpretation of neuronal populations based on this approach offers insights on the role of pharmacological and endocrinal signaling in the homeostatic regulation of neuronal population activity.

On the Importance of Being Flexible: Dynamic Brain Networks and Their Potential Functional Significances

In this theoretical review, we begin by discussing brains and minds from a dynamical systems perspective, and then go on to describe methods for characterizing the flexibility of dynamic networks. We discuss how varying degrees and kinds of flexibility may be adaptive (or maladaptive) in different contexts, specifically focusing on measures related to either more disjoint or cohesive dynamics. While disjointed flexibility may be useful for assessing neural entropy, cohesive flexibility may potentially serve as a proxy for self-organized criticality as a fundamental property enabling adaptive behavior in complex systems. Particular attention is given to recent studies in which flexibility methods have been used to investigate neurological and cognitive maturation, as well as the breakdown of conscious processing under varying levels of anesthesia. We further discuss how these findings and methods might be contextualized within the Free Energy Principle with respect to the fundamentals of brain organization and biological functioning more generally, and describe potential methodological advances from this paradigm. Finally, with relevance to computational psychiatry, we propose a research program for obtaining a better understanding of ways that dynamic networks may relate to different forms of psychological flexibility, which may be the single most important factor for ensuring human flourishing.

Eco-epidemiological scaling of Leptospirosis: Vulnerability mapping and early warning forecasts

Infectious disease epidemics are plaguing the world and a lot of research is focused on the development of models to reproduce disease dynamics for eco-environmental and biological investigation, and disease management. Leptospirosis is an example of a neglected zoonosis strongly mediated by ecohydrological dynamics with emerging endemic and epidemic patterns worldwide in both animal and human populations. By accounting for large heterogeneities of affected areas we show how exponential endemics and scale-free epidemics are largely predictable and linked to common socio-environmental features via scaling laws with different exponents that inform about vulnerability factors. This led to the development of a novel pattern-oriented integrated model that can be used as an early-warning signal (EWS) tool for endemic-epidemic regime classification, risk determinant attribution, and near real-time forecast of outbreaks. Forecasts are grounded on expected outbreak recurrence time dependent on exceedance probabilities and statistical EWS that sense outbreak onset. A stochastic spatially-explicit model is shown to comprehensively predict outbreak dynamics (early sensing, timing, magnitude, decay, and eco-environmental determinants) and derive a spreading factor characterizing endemics and epidemics, where average over maximum rainfall is the critical factor characterizing disease transitions. Dynamically, case cross-correlation considering neighboring communities senses 2-weeks in advance outbreaks. Eco-environmental scaling relationships highlight how predicted host suitability and topographic index can be used as epidemiological footprints to effectively distinguish and control Leptospirosis regimes and areas dependent on hydro-climatological dynamics as the main trigger. The spatio-temporal scale-invariance of epidemics - underpinning persistent criticality and neutrality or independence among areas - is emphasized by the high accuracy in reproducing sequence and magnitude of cases via reliable surveillance. Further investigations of robustness and universality of eco-environmental determinants are required; nonetheless a comprehensive and computationally simple EWS method for the full characterization of Leptospirosis is provided. The tool is extendable to other climate-sensitive zoonoses to define vulnerability factors and predict outbreaks useful for optimal disease risk prevention and control.

Mutual Information Disentangles Interactions from Changing Environments

Real-world systems are characterized by complex interactions of their internal degrees of freedom, while living in ever-changing environments whose net effect is to act as additional couplings. Here, we introduce a paradigmatic interacting model in a switching, but unobserved, environment. We show that the limiting properties of the mutual information of the system allow for a disentangling of these two sources of couplings. Further, our approach might stand as a general method to discriminate complex internal interactions from equally complex changing environments.

Neuronal Avalanches Across the Rat Somatosensory Barrel Cortex and the Effect of Single Whisker Stimulation

Since its first experimental signatures, the so called "critical brain hypothesis" has been extensively studied. Yet, its actual foundations remain elusive. According to a widely accepted teleological reasoning, the brain would be poised to a critical state to optimize the mapping of the noisy and ever changing real-world inputs, thus suggesting that primary sensory cortical areas should be critical. We investigated whether a single barrel column of the somatosensory cortex of the anesthetized rat displays a critical behavior. Neuronal avalanches were recorded across all cortical layers in terms of both multi-unit activities and population local field potentials, and their behavior during spontaneous activity compared to the one evoked by a controlled single whisker deflection. By applying a maximum likelihood statistical method based on timeseries undersampling to fit the avalanches distributions, we show that neuronal avalanches are power law distributed for both multi-unit activities and local field potentials during spontaneous activity, with exponents that are spread along a scaling line. Instead, after the tactile stimulus, activity switches to a transient across-layers synchronization mode that appears to dominate the cortical representation of the single sensory input.

Optimal Input Representation in Neural Systems at the Edge of Chaos

Shedding light on how biological systems represent, process and store information in noisy environments is a key and challenging goal. A stimulating, though controversial, hypothesis poses that operating in dynamical regimes near the edge of a phase transition, i.e., at criticality or the "edge of chaos", can provide information-processing living systems with important operational advantages, creating, e.g., an optimal trade-off between robustness and flexibility. Here, we elaborate on a recent theoretical result, which establishes that the spectrum of covariance matrices of neural networks representing complex inputs in a robust way needs to decay as a power-law of the rank, with an exponent close to unity, a result that has been indeed experimentally verified in neurons of the mouse visual cortex. Aimed at understanding and mimicking these results, we construct an artificial neural network and train it to classify images. We find that the best performance in such a task is obtained when the network operates near the critical point, at which the eigenspectrum of the covariance matrix follows the very same statistics as actual neurons do. Thus, we conclude that operating near criticality can also have-besides the usually alleged virtues-the advantage of allowing for flexible, robust and efficient input representations.

Avalanches and edge-of-chaos learning in neuromorphic nanowire networks

The brain's efficient information processing is enabled by the interplay between its neuro-synaptic elements and complex network structure. This work reports on the neuromorphic dynamics of nanowire networks (NWNs), a unique brain-inspired system with synapse-like memristive junctions embedded within a recurrent neural network-like structure. Simulation and experiment elucidate how collective memristive switching gives rise to long-range transport pathways, drastically altering the network's global state via a discontinuous phase transition. The spatio-temporal properties of switching dynamics are found to be consistent with avalanches displaying power-law size and life-time distributions, with exponents obeying the crackling noise relationship, thus satisfying criteria for criticality, as observed in cortical neuronal cultures. Furthermore, NWNs adaptively respond to time varying stimuli, exhibiting diverse dynamics tunable from order to chaos. Dynamical states at the edge-of-chaos are found to optimise information processing for increasingly complex learning tasks. Overall, these results reveal a rich repertoire of emergent, collective neural-like dynamics in NWNs, thus demonstrating the potential for a neuromorphic advantage in information processing.

Signatures of brain criticality unveiled by maximum entropy analysis across cortical states

It has recently been reported that statistical signatures of brain criticality, obtained from distributions of neuronal avalanches, can depend on the cortical state. We revisit these claims with a completely different and independent approach, employing a maximum entropy model to test whether signatures of criticality appear in urethane-anesthetized rats. To account for the spontaneous variation of cortical states, we parse the time series and perform the maximum entropy analysis as a function of the variability of the population spiking activity. To compare data sets with different numbers of neurons, we define a normalized distance to criticality that takes into account the peak and width of the specific heat curve. We found a universal collapse of the normalized distance to criticality dependence on the cortical state, on an animal by animal basis. This indicates a universal dynamics and a critical point at an intermediate value of spiking variability.

Robust cortical criticality and diverse dynamics resulting from functional specification

Despite the recognition of the layered structure and evident criticality in the cortex, how the specification of input, output, and computational layers affects the self-organized criticality has not been much explored. By constructing heterogeneous structures with a well-accepted model of leaky neurons, we find that the specification can lead to robust criticality rather insensitive to the strength of external stimuli. This naturally unifies the adaptation to strong inputs without extra synaptic plasticity mechanisms. Low degree of recurrence constitutes an alternative explanation to subcriticality other than the high-frequency inputs. Unlike fully recurrent networks where external stimuli always render subcriticality, the dynamics of networks with sufficient feedforward connections can be driven to criticality and supercriticality. These findings indicate that functional and structural specification and their interplay with external stimuli are of crucial importance for the network dynamics. The robust criticality puts forward networks of the leaky neurons as promising platforms for realizing artificial neural networks that work in the vicinity of critical points.