How critical is brain criticality?

Adolescent-to-adult gains in cognitive flexibility are adaptively supported by reward sensitivity, exploration, and neural variability

Cognitive flexibility exhibits dynamic changes throughout development, with different forms of flexibility showing dissociable developmental trajectories. In this review, we propose that an adolescent-specific mode of flexibility in the face of changing environmental contingencies supports the emergence of adolescent-to-adult gains in cognitive shifting efficiency. We first describe how cognitive shifting abilities monotonically improve from childhood to adulthood, accompanied by increases in brain state flexibility, neural variability, and excitatory/inhibitory balance. We next summarize evidence supporting the existence of a dopamine-driven, adolescent peak in flexible behavior that results in reward seeking, undirected exploration, and environmental sampling. We propose a neurodevelopmental framework that relates these adolescent behaviors to the refinement of neural phenotypes relevant to mature cognitive flexibility, and thus highlight the importance of the adolescent period in fostering healthy neurocognitive trajectories.

Optimal input reverberation and homeostatic self-organization toward the edge of synchronization

Transient or partial synchronization can be used to do computations, although a fully synchronized network is sometimes related to the onset of epileptic seizures. Here, we propose a homeostatic mechanism that is capable of maintaining a neuronal network at the edge of a synchronization transition, thereby avoiding the harmful consequences of a fully synchronized network. We model neurons by maps since they are dynamically richer than integrate-and-fire models and more computationally efficient than conductance-based approaches. We first describe the synchronization phase transition of a dense network of neurons with different tonic spiking frequencies coupled by gap junctions. We show that at the transition critical point, inputs optimally reverberate through the network activity through transient synchronization. Then, we introduce a local homeostatic dynamic in the synaptic coupling and show that it produces a robust self-organization toward the edge of this phase transition. We discuss the potential biological consequences of this self-organization process, such as its relation to the Brain Criticality hypothesis, its input processing capacity, and how its malfunction could lead to pathological synchronization and the onset of seizure-like activity.

Low-dimensional criticality embedded in high-dimensional awake brain dynamics

Whether cortical neurons operate in a strongly or weakly correlated dynamical regime determines fundamental information processing capabilities and has fueled decades of debate. We offer a resolution of this debate; we show that two important dynamical regimes, typically considered incompatible, can coexist in the same local cortical circuit by separating them into two different subspaces. In awake mouse motor cortex, we find a low-dimensional subspace with large fluctuations consistent with criticality-a dynamical regime with moderate correlations and multi-scale information capacity and transmission. Orthogonal to this critical subspace, we find a high-dimensional subspace containing a desynchronized dynamical regime, which may optimize input discrimination. The critical subspace is apparent only at long timescales, which explains discrepancies among some previous studies. Using a computational model, we show that the emergence of a low-dimensional critical subspace at large timescales agrees with established theory of critical dynamics. Our results suggest that the cortex leverages its high dimensionality to multiplex dynamical regimes across different subspaces.

Directed and acyclic synaptic connectivity in the human layer 2-3 cortical microcircuit

The computational capabilities of neuronal networks are fundamentally constrained by their specific connectivity. Previous studies of cortical connectivity have mostly been carried out in rodents; whether the principles established therein also apply to the evolutionarily expanded human cortex is unclear. We studied network properties within the human temporal cortex using samples obtained from brain surgery. We analyzed multineuron patch-clamp recordings in layer 2-3 pyramidal neurons and identified substantial differences compared with rodents. Reciprocity showed random distribution, synaptic strength was independent from connection probability, and connectivity of the supragranular temporal cortex followed a directed and mostly acyclic graph topology. Application of these principles in neuronal models increased dimensionality of network dynamics, suggesting a critical role for cortical computation.

Local signal variability and functional connectivity: Sensitive measures of the excitation-inhibition ratio?

A novel network version of permutation entropy, the inverted joint permutation entropy (JPEinv), holds potential as non-invasive biomarker of abnormal excitation-inhibition (E-I) ratio in Alzheimer's disease (AD). In this computational modelling study, we test the hypotheses that this metric, and related measures of signal variability and functional connectivity, are sensitive to altered E-I ratios. The E-I ratio in each neural mass of a whole-brain computational network model was systematically varied. We evaluated whether JPEinv, local signal variability (by permutation entropy) and functional connectivity (by weighted symbolic mutual information (wsMI)) were related to E-I ratio, on whole-brain and regional level. The hub disruption index can identify regions primarily affected in terms of functional connectivity strength (or: degree) by the altered E-I ratios. Analyses were performed for a range of coupling strengths, filter and time-delay settings. On whole-brain level, higher E-I ratios were associated with higher functional connectivity (by JPEinv and wsMI) and lower local signal variability. These relationships were nonlinear and depended on the coupling strength, filter and time-delay settings. On regional level, hub-like regions showed a selective decrease in functional degree (by JPEinv and wsMI) upon a lower E-I ratio, and non-hub-like regions showed a selective increase in degree upon a higher E-I ratio. These results suggest that abnormal functional connectivity and signal variability, as previously reported in patients across the AD continuum, can inform us about altered E-I ratios.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11571-023-10003-x.

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.

Fractional order memcapacitive neuromorphic elements reproduce and predict neuronal function

There is an increasing need to implement neuromorphic systems that are both energetically and computationally efficient. There is also great interest in using electric elements with memory, memelements, that can implement complex neuronal functions intrinsically. A feature not widely incorporated in neuromorphic systems is history-dependent action potential time adaptation which is widely seen in real cells. Previous theoretical work shows that power-law history dependent spike time adaptation, seen in several brain areas and species, can be modeled with fractional order differential equations. Here, we show that fractional order spiking neurons can be implemented using super-capacitors. The super-capacitors have fractional order derivative and memcapacitive properties. We implemented two circuits, a leaky integrate and fire and a Hodgkin-Huxley. Both circuits show power-law spiking time adaptation and optimal coding properties. The spiking dynamics reproduced previously published computer simulations. However, the fractional order Hodgkin-Huxley circuit showed novel dynamics consistent with criticality. We compared the responses of this circuit to recordings from neurons in the weakly-electric fish that have previously been shown to perform fractional order differentiation of their sensory input. The criticality seen in the circuit was confirmed in spontaneous recordings in the live fish. Furthermore, the circuit also predicted long-lasting stimulation that was also corroborated experimentally. Our work shows that fractional order memcapacitors provide intrinsic memory dependence that could allow implementation of computationally efficient neuromorphic devices. Memcapacitors are static elements that consume less energy than the most widely studied memristors, thus allowing the realization of energetically efficient neuromorphic devices.

Does the brain behave like a (complex) network? I. Dynamics

Graph theory is now becoming a standard tool in system-level neuroscience. However, endowing observed brain anatomy and dynamics with a complex network structure does not entail that the brain actually works as a network. Asking whether the brain behaves as a network means asking whether network properties count. From the viewpoint of neurophysiology and, possibly, of brain physics, the most substantial issues a network structure may be instrumental in addressing relate to the influence of network properties on brain dynamics and to whether these properties ultimately explain some aspects of brain function. Here, we address the dynamical implications of complex network, examining which aspects and scales of brain activity may be understood to genuinely behave as a network. To do so, we first define the meaning of networkness, and analyse some of its implications. We then examine ways in which brain anatomy and dynamics can be endowed with a network structure and discuss possible ways in which network structure may be shown to represent a genuine organisational principle of brain activity, rather than just a convenient description of its anatomy and dynamics.

Landau model for illustrating the learning and unlearning process of nociplastic pain

Recent advancements in understanding the consolidation of nociplastic pain point to a complex, non-conscious learned process of threat perception. Neurobiological pain education is emerging as a promising approach to unlearn nociplastic pain, supported by biopsychosocial tools such as exposure to movement, mindfulness, and group sharing formats. However, this approach is still not well-known among clinicians and the society at large, creating a communication problem that unfortunately perpetuates the suffering of patients. Herein, we propose a Landau model to describe the learning and unlearning process of nociplastic pain, aiming to clarify this complex situation and facilitate communication across different sectors of the society. Nociplastic pain corresponds to a first-order transition, with attention more likely in the alert-protection state than in the trust-explore state. Two appealing results of the model are that the perception of the critical context depends on personal history regarding the symptom and that biopsychosocial loops are formed when there is alarming learned historical information about the symptom, along with confused and contradictory expert information, as seen in nocebo messages. Learning and unlearning in the model correspond to a chang in control parametrs that can weigh more on the alert-protection state, trust-explore state, uncertain state or neutral state. This description clarifies why neurobiological education is the foundational therapy from which others must be built to embody the accessible, clear, and trustworthy information.

A bifurcation integrates information from many noisy ion channels and allows for milli-Kelvin thermal sensitivity in the snake pit organ

In various biological systems, information from many noisy molecular receptors must be integrated into a collective response. A striking example is the thermal imaging organ of pit vipers. Single nerve fibers in the organ reliably respond to milli-Kelvin (mK) temperature increases, a thousand times more sensitive than their molecular sensors, thermo-transient receptor potential (TRP) ion channels. Here, we propose a mechanism for the integration of this molecular information. In our model, amplification arises due to proximity to a dynamical bifurcation, separating a regime with frequent and regular action potentials (APs), from a regime where APs are irregular and infrequent. Near the transition, AP frequency can have an extremely sharp dependence on temperature, naturally accounting for the thousand-fold amplification. Furthermore, close to the bifurcation, most of the information about temperature available in the TRP channels' kinetics can be read out from the times between consecutive APs even in the presence of readout noise. A key model prediction is that the coefficient of variation in the distribution of interspike times decreases with AP frequency, and quantitative comparison with experiments indeed suggests that nerve fibers of snakes are located very close to the bifurcation. While proximity to such bifurcation points typically requires fine-tuning of parameters, we propose that having feedback act from the order parameter (AP frequency) onto the control parameter robustly maintains the system in the vicinity of the bifurcation. This robustness suggests that similar feedback mechanisms might be found in other sensory systems which also need to detect tiny signals in a varying environment.

Mindfulness meditation styles differently modulate source-level MEG microstate dynamics and complexity

Background

The investigation of mindfulness meditation practice, classically divided into focused attention meditation (FAM), and open monitoring meditation (OMM) styles, has seen a long tradition of theoretical, affective, neurophysiological and clinical studies. In particular, the high temporal resolution of magnetoencephalography (MEG) or electroencephalography (EEG) has been exploited to fill the gap between the personal experience of meditation practice and its neural correlates. Mounting evidence, in fact, shows that human brain activity is highly dynamic, transiting between different brain states (microstates). In this study, we aimed at exploring MEG microstates at source-level during FAM, OMM and in the resting state, as well as the complexity and criticality of dynamic transitions between microstates.

Methods

Ten right-handed Theravada Buddhist monks with a meditative expertise of minimum 2,265 h participated in the experiment. MEG data were acquired during a randomized block design task (6 min FAM, 6 min OMM, with each meditative block preceded and followed by 3 min resting state). Source reconstruction was performed using eLORETA on individual cortical space, and then parcellated according to the Human Connect Project atlas. Microstate analysis was then applied to parcel level signals in order to derive microstate topographies and indices. In addition, from microstate sequences, the Hurst exponent and the Lempel-Ziv complexity (LZC) were computed.

Results

Our results show that the coverage and occurrence of specific microstates are modulated either by being in a meditative state or by performing a specific meditation style. Hurst exponent values in both meditation conditions are reduced with respect to the value observed during rest, LZC shows significant differences between OMM, FAM, and REST, with a progressive increase from REST to FAM to OMM.

Discussion

Importantly, we report changes in brain criticality indices during meditation and between meditation styles, in line with a state-like effect of meditation on cognitive performance. In line with previous reports, we suggest that the change in cognitive state experienced in meditation is paralleled by a shift with respect to critical points in brain dynamics.

Reconstruction of phase dynamics from macroscopic observations based on linear and nonlinear response theories

We propose a method to reconstruct the phase dynamics in rhythmical interacting systems from macroscopic responses to weak inputs by developing linear and nonlinear response theories, which predict the responses in a given system. By solving an inverse problem, the method infers an unknown system: the natural frequency distribution, the coupling function, and the time delay which is inevitable in real systems. In contrast to previous methods, our method requires neither strong invasiveness nor microscopic observations. We demonstrate that the method reconstructs two phase systems from observed responses accurately. The qualitative methodological advantages demonstrated by our quantitative numerical examinations suggest its broad applicability in various fields, including brain systems, which are often observed through macroscopic signals such as electroencephalograms and functional magnetic response imaging.

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.

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.

Connectome-based reservoir computing with the conn2res toolbox

The connection patterns of neural circuits form a complex network. How signaling in these circuits manifests as complex cognition and adaptive behaviour remains the central question in neuroscience. Concomitant advances in connectomics and artificial intelligence open fundamentally new opportunities to understand how connection patterns shape computational capacity in biological brain networks. Reservoir computing is a versatile paradigm that uses high-dimensional, nonlinear dynamical systems to perform computations and approximate cognitive functions. Here we present conn2res: an open-source Python toolbox for implementing biological neural networks as artificial neural networks. conn2res is modular, allowing arbitrary network architecture and dynamics to be imposed. The toolbox allows researchers to input connectomes reconstructed using multiple techniques, from tract tracing to noninvasive diffusion imaging, and to impose multiple dynamical systems, from spiking neurons to memristive dynamics. The versatility of the conn2res toolbox allows us to ask new questions at the confluence of neuroscience and artificial intelligence. By reconceptualizing function as computation, conn2res sets the stage for a more mechanistic understanding of structure-function relationships in brain networks.

Dynamical models reveal anatomically reliable attractor landscapes embedded in resting state brain networks

Analyses of functional connectivity (FC) in resting-state brain networks (RSNs) have generated many insights into cognition. However, the mechanistic underpinnings of FC and RSNs are still not well-understood. It remains debated whether resting state activity is best characterized as noise-driven fluctuations around a single stable state, or instead, as a nonlinear dynamical system with nontrivial attractors embedded in the RSNs. Here, we provide evidence for the latter, by constructing whole-brain dynamical systems models from individual resting-state fMRI (rfMRI) recordings, using the Mesoscale Individualized NeuroDynamic (MINDy) platform. The MINDy models consist of hundreds of neural masses representing brain parcels, connected by fully trainable, individualized weights. We found that our models manifested a diverse taxonomy of nontrivial attractor landscapes including multiple equilibria and limit cycles. However, when projected into anatomical space, these attractors mapped onto a limited set of canonical RSNs, including the default mode network (DMN) and frontoparietal control network (FPN), which were reliable at the individual level. Further, by creating convex combinations of models, bifurcations were induced that recapitulated the full spectrum of dynamics found via fitting. These findings suggest that the resting brain traverses a diverse set of dynamics, which generates several distinct but anatomically overlapping attractor landscapes. Treating rfMRI as a unimodal stationary process (i.e., conventional FC) may miss critical attractor properties and structure within the resting brain. Instead, these may be better captured through neural dynamical modeling and analytic approaches. The results provide new insights into the generative mechanisms and intrinsic spatiotemporal organization of brain networks.

Criticality supports cross-frequency cortical-thalamic information transfer during conscious states

Consciousness is thought to be regulated by bidirectional information transfer between the cortex and thalamus, but the nature of this bidirectional communication - and its possible disruption in unconsciousness - remains poorly understood. Here, we present two main findings elucidating mechanisms of corticothalamic information transfer during conscious states. First, we identify a highly preserved spectral channel of cortical-thalamic communication that is present during conscious states, but which is diminished during the loss of consciousness and enhanced during psychedelic states. Specifically, we show that in humans, mice, and rats, information sent from either the cortex or thalamus via δ/θ/α waves (∼1-13 Hz) is consistently encoded by the other brain region by high γ waves (52-104 Hz); moreover, unconsciousness induced by propofol anesthesia or generalized spike-and-wave seizures diminishes this cross-frequency communication, whereas the psychedelic 5-methoxy-N,N-dimethyltryptamine (5-MeO-DMT) enhances this low-to-high frequency interregional communication. Second, we leverage numerical simulations and neural electrophysiology recordings from the thalamus and cortex of human patients, rats, and mice to show that these changes in cross-frequency cortical-thalamic information transfer may be mediated by excursions of low-frequency thalamocortical electrodynamics toward/away from edge-of-chaos criticality, or the phase transition from stability to chaos. Overall, our findings link thalamic-cortical communication to consciousness, and further offer a novel, mathematically well-defined framework to explain the disruption to thalamic-cortical information transfer during unconscious states.

Protosequences in human cortical organoids model intrinsic states in the developing cortex

Neuronal firing sequences are thought to be the basic building blocks of neural coding and information broadcasting within the brain. However, when sequences emerge during neurodevelopment remains unknown. We demonstrate that structured firing sequences are present in spontaneous activity of human brain organoids and ex vivo neonatal brain slices from the murine somatosensory cortex. We observed a balance between temporally rigid and flexible firing patterns that are emergent phenomena in human brain organoids and early postnatal murine somatosensory cortex, but not in primary dissociated cortical cultures. Our findings suggest that temporal sequences do not arise in an experience-dependent manner, but are rather constrained by an innate preconfigured architecture established during neurogenesis. These findings highlight the potential for brain organoids to further explore how exogenous inputs can be used to refine neuronal circuits and enable new studies into the genetic mechanisms that govern assembly of functional circuitry during early human brain development.

Altered spread of waves of activities at large scale is influenced by cortical thickness organization in temporal lobe epilepsy: a magnetic resonance imaging-high-density electroencephalography study

Temporal lobe epilepsy is a brain network disorder characterized by alterations at both the structural and the functional levels. It remains unclear how structure and function are related and whether this has any clinical relevance. In the present work, we adopted a novel methodological approach investigating how network structural features influence the large-scale dynamics. The functional network was defined by the spatio-temporal spreading of aperiodic bursts of activations (neuronal avalanches), as observed utilizing high-density electroencephalography in patients with temporal lobe epilepsy. The structural network was modelled as the region-based thickness covariance. Loosely speaking, we quantified the similarity of the cortical thickness of any two brain regions, both across groups and at the individual level, the latter utilizing a novel approach to define the subject-wise structural covariance network. In order to compare the structural and functional networks (at the nodal level), we studied the correlation between the probability that a wave of activity would propagate from a source to a target region and the similarity of the source region thickness as compared with other target brain regions. Building on the recent evidence that large-waves of activities pathologically spread through the epileptogenic network in temporal lobe epilepsy, also during resting state, we hypothesize that the structural cortical organization might influence such altered spatio-temporal dynamics. We observed a stable cluster of structure-function correlation in the bilateral limbic areas across subjects, highlighting group-specific features for left, right and bilateral temporal epilepsy. The involvement of contralateral areas was observed in unilateral temporal lobe epilepsy. We showed that in temporal lobe epilepsy, alterations of structural and functional networks pair in the regions where seizures propagate and are linked to disease severity. In this study, we leveraged on a well-defined model of neurological disease and pushed forward personalization approaches potentially useful in clinical practice. Finally, the methods developed here could be exploited to investigate the relationship between structure-function networks at subject level in other neurological conditions.

Cortical ripples mediate top-down suppression of hippocampal reactivation during sleep memory consolidation

Consolidation of initially encoded hippocampal representations in the neocortex through reactivation is crucial for long-term memory formation, and is facilitated by the coordination of hippocampal sharp-wave ripples (SWRs) with cortical oscillations during non-REM sleep. However, the contribution of high-frequency cortical ripples to consolidation is still unclear. We used continuous recordings in the hippocampus and prefrontal cortex (PFC) over the course of spatial learning and show that independent PFC ripples, when dissociated from SWRs, predominantly suppress hippocampal activity in non-REM sleep. PFC ripples paradoxically mediate top-down suppression of hippocampal reactivation, which is inversely related to reactivation strength during coordinated CA1-PFC ripples. Further, we show non-canonical, serial coordination of ripples with cortical slow and spindle oscillations. These results establish a role for cortical ripples in regulating consolidation.

Metabolically regulated spiking could serve neuronal energy homeostasis and protect from reactive oxygen species

So-called spontaneous activity is a central hallmark of most nervous systems. Such non-causal firing is contrary to the tenet of spikes as a means of communication, and its purpose remains unclear. We propose that self-initiated firing can serve as a release valve to protect neurons from the toxic conditions arising in mitochondria from lower-than-baseline energy consumption. To demonstrate the viability of our hypothesis, we built a set of models that incorporate recent experimental results indicating homeostatic control of metabolic products-Adenosine triphosphate (ATP), adenosine diphosphate (ADP), and reactive oxygen species (ROS)-by changes in firing. We explore the relationship of metabolic cost of spiking with its effect on the temporal patterning of spikes and reproduce experimentally observed changes in intrinsic firing in the fruitfly dorsal fan-shaped body neuron in a model with ROS-modulated potassium channels. We also show that metabolic spiking homeostasis can produce indefinitely sustained avalanche dynamics in cortical circuits. Our theory can account for key features of neuronal activity observed in many studies ranging from ion channel function all the way to resting state dynamics. We finish with a set of experimental predictions that would confirm an integrated, crucial role for metabolically regulated spiking and firmly link metabolic homeostasis and neuronal function.

Computation via Neuron-like Spiking in Percolating Networks of Nanoparticles

The biological brain is a highly efficient computational system in which information processing is performed via electrical spikes. Neuromorphic computing systems that work on similar principles could support the development of the next generation of artificial intelligence and, in particular, enable low-power edge computing. Percolating networks of nanoparticles (PNNs) have previously been shown to exhibit critical spiking behavior, with promise for highly efficient natural computation. Here we employ a rate coding scheme to show that PNNs can perform Boolean operations and image classification. Near perfect accuracy is achieved in both tasks by manipulating the spiking activity using certain control voltages. We demonstrate that the key to successful computation is that nanoscale tunnel gaps within the percolating networks transform input data through a powerful modulus-like nonlinearity. These results provide a basis for implementation of further computational schemes that exploit the brain-like criticality of these networks.

Criticality of resting-state EEG predicts perturbational complexity and level of consciousness during anesthesia

Consciousness has been proposed to be supported by electrophysiological patterns poised at criticality, a dynamical regime which exhibits adaptive computational properties, maximally complex patterns and divergent sensitivity to perturbation. Here, we investigated dynamical properties of the resting-state electroencephalogram of healthy subjects undergoing general anesthesia with propofol, xenon or ketamine. We then studied the relation of these dynamic properties with the perturbational complexity index (PCI), which has shown remarkably high sensitivity in detecting consciousness independent of behavior. All participants were unresponsive under anesthesia, while consciousness was retained only during ketamine anesthesia (in the form of vivid dreams)., enabling an experimental dissociation between unresponsiveness and unconsciousness. We estimated (i) avalanche criticality, (ii) chaoticity, and (iii) criticality-related measures, and found that states of unconsciousness were characterized by a distancing from both the edge of activity propagation and the edge of chaos. We were then able to predict individual subjects' PCI (i.e., PCImax) with a mean absolute error below 7%. Our results establish a firm link between the PCI and criticality and provide further evidence for the role of criticality in the emergence of consciousness.

Criticality of neuronal avalanches in human sleep and their relationship with sleep macro- and micro-architecture

Sleep plays a key role in preserving brain function, keeping brain networks in a state that ensures optimal computation. Empirical evidence indicates that this state is consistent with criticality, where scale-free neuronal avalanches emerge. However, the connection between sleep architecture and brain tuning to criticality remains poorly understood. Here, we characterize the critical behavior of avalanches and study their relationship with sleep macro- and micro-architectures, in particular, the cyclic alternating pattern (CAP). We show that avalanches exhibit robust scaling behaviors, with exponents obeying scaling relations consistent with the mean-field directed percolation universality class. We demonstrate that avalanche dynamics is modulated by the NREM-REM cycles and that, within NREM sleep, avalanche occurrence correlates with CAP activation phases-indicating a potential link between CAP and brain tuning to criticality. The results open new perspectives on the collective dynamics underlying CAP function, and on the relationship between sleep architecture, avalanches, and self-organization to criticality.

Improvements in task performance after practice are associated with scale-free dynamics of brain activity

Although practicing a task generally benefits later performance on that same task, there are individual differences in practice effects. One avenue to model such differences comes from research showing that brain networks extract functional advantages from operating in the vicinity of criticality, a state in which brain network activity is more scale-free. We hypothesized that higher scale-free signal from fMRI data, measured with the Hurst exponent (H), indicates closer proximity to critical states. We tested whether individuals with higher H during repeated task performance would show greater practice effects. In Study 1, participants performed a dual-n-back task (DNB) twice during MRI (n = 56). In Study 2, we used two runs of n-back task (NBK) data from the Human Connectome Project sample (n = 599). In Study 3, participants performed a word completion task (CAST) across six runs (n = 44). In all three studies, multivariate analysis was used to test whether higher H was related to greater practice-related performance improvement. Supporting our hypothesis, we found patterns of higher H that reliably correlated with greater performance improvement across participants in all three studies. However, the predictive brain regions were distinct, suggesting that the specific spatial H↑ patterns are not task-general.

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.

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

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

Beyond rhythm - a framework for understanding the frequency spectrum of neural activity

Cognitive and behavioral processes are often accompanied by changes within well-defined frequency bands of the local field potential (LFP i.e., the voltage induced by neuronal activity). These changes are detectable in the frequency domain using the Fourier transform and are often interpreted as neuronal oscillations. However, aside some well-known exceptions, the processes underlying such changes are difficult to track in time, making their oscillatory nature hard to verify. In addition, many non-periodic neural processes can also have spectra that emphasize specific frequencies. Thus, the notion that spectral changes reflect oscillations is likely too restrictive. In this study, we use a simple yet versatile framework to understand the frequency spectra of neural recordings. Using simulations, we derive the Fourier spectra of periodic, quasi-periodic and non-periodic neural processes having diverse waveforms, illustrating how these attributes shape their spectral signatures. We then show how neural processes sum their energy in the local field potential in simulated and real-world recording scenarios. We find that the spectral power of neural processes is essentially determined by two aspects: (1) the distribution of neural events in time and (2) the waveform of the voltage induced by single neural events. Taken together, this work guides the interpretation of the Fourier spectrum of neural recordings and indicates that power increases in specific frequency bands do not necessarily reflect periodic neural activity.

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.

The 1/f-like behavior of neural field spectra are a natural consequence of noise driven brain dynamics

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

Brain criticality predicts individual levels of inter-areal synchronization in human electrophysiological data

Neuronal oscillations and their synchronization between brain areas are fundamental for healthy brain function. Yet, synchronization levels exhibit large inter-individual variability that is associated with behavioral variability. We test whether individual synchronization levels are predicted by individual brain states along an extended regime of critical-like dynamics - the Griffiths phase (GP). We use computational modelling to assess how synchronization is dependent on brain criticality indexed by long-range temporal correlations (LRTCs). We analyze LRTCs and synchronization of oscillations from resting-state magnetoencephalography and stereo-electroencephalography data. Synchronization and LRTCs are both positively linearly and quadratically correlated among healthy subjects, while in epileptogenic areas they are negatively linearly correlated. These results show that variability in synchronization levels is explained by the individual position along the GP with healthy brain areas operating in its subcritical and epileptogenic areas in its supercritical side. We suggest that the GP is fundamental for brain function allowing individual variability while retaining functional advantages of criticality.

Event-related variability is modulated by task and development

In carefully designed experimental paradigms, cognitive scientists interpret the mean event-related potentials (ERP) in terms of cognitive operations. However, the huge signal variability from one trial to the next, questions the representability of such mean events. We explored here whether this variability is an unwanted noise, or an informative part of the neural response. We took advantage of the rapid changes in the visual system during human infancy and analyzed the variability of visual responses to central and lateralized faces in 2-to 6-month-old infants compared to adults using high-density electroencephalography (EEG). We observed that neural trajectories of individual trials always remain very far from ERP components, only moderately bending their direction with a substantial temporal jitter across trials. However, single trial trajectories displayed characteristic patterns of acceleration and deceleration when approaching ERP components, as if they were under the active influence of steering forces causing transient attraction and stabilization. These dynamic events could only partly be accounted for by induced microstate transitions or phase reset phenomena. Importantly, these structured modulations of response variability, both between and within trials, had a rich sequential organization, which in infants, was modulated by the task difficulty and age. Our approaches to characterize Event Related Variability (ERV) expand on classic ERP analyses and provide the first evidence for the functional role of ongoing neural variability in human infants.

Low dimensional criticality embedded in high dimensional awake brain dynamics

Whether cortical neurons operate in a strongly or weakly correlated dynamical regime determines fundamental information processing capabilities and has fueled decades of debate. Here we offer a resolution of this debate; we show that two important dynamical regimes, typically considered incompatible, can coexist in the same local cortical circuit by separating them into two different subspaces. In awake mouse motor cortex, we find a low-dimensional subspace with large fluctuations consistent with criticality - a dynamical regime with moderate correlations and multi-scale information capacity and transmission. Orthogonal to this critical subspace, we find a high-dimensional subspace containing a desynchronized dynamical regime, which may optimize input discrimination. The critical subspace is apparent only at long timescales, which explains discrepancies among some previous studies. Using a computational model, we show that the emergence of a low-dimensional critical subspace at large timescale agrees with established theory of critical dynamics. Our results suggest that cortex leverages its high dimensionality to multiplex dynamical regimes across different subspaces.

Biases in inverse Ising estimates of near-critical behavior

Inverse Ising inference allows pairwise interactions of complex binary systems to be reconstructed from empirical correlations. Typical estimators used for this inference, such as pseudo-likelihood maximization (PLM), are biased. Using the Sherrington-Kirkpatrick model as a benchmark, we show that these biases are large in critical regimes close to phase boundaries, and they may alter the qualitative interpretation of the inferred model. In particular, we show that the small-sample bias causes models inferred through PLM to appear closer to criticality than one would expect from the data. Data-driven methods to correct this bias are explored and applied to a functional magnetic resonance imaging data set from neuroscience. Our results indicate that additional care should be taken when attributing criticality to real-world data sets.

Distance to Criticality Undergoes Critical Transition Before Epileptic Seizure Attacks

Epilepsy is a common neurological disorder characterized by recurring seizures, but its underlying mechanisms remain poorly understood. Despite extensive research, there are still gaps in our knowledge about the relationship between brain dynamics and seizures. In this study, our aim is to address these gaps by proposing a novel approach to assess the role of brain network dynamics in the onset of seizures. Specifically, we investigate the relationship between brain dynamics and seizures by tracking the distance to criticality. Our hypothesis is that this distance plays a crucial role in brain state changes and that seizures may be related to critical transitions of this distance. To test this hypothesis, we develop a method to measure the evolution of the brain network's distance to the critical dynamic systems (i.e., the distance to the tipping point, DTP) using dynamic network biomarker theory and random matrix theory. The results show that the DTP of the brain decreases significantly immediately after onset of an epileptic seizure, suggesting that the brain loses its well-defined quasi-critical state during seizures. We refer to this phenomenon as the "criticality of the criticality" (COC). Furthermore, we observe that DTP exhibits a shape transition before and after the onset of the seizures. This phenomenon suggests the possibility of early warning signal (EWS) identification in the dynamic sequence of DTP, which could be utilized for seizure prediction. Our results show that the Hurst exponent, skewness, kurtosis, autocorrelation, and variance of the DTP sequence are potential EWS features. This study advances our understanding of the relationship between brain dynamics and seizures and highlights the potential for using criticality-based measures to predict and prevent seizures.

Bifurcation in space: Emergence of function modularity in the neocortex

How does functional modularity emerge in a multiregional cortex made with repeats of a canonical local circuit architecture? We investigated this question by focusing on neural coding of working memory, a core cognitive function. Here we report a mechanism dubbed "bifurcation in space", and show that its salient signature is spatially localized "critical slowing down" leading to an inverted V-shaped profile of neuronal time constants along the cortical hierarchy during working memory. The phenomenon is confirmed in connectome-based large-scale models of mouse and monkey cortices, offering an experimentally testable prediction to assess whether working memory representation is modular. Many bifurcations in space could explain the emergence of different activity patterns potentially deployed for distinct cognitive functions, This work demonstrates that a distributed mental representation is compatible with functional specificity as a consequence of macroscopic gradients of neurobiological properties across the cortex, suggesting a general principle for understanding brain's modular organization.

A scoping review for building a criticality-based conceptual framework of altered states of consciousness

The healthy conscious brain is thought to operate near a critical state, reflecting optimal information processing and high susceptibility to external stimuli. Conversely, deviations from the critical state are hypothesized to give rise to altered states of consciousness (ASC). Measures of criticality could therefore be an effective way of establishing the conscious state of an individual. Furthermore, characterizing the direction of a deviation from criticality may enable the development of treatment strategies for pathological ASC. The aim of this scoping review is to assess the current evidence supporting the criticality hypothesis, and the use of criticality as a conceptual framework for ASC. Using the PRISMA guidelines, Web of Science and PubMed were searched from inception to February 7th 2022 to find articles relating to measures of criticality across ASC. N = 427 independent papers were initially found on the subject. N = 378 were excluded because they were either: not related to criticality; not related to consciousness; not presenting results from a primary study; presenting model data. N = 49 independent papers were included in the present research, separated in 7 sub-categories of ASC: disorders of consciousness (DOC) (n = 5); sleep (n = 13); anesthesia (n = 18); epilepsy (n = 12); psychedelics and shamanic state of consciousness (n = 4); delirium (n = 1); meditative state (n = 2). Each category included articles suggesting a deviation of the critical state. While most studies were only able to identify a deviation from criticality without being certain of its direction, the preliminary consensus arising from the literature is that non-rapid eye movement (NREM) sleep reflects a subcritical state, epileptic seizures reflect a supercritical state, and psychedelics are closer to the critical state than normal consciousness. This scoping review suggests that, though the literature is limited and methodologically inhomogeneous, ASC are characterized by a deviation from criticality, though its direction is not clearly reported in a majority of studies. Criticality could become, with more extensive research, an effective and objective way to characterize ASC, and help identify therapeutic avenues to improve criticality in pathological brain states. Furthermore, we suggest how anesthesia and psychedelics could potentially be used as neuromodulation techniques to restore criticality in DOC.

Pupil size predicts the onset of exploration in brain and behavior

In uncertain environments, intelligent decision-makers exploit actions that have been rewarding in the past, but also explore actions that could be even better. Several neuromodulatory systems are implicated in exploration, based, in part, on work linking exploration to pupil size-a peripheral correlate of neuromodulatory tone and index of arousal. However, pupil size could instead track variables that make exploration more likely, like volatility or reward, without directly predicting either exploration or its neural bases. Here, we simultaneously measured pupil size, exploration, and neural population activity in the prefrontal cortex while two rhesus macaques explored and exploited in a dynamic environment. We found that pupil size under constant luminance specifically predicted the onset of exploration, beyond what could be explained by reward history. Pupil size also predicted disorganized patterns of prefrontal neural activity at both the single neuron and population levels, even within periods of exploitation. Ultimately, our results support a model in which pupil-linked mechanisms promote the onset of exploration via driving the prefrontal cortex through a critical tipping point where prefrontal control dynamics become disorganized and exploratory decisions are possible.

Aperiodic brain activity and response to anesthesia vary in disorders of consciousness

In the human electroencephalogram (EEG), oscillatory power peaks co-exist with non-oscillatory, aperiodic activity. Although EEG analysis has traditionally focused exclusively on oscillatory power, recent investigations have shown that the aperiodic EEG component can distinguish conscious wakefulness from sleep and anesthetic-induced unconsciousness. This study investigates the aperiodic EEG component of individuals in a disorder of consciousness (DOC); how it changes in response to exposure to anesthesia; and how it relates to the brain's information richness and criticality. High-density EEG was recorded from 43 individuals in a DOC, with 16 of these individuals undergoing a protocol of propofol anesthesia. The aperiodic component was defined by the spectral slope of the power spectral density. Our results demonstrate that the EEG aperiodic component is more informative about the participants' level of consciousness than the oscillatory component, especially for patients that suffered from a stroke. Importantly, the pharmacologically induced change in the spectral slope from 30-45 Hz positively correlated with individual's pre-anesthetic level of consciousness. The pharmacologically induced loss of information-richness and criticality was associated with individual's pre-anesthetic aperiodic component. During exposure to anesthesia, the aperiodic component was correlated with 3-month recovery status for individuals with DOC. The aperiodic EEG component has been historically neglected; this research highlights the necessity of considering this measure for the assessment of individuals in DOC and future research that seeks to understand the neurophysiological underpinnings of consciousness.

Resting-state oscillations reveal disturbed excitation-inhibition ratio in Alzheimer's disease patients

An early disruption of neuronal excitation-inhibition (E-I) balance in preclinical animal models of Alzheimer's disease (AD) has been frequently reported, but is difficult to measure directly and non-invasively in humans. Here, we examined known and novel neurophysiological measures sensitive to E-I in patients across the AD continuum. Resting-state magnetoencephalography (MEG) data of 86 amyloid-biomarker-confirmed subjects across the AD continuum (17 patients diagnosed with subjective cognitive decline, 18 with mild cognitive impairment (MCI) and 51 with dementia due to probable AD (AD dementia)), 46 healthy elderly and 20 young control subjects were reconstructed to source-space. E-I balance was investigated by detrended fluctuation analysis (DFA), a functional E/I (fE/I) algorithm, and the aperiodic exponent of the power spectrum. We found a disrupted E-I ratio in AD dementia patients specifically, by a lower DFA, and a shift towards higher excitation, by a higher fE/I and a lower aperiodic exponent. Healthy subjects showed lower fE/I ratios (< 1.0) than reported in previous literature, not explained by age or choice of an arbitrary threshold parameter, which warrants caution in interpretation of fE/I results. Correlation analyses showed that a lower DFA (E-I imbalance) and a lower aperiodic exponent (more excitation) was associated with a worse cognitive score in AD dementia patients. In contrast, a higher DFA in the hippocampi of MCI patients was associated with a worse cognitive score. This MEG-study showed E-I imbalance, likely due to increased excitation, in AD dementia, but not in early stage AD patients. To accurately determine the direction of shift in E-I balance, validations of the currently used markers and additional in vivo markers of E-I are required.

Structure and function in artificial, zebrafish and human neural networks

Network science provides a set of tools for the characterization of the structure and functional behavior of complex systems. Yet a major problem is to quantify how the structural domain is related to the dynamical one. In other words, how the diversity of dynamical states of a system can be predicted from the static network structure? Or the reverse problem: starting from a set of signals derived from experimental recordings, how can one discover the network connections or the causal relations behind the observed dynamics? Despite the advances achieved over the last two decades, many challenges remain concerning the study of the structure-dynamics interplay of complex systems. In neuroscience, progress is typically constrained by the low spatio-temporal resolution of experiments and by the lack of a universal inferring framework for empirical systems. To address these issues, applications of network science and artificial intelligence to neural data have been rapidly growing. In this article, we review important recent applications of methods from those fields to the study of the interplay between structure and functional dynamics of human and zebrafish brain. We cover the selection of topological features for the characterization of brain networks, inference of functional connections, dynamical modeling, and close with applications to both the human and zebrafish brain. This review is intended to neuroscientists who want to become acquainted with techniques from network science, as well as to researchers from the latter field who are interested in exploring novel application scenarios in neuroscience.

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.

Bistability and criticality in the stochastic Wilson-Cowan model

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

Altered spreading of neuronal avalanches in temporal lobe epilepsy relates to cognitive performance: A resting-state hdEEG study

OBJECTIVE

Large aperiodic bursts of activations named neuronal avalanches have been used to characterize whole-brain activity, as their presence typically relates to optimal dynamics. Epilepsy is characterized by alterations in large-scale brain network dynamics. Here we exploited neuronal avalanches to characterize differences in electroencephalography (EEG) basal activity, free from seizures and/or interictal spikes, between patients with temporal lobe epilepsy (TLE) and matched controls.

METHOD

We defined neuronal avalanches as starting when the z-scored source-reconstructed EEG signals crossed a specific threshold in any region and ending when all regions returned to baseline. This technique avoids data manipulation or assumptions of signal stationarity, focusing on the aperiodic, scale-free components of the signals. We computed individual avalanche transition matrices to track the probability of avalanche spreading across any two regions, compared them between patients and controls, and related them to memory performance in patients.

RESULTS

We observed a robust topography of significant edges clustering in regions functionally and structurally relevant for the TLE, such as the entorhinal cortex, the inferior parietal and fusiform area, the inferior temporal gyrus, and the anterior cingulate cortex. We detected a significant correlation between the centrality of the entorhinal cortex in the transition matrix and the long-term memory performance (delay recall Rey-Osterrieth Complex Figure Test).

SIGNIFICANCE

Our results show that the propagation patterns of large-scale neuronal avalanches are altered in TLE during the resting state, suggesting a potential diagnostic application in epilepsy. Furthermore, the relationship between specific patterns of propagation and memory performance support the neurophysiological relevance of neuronal avalanches.

Spontaneous neuronal avalanches as a correlate of access consciousness

Decades of research have advanced our understanding of the biophysical mechanisms underlying consciousness. However, an overarching framework bridging between models of consciousness and the large-scale organization of spontaneous brain activity is still missing. Based on the observation that spontaneous brain activity dynamically switches between epochs of segregation and large-scale integration of information, we hypothesize a brain-state dependence of conscious access, whereby the presence of either segregated or integrated states marks distinct modes of information processing. We first review influential works on the neuronal correlates of consciousness, spontaneous resting-state brain activity and dynamical system theory. Then, we propose a test experiment to validate our hypothesis that conscious access occurs in aperiodic cycles, alternating windows where new incoming information is collected but not experienced, to punctuated short-lived integration events, where conscious access to previously collected content occurs. In particular, we suggest that the integration events correspond to neuronal avalanches, which are collective bursts of neuronal activity ubiquitously observed in electrophysiological recordings. If confirmed, the proposed framework would link the physics of spontaneous cortical dynamics, to the concept of ignition within the global neuronal workspace theory, whereby conscious access manifest itself as a burst of neuronal activity.