Fractal complexity in spontaneous EEG metastable-state transitions: new vistas on integrated neural dynamics

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.

EEG Complexity Analysis of Brain States, Tasks and ASD Risk

Autism spectrum disorder is an increasingly prevalent and debilitating neurodevelopmental condition and an electroencephalogram (EEG) diagnostic challenge. Despite large amounts of electrophysiological research over many decades, an EEG biomarker for autism spectrum disorder (ASD) has not been found. We hypothesized that reductions in complex dynamical system behaviour in the human central nervous system as part of the macroscale neuronal function during cognitive processes might be detectable in whole EEG for higher-risk ASD adults. In three studies, we compared the medians of correlation dimension, largest Lyapunov exponent, Higuchi's fractal dimension, multiscale entropy, multifractal detrended fluctuation analysis and Kolmogorov complexity during resting, cognitive and social skill tasks in 20 EEG channels of 39 adults over a range of ASD risk. We found heterogeneous complexity distribution with clusters of hierarchical sequences pointing to potential cognitive processing differences, but no clear distinction based on ASD risk. We suggest that there is indication of statistically significant differences between complexity measures of brain states and tasks. Though replication of our studies is needed with a larger sample, we believe that our electrophysiological and analytic approach has potential as a biomarker for earlier ASD diagnosis.

The time scales of irreversibility in spontaneous brain activity are altered in obsessive compulsive disorder

We study how obsessive-compulsive disorder (OCD) affects the complexity and time-reversal symmetry-breaking (irreversibility) of the brain resting-state activity as measured by magnetoencephalography (MEG). Comparing MEG recordings from OCD patients and age/sex matched control subjects, we find that irreversibility is more concentrated at faster time scales and more uniformly distributed across different channels of the same hemisphere in OCD patients than in control subjects. Furthermore, the interhemispheric asymmetry between homologous areas of OCD patients and controls is also markedly different. Some of these differences were reduced by 1-year of Kundalini Yoga meditation treatment. Taken together, these results suggest that OCD alters the dynamic attractor of the brain's resting state and hint at a possible novel neurophysiological characterization of this psychiatric disorder and how this therapy can possibly modulate brain function.

Neural and behavioral adaptations to frontal theta neurofeedback training: A proof of concept study

Previous neurofeedback research has shown training-related frontal theta increases and performance improvements on some executive tasks in real feedback versus sham control groups. However, typical sham control groups receive false or non-contingent feedback, making it difficult to know whether observed differences between groups are associated with accurate contingent feedback or other cognitive mechanisms (motivation, control strategies, attentional engagement, fatigue, etc.). To address this question, we investigated differences between two frontal theta training groups, each receiving accurate contingent feedback, but with different top-down goals: (1) increase and (2) alternate increase/decrease. We hypothesized that the increase group would exhibit greater increases in frontal theta compared to the alternate group, which would exhibit lower frontal theta during down- versus up-modulation blocks over sessions. We also hypothesized that the alternate group would exhibit greater performance improvements on a Go-NoGo shooting task requiring alterations in behavioral activation and inhibition, as the alternate group would be trained with greater task specificity, suggesting that receiving accurate contingent feedback may be the more salient learning mechanism underlying frontal theta neurofeedback training gains. Thirty young healthy volunteers were randomly assigned to increase or alternate groups. Training consisted of an orientation session, five neurofeedback training sessions (six blocks of six 30-s trials of FCz theta modulation (4-7 Hz) separated by 10-s rest intervals), and six Go-NoGo testing sessions (four blocks of 90 trials in both Low and High time-stress conditions). Multilevel modeling revealed greater frontal theta increases in the alternate group over training sessions. Further, Go-NoGo task performance increased at a greater rate in the increase group (accuracy and reaction time, but not commission errors). Overall, these results reject our hypotheses and suggest that changes in frontal theta and performance outcomes were not explained by reinforcement learning afforded by accurate contingent feedback. We discuss our findings in terms of alternative conceptual and methodological considerations, as well as limitations of this research.

Attaining the recesses of the cognitive space

Existing neuropsychological tests of executive function often manifest a difficulty pinpointing cognitive deficits when these are intermittent and come in the form of omissions. We discuss the hypothesis that two partially interrelated reasons for this failure stem from relative inability of neuropsychological tests to explore the cognitive space and to explicitly take into account strategic and opportunistic resource allocation decisions, and to address the temporal aspects of both behaviour and task-related brain function in data analysis. Criteria for tasks suitable for neuropsychological assessment of executive function, as well as appropriate ways to analyse and interpret observed behavioural data are suggested. It is proposed that experimental tasks should be devised which emphasize typical rather than optimal performance, and that analyses should quantify path-dependent fluctuations in performance levels rather than averaged behaviour. Some implications for experimental neuropsychology are illustrated for the case of planning and problem-solving abilities and with particular reference to cognitive impairment in closed-head injury.

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.

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.

Functional Complex Networks Based on Operational Architectonics: Application on Electroencephalography-Brain-computer Interface for Imagined Speech

A new method for analyzing brain complex dynamics and states is presented. This method constructs functional brain graphs and is comprised of two pylons: (a) Operational architectonics (OA) concept of brain and mind functioning. (b) Network neuroscience. In particular, the algorithm utilizes OA framework for a non-parametric segmentation of EEGs, which leads to the identification of change points, namely abrupt jumps in EEG amplitude, called Rapid Transition Processes (RTPs). Subsequently, the time coordinates of RTPs are used for the generation of undirected weighted complex networks fulfilling a scale-free topology criterion, from which various network metrics of brain connectivity are estimated. These metrics form feature vectors, which can be used in machine learning algorithms for classification and/or prediction. The method is tested in classification problems on an EEG-based BCI data set, acquired from individuals during imagery pronunciation tasks of various words/vowels. The classification results, based on a Naïve Bayes classifier, show that the overall accuracies were found to be above chance level in all tested cases. This method was also compared with other state-of-the-art computational approaches commonly used for functional network generation, exhibiting competitive performance. The method can be useful to neuroscientists wishing to enhance their repository of brain research algorithms.

Characterization of the Functional Dynamics in the Neonatal Brain during REM and NREM Sleep States by means of Microstate Analysis

Neonates spend most of their life sleeping. During sleep, their brain experiences fast changes in its functional organization. Microstate analysis permits to capture the rapid dynamical changes occurring in the functional organization of the brain by representing the changing spatio-temporal features of the electroencephalogram (EEG) as a sequence of short-lasting scalp topographies—the microstates. In this study, we modeled the ongoing neonatal EEG into sequences of a limited number of microstates and investigated whether the extracted microstate features are altered in REM and NREM sleep (usually known as active and quiet sleep states—AS and QS—in the newborn) and depend on the EEG frequency band. 19-channel EEG recordings from 60 full-term healthy infants were analyzed using a modified version of the k-means clustering algorithm. The results show that ~ 70% of the variance in the datasets can be described using 7 dominant microstate templates. The mean duration and mean occurrence of the dominant microstates were significantly different in the two sleep states. Microstate syntax analysis demonstrated that the microstate sequences characterizing AS and QS had specific non-casual structures that differed in the two sleep states. Microstate analysis of the neonatal EEG in specific frequency bands showed a clear dependence of the explained variance on frequency. Overall, our findings demonstrate that (1) the spatio-temporal dynamics of the neonatal EEG can be described by non-casual sequences of a limited number of microstate templates; (2) the brain dynamics described by these microstate templates depends on frequency; (3) the features of the microstate sequences can well differentiate the physiological conditions characterizing AS and QS.

Neural field theory of neural avalanche exponents

The power-law exponents of observed size and lifetime distributions of near-critical neural avalanches are calculated from neural field theory using diagrammatic methods. This brings neural avalanches within the ambit of neural field theory, which has also previously explained near-critical 1/f spectra and many other observed features of neural activity. This strengthens the case for near-criticality of the brain and opens the way for these other phenomena to be interrelated with avalanches and their dynamics.

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

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

Power spectrum of resting-state blood-oxygen-level-dependent signal

Hemodynamic modeling is used to explore the origin, predict, and analyze the power spectrum of the resting-state blood-oxygen-level-dependent (BOLD) signal measured by functional magnetic resonance imaging (fMRI), which has been reported to have a power-law form, i.e., P(f)∝f^{-s}, where P(f) is the power, f is the frequency, and s>0 is the power-law exponent. However, current fMRI experimental paradigms have limited acquisition durations, affecting the spectral resolution of fMRI data at the low-frequency regime. Here, the claimed power-law spectrum is investigated by using a recent hemodynamic model to analytically derive the BOLD power spectrum, with parameters that are related to neurophysiology. The theoretical results show that, for all realistic parameter combinations, the BOLD power spectrum is flat at f≲0.01Hz, has a weak resonance originating from intrinsic oscillations of vasodilatory response, and becomes a power law for high frequencies, all of which is in agreement with an empirical data set that describes the spectrum of one subject and brain region. However, the results are contrary to studies reporting a pure power-law spectrum at f≲0.2Hz. The discrepancy is attributed largely to data averaging employed by current approaches that averages together important properties of the BOLD power spectrum, such as its resonance, that biases the spectrum to only show a power law. Data averaging also reduces the high-frequency power-law exponent relative to individual cases. Overall, this work demonstrates how the model can reproduce BOLD dynamics and further analyze its low-frequency behavior. Moreover, it also uses the model to explain the impact of procedures, such as data averaging, on the reported features of the BOLD power spectrum.

The avalanche-like behaviour of large-scale haemodynamic activity from wakefulness to deep sleep

Increasing evidence suggests that responsiveness is associated with critical or near-critical cortical dynamics, which exhibit scale-free cascades of spatio-temporal activity. These cascades, or 'avalanches', have been detected at multiple scales, from in vitro and in vivo microcircuits to voltage imaging and brain-wide functional magnetic resonance imaging (fMRI) recordings. Criticality endows the cortex with certain information-processing capacities postulated as necessary for conscious wakefulness, yet it remains unknown how unresponsiveness impacts on the avalanche-like behaviour of large-scale human haemodynamic activity. We observed a scale-free hierarchy of co-activated connected clusters by applying a point-process transformation to fMRI data recorded during wakefulness and non-rapid eye movement (NREM) sleep. Maximum-likelihood estimates revealed a significant effect of sleep stage on the scaling parameters of the cluster size power-law distributions. Post hoc statistical tests showed that differences were maximal between wakefulness and N2 sleep. These results were robust against spatial coarse graining, fitting alternative statistical models and different point-process thresholds, and disappeared upon phase shuffling the fMRI time series. Evoked neural bistabilities preventing arousals during N2 sleep do not suffice to explain these differences, which point towards changes in the intrinsic dynamics of the brain that could be necessary to consolidate a state of deep unresponsiveness.

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

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

Gauging Functional Brain Activity: From Distinguishability to Accessibility

Standard neuroimaging techniques provide non-invasive access not only to human brain anatomy but also to its physiology. The activity recorded with these techniques is generally called functional imaging, but what is observed per se is an instance of dynamics, from which functional brain activity should be extracted. Distinguishing between bare dynamics and genuine function is a highly non-trivial task, but a crucially important one when comparing experimental observations and interpreting their significance. Here we illustrate how neuroimaging's ability to extract genuine functional brain activity is bounded by functional representations' structure. To do so, we first provide a simple definition of functional brain activity from a system-level brain imaging perspective. We then review how the properties of the space on which brain activity is represented induce relations on observed imaging data which allow determining the extent to which two observations are functionally distinguishable and quantifying how far apart they are. It is also proposed that genuine functional distances would require defining accessibility, i.e., how a given observed condition can be accessed from another given one, under the dynamics of some neurophysiological process. We show how these properties result from the structure defined on dynamical data and dynamics-to-function projections, and consider some implications that the way and extent to which these are defined have for the interpretation of experimental data from standard system-level brain recording techniques.

Neurofeedback: Principles, appraisal, and outstanding issues

Neurofeedback is a form of brain training in which subjects are fed back information about some measure of their brain activity which they are instructed to modify in a way thought to be functionally advantageous. Over the last 20 years, neurofeedback has been used to treat various neurological and psychiatric conditions, and to improve cognitive function in various contexts. However, in spite of a growing popularity, neurofeedback protocols typically make (often covert) assumptions on what aspects of brain activity to target, where in the brain to act and how, which have far-reaching implications for the assessment of its potential and efficacy. Here we critically examine some conceptual and methodological issues associated with the way neurofeedback's general objectives and neural targets are defined. The neural mechanisms through which neurofeedback may act at various spatial and temporal scales, and the way its efficacy is appraised are reviewed, and the extent to which neurofeedback may be used to control functional brain activity discussed. Finally, it is proposed that gauging neurofeedback's potential, as well as assessing and improving its efficacy will require better understanding of various fundamental aspects of brain dynamics and a more precise definition of functional brain activity and brain-behaviour relationships.

K-th Nearest Neighbor (KNN) Entropy Estimates of Complexity and Integration from Ongoing and Stimulus-Evoked Electroencephalographic (EEG) Recordings of the Human Brain

Information-theoretic measures for quantifying multivariate statistical dependence have proven useful for the study of the unity and diversity of the human brain. Two such measures–integration, I(X), and interaction complexity, C_I(X)–have been previously applied to electroencephalographic (EEG) signals recorded during ongoing wakeful brain states. Here, I(X) and C_I(X) were computed for empirical and simulated visually-elicited alpha-range (8–13 Hz) EEG signals. Integration and complexity of evoked (stimulus-locked) and induced (non-stimulus-locked) EEG responses were assessed using nonparametric k-th nearest neighbor (KNN) entropy estimation, which is robust to the nonstationarity of stimulus-elicited EEG signals. KNN-based I(X) and C_I(X) were also computed for the alpha-range EEG of ongoing wakeful brain states. I(X) and C_I(X) patterns differentiated between induced and evoked EEG signals and replicated previous wakeful EEG findings obtained using Gaussian-based entropy estimators. Absolute levels of I(X) and C_I(X) were related to absolute levels of alpha-range EEG power and phase synchronization, but stimulus-related changes in the information-theoretic and other EEG properties were independent. These findings support the hypothesis that visual perception and ongoing wakeful mental states emerge from complex, dynamical interaction among segregated and integrated brain networks operating near an optimal balance between order and disorder.

Bridging Waves and Crucial Events in the Dynamics of the Brain

Earlier research work on the dynamics of the brain, disclosing the existence of crucial events, is revisited for the purpose of making the action of crucial events, responsible for the 1/f −noise in the brain, compatible with the wave-like nature of the brain processes. We review the relevant neurophysiological literature to make clear that crucial events are generated by criticality. We also show that although criticality generates a strong deviation from the regular wave-like behavior, under the form of Rapid Transition Processes, the brain dynamics also host crucial events in regions of nearly coherent oscillations, thereby making many crucial events virtually invisible. Furthermore, the anomalous scaling generated by the crucial events can be established with high accuracy by means of direct analysis of raw data, suggested by a theoretical perspective not requiring the crucial events to yield a visible physical effect. The latter follows from the fact that periodicity, waves and crucial events are the consequences of a spontaneous process of self-organization. We obtain three main results: (a) the important role of crucial events is confirmed and established with greater accuracy than previously; (b) we demonstrate the theoretical tools necessary to understand the joint action of crucial events and periodicity; (c) we argue that the results of this paper can be used to shed light on the nature of this important process of self-organization, thereby contributing to the understanding of cognition.

Negative mood and mind wandering increase long-range temporal correlations in attention fluctuations

There is growing evidence that the intermittent nature of mind wandering episodes and mood have a pronounced influence on trial-to-trial variability in performance. Nevertheless, the temporal dynamics and significance of such lapses in attention remains inadequately understood. Here, we hypothesize that the dynamics of fluctuations in sustained attention between external and internal sources of information obey so-called critical-state dynamics, characterized by trial-to-trial dependencies with long-range temporal correlations. To test this, we performed behavioral investigations measuring reaction times in a visual sustained attention task and cued introspection in probe-caught reports of mind wandering. We show that trial-to-trial variability in reaction times exhibit long-range temporal correlations in agreement with the criticality hypothesis. Interestingly, we observed the fastest responses in subjects with the weakest long-range temporal correlations and show the vital effect of mind wandering and bad mood on this response variability. The implications of these results stress the importance of future research to increase focus on behavioral variability.

Neuronal avalanche dynamics indicates different universality classes in neuronal cultures

Neuronal avalanches have become an ubiquitous tool to describe the activity of large neuronal assemblies. The emergence of scale-free statistics with well-defined exponents has led to the belief that the brain might operate near a critical point. Yet not much is known in terms of how the different exponents arise or how robust they are. Using calcium imaging recordings of dissociated neuronal cultures we show that the exponents are not universal, and that significantly different exponents arise with different culture preparations, leading to the existence of different universality classes. Naturally developing cultures show avalanche statistics consistent with those of a mean-field branching process, however, cultures grown in the presence of folic acid metabolites appear to be in a distinct universality class with significantly different critical exponents. Given the increased synaptic density and number of feedback loops in folate reared cultures, our results suggest that network topology plays a leading role in shaping the avalanche dynamics. We also show that for both types of cultures pronounced correlations exist in the sizes of neuronal avalanches indicating size clustering, being much stronger in folate reared cultures.

Connectome-harmonic decomposition of human brain activity reveals dynamical repertoire re-organization under LSD

Recent studies have started to elucidate the effects of lysergic acid diethylamide (LSD) on the human brain but the underlying dynamics are not yet fully understood. Here we used 'connectome-harmonic decomposition', a novel method to investigate the dynamical changes in brain states. We found that LSD alters the energy and the power of individual harmonic brain states in a frequency-selective manner. Remarkably, this leads to an expansion of the repertoire of active brain states, suggestive of a general re-organization of brain dynamics given the non-random increase in co-activation across frequencies. Interestingly, the frequency distribution of the active repertoire of brain states under LSD closely follows power-laws indicating a re-organization of the dynamics at the edge of criticality. Beyond the present findings, these methods open up for a better understanding of the complex brain dynamics in health and disease.

Topodynamics of metastable brains

The brain displays both the anatomical features of a vast amount of interconnected topological mappings as well as the functional features of a nonlinear, metastable system at the edge of chaos, equipped with a phase space where mental random walks tend towards lower energetic basins. Nevertheless, with the exception of some advanced neuro-anatomic descriptions and present-day connectomic research, very few studies have been addressing the topological path of a brain embedded or embodied in its external and internal environment. Herein, by using new formal tools derived from algebraic topology, we provide an account of the metastable brain, based on the neuro-scientific model of Operational Architectonics of brain-mind functioning. We introduce a "topodynamic" description that shows how the relationships among the countless intertwined spatio-temporal levels of brain functioning can be assessed in terms of projections and mappings that take place on abstract structures, equipped with different dimensions, curvatures and energetic constraints. Such a topodynamical approach, apart from providing a biologically plausible model of brain function that can be operationalized, is also able to tackle the issue of a long-standing dichotomy: it throws indeed a bridge between the subjective, immediate datum of the naïve complex of sensations and mentations and the objective, quantitative, data extracted from experimental neuro-scientific procedures. Importantly, it opens the door to a series of new predictions and future directions of advancement for neuroscientific research.

Self-organized dynamical complexity in human wakefulness and sleep: different critical brain-activity feedback for conscious and unconscious states

Criticality reportedly describes brain dynamics. The main critical feature is the presence of scale-free neural avalanches, whose auto-organization is determined by a critical branching ratio of neural-excitation spreading. Other features, directly associated to second-order phase transitions, are: (i) scale-free-network topology of functional connectivity, stemming from suprathreshold pairwise correlations, superimposable, in waking brain activity, with that of ferromagnets at Curie temperature; (ii) temporal long-range memory associated to renewal intermittency driven by abrupt fluctuations in the order parameters, detectable in human brain via spatially distributed phase or amplitude changes in EEG activity. Herein we study intermittent events, extracted from 29 night EEG recordings, including presleep wakefulness and all phases of sleep, where different levels of mentation and consciousness are present. We show that while critical avalanching is unchanged, at least qualitatively, intermittency and functional connectivity, present during conscious phases (wakefulness and REM sleep), break down during both shallow and deep non-REM sleep. We provide a theory for fragmentation-induced intermittency breakdown and suggest that the main difference between conscious and unconscious states resides in the backwards causation, namely on the constraints that the emerging properties at large scale induce to the lower scales. In particular, while in conscious states this backwards causation induces a critical slowing down, preserving spatiotemporal correlations, in dreamless sleep we see a self-organized maintenance of moduli working in parallel. Critical avalanches are still present, and establish transient auto-organization, whose enhanced fluctuations are able to trigger sleep-protecting mechanisms that reinstate parallel activity. The plausible role of critical avalanches in dreamless sleep is to provide a rapid recovery of consciousness, if stimuli are highly arousing.

How can we study reasoning in the brain?

The brain did not develop a dedicated device for reasoning. This fact bears dramatic consequences. While for perceptuo-motor functions neural activity is shaped by the input's statistical properties, and processing is carried out at high speed in hardwired spatially segregated modules, in reasoning, neural activity is driven by internal dynamics and processing times, stages, and functional brain geometry are largely unconstrained a priori. Here, it is shown that the complex properties of spontaneous activity, which can be ignored in a short-lived event-related world, become prominent at the long time scales of certain forms of reasoning. It is argued that the neural correlates of reasoning should in fact be defined in terms of non-trivial generic properties of spontaneous brain activity, and that this implies resorting to concepts, analytical tools, and ways of designing experiments that are as yet non-standard in cognitive neuroscience. The implications in terms of models of brain activity, shape of the neural correlates, methods of data analysis, observability of the phenomenon, and experimental designs are discussed.

Relationship of fast- and slow-timescale neuronal dynamics in human MEG and SEEG

A growing body of evidence suggests that the neuronal dynamics are poised at criticality. Neuronal avalanches and long-range temporal correlations (LRTCs) are hallmarks of such critical dynamics in neuronal activity and occur at fast (subsecond) and slow (seconds to hours) timescales, respectively. The critical dynamics at different timescales can be characterized by their power-law scaling exponents. However, insight into the avalanche dynamics and LRTCs in the human brain has been largely obtained with sensor-level MEG and EEG recordings, which yield only limited anatomical insight and results confounded by signal mixing. We investigated here the relationship between the human neuronal dynamics at fast and slow timescales using both source-reconstructed MEG and intracranial stereotactical electroencephalography (SEEG). Both MEG and SEEG revealed avalanche dynamics that were characterized parameter-dependently by power-law or truncated-power-law size distributions. Both methods also revealed robust LRTCs throughout the neocortex with distinct scaling exponents in different functional brain systems and frequency bands. The exponents of power-law regimen neuronal avalanches and LRTCs were strongly correlated across subjects. Qualitatively similar power-law correlations were also observed in surrogate data without spatial correlations but with scaling exponents distinct from those of original data. Furthermore, we found that LRTCs in the autonomous nervous system, as indexed by heart-rate variability, were correlated in a complex manner with cortical neuronal avalanches and LRTCs in MEG but not SEEG. These scalp and intracranial data hence show that power-law scaling behavior is a pervasive but neuroanatomically inhomogeneous property of neuronal dynamics in central and autonomous nervous systems.

Self-organized criticality as a fundamental property of neural systems

The neural criticality hypothesis states that the brain may be poised in a critical state at a boundary between different types of dynamics. Theoretical and experimental studies show that critical systems often exhibit optimal computational properties, suggesting the possibility that criticality has been evolutionarily selected as a useful trait for our nervous system. Evidence for criticality has been found in cell cultures, brain slices, and anesthetized animals. Yet, inconsistent results were reported for recordings in awake animals and humans, and current results point to open questions about the exact nature and mechanism of criticality, as well as its functional role. Therefore, the criticality hypothesis has remained a controversial proposition. Here, we provide an account of the mathematical and physical foundations of criticality. In the light of this conceptual framework, we then review and discuss recent experimental studies with the aim of identifying important next steps to be taken and connections to other fields that should be explored.

Functional significance of complex fluctuations in brain activity: from resting state to cognitive neuroscience

Behavioral studies have shown that human cognition is characterized by properties such as temporal scale invariance, heavy-tailed non-Gaussian distributions, and long-range correlations at long time scales, suggesting models of how (non observable) components of cognition interact. On the other hand, results from functional neuroimaging studies show that complex scaling and intermittency may be generic spatio-temporal properties of the brain at rest. Somehow surprisingly, though, hardly ever have the neural correlates of cognition been studied at time scales comparable to those at which cognition shows scaling properties. Here, we analyze the meanings of scaling properties and the significance of their task-related modulations for cognitive neuroscience. It is proposed that cognitive processes can be framed in terms of complex generic properties of brain activity at rest and, ultimately, of functional equations, limiting distributions, symmetries, and possibly universality classes characterizing them.

Universal organization of resting brain activity at the thermodynamic critical point

Thermodynamic criticality describes emergent phenomena in a wide variety of complex systems. In the mammalian cortex, one type of complex dynamics that spontaneously emerges from neuronal interactions has been characterized as neuronal avalanches. Several aspects of neuronal avalanches such as their size and life time distributions are described by power laws with unique exponents, indicating an underlying critical branching process that governs avalanche formation. Here, we show that neuronal avalanches also reflect an organization of brain dynamics close to a thermodynamic critical point. We recorded spontaneous cortical activity in monkeys and humans at rest using high-density intracranial microelectrode arrays and magnetoencephalography, respectively. By numerically changing a control parameter equivalent to thermodynamic temperature, we observed typical critical behavior in cortical activities near the actual physiological condition, including the phase transition of an order parameter, as well as the divergence of susceptibility and specific heat. Finite-size scaling of these quantities allowed us to derive robust critical exponents highly consistent across monkey and humans that uncover a distinct, yet universal organization of brain dynamics. Our results demonstrate that normal brain dynamics at rest resides near or at criticality, which maximizes several aspects of information processing such as input sensitivity and dynamic range.

Time scales in cognitive neuroscience

Cognitive neuroscience boils down to describing the ways in which cognitive function results from brain activity. In turn, brain activity shows complex fluctuations, with structure at many spatio-temporal scales. Exactly how cognitive function inherits the physical dimensions of neural activity, though, is highly non-trivial, and so are generally the corresponding dimensions of cognitive phenomena. As for any physical phenomenon, when studying cognitive function, the first conceptual step should be that of establishing its dimensions. Here, we provide a systematic presentation of the temporal aspects of task-related brain activity, from the smallest scale of the brain imaging technique's resolution, to the observation time of a given experiment, through the characteristic time scales of the process under study. We first review some standard assumptions on the temporal scales of cognitive function. In spite of their general use, these assumptions hold true to a high degree of approximation for many cognitive (viz. fast perceptual) processes, but have their limitations for other ones (e.g., thinking or reasoning). We define in a rigorous way the temporal quantifiers of cognition at all scales, and illustrate how they qualitatively vary as a function of the properties of the cognitive process under study. We propose that each phenomenon should be approached with its own set of theoretical, methodological and analytical tools. In particular, we show that when treating cognitive processes such as thinking or reasoning, complex properties of ongoing brain activity, which can be drastically simplified when considering fast (e.g., perceptual) processes, start playing a major role, and not only characterize the temporal properties of task-related brain activity, but also determine the conditions for proper observation of the phenomena. Finally, some implications on the design of experiments, data analyses, and the choice of recording parameters are discussed.

Neuronal avalanches in the resting MEG of the human brain

What constitutes normal cortical dynamics in healthy human subjects is a major question in systems neuroscience. Numerous in vitro and in vivo animal studies have shown that ongoing or resting cortical dynamics are characterized by cascades of activity across many spatial scales, termed neuronal avalanches. In experiment and theory, avalanche dynamics are identified by two measures: (1) a power law in the size distribution of activity cascades with an exponent of -3/2 and (2) a branching parameter of the critical value of 1, reflecting balanced propagation of activity at the border of premature termination and potential blowup. Here we analyzed resting-state brain activity recorded using noninvasive magnetoencephalography (MEG) from 124 healthy human subjects and two different MEG facilities using different sensor technologies. We identified large deflections at single MEG sensors and combined them into spatiotemporal cascades on the sensor array using multiple timescales. Cascade size distributions obeyed power laws. For the timescale at which the branching parameter was close to 1, the power law exponent was -3/2. This relationship was robust to scaling and coarse graining of the sensor array. It was absent in phase-shuffled controls with the same power spectrum or empty scanner data. Our results demonstrate that normal cortical activity in healthy human subjects at rest organizes as neuronal avalanches and is well described by a critical branching process. Theory and experiment have shown that such critical, scale-free dynamics optimize information processing. Therefore, our findings imply that the human brain attains an optimal dynamical regime for information processing.

The role of the circadian system in fractal neurophysiological control

Many neurophysiological variables such as heart rate, motor activity, and neural activity are known to exhibit intrinsic fractal fluctuations - similar temporal fluctuation patterns at different time scales. These fractal patterns contain information about health, as many pathological conditions are accompanied by their alteration or absence. In physical systems, such fluctuations are characteristic of critical states on the border between randomness and order, frequently arising from nonlinear feedback interactions between mechanisms operating on multiple scales. Thus, the existence of fractal fluctuations in physiology challenges traditional conceptions of health and disease, suggesting that high levels of integrity and adaptability are marked by complex variability, not constancy, and are properties of a neurophysiological network, not individual components. Despite the subject's theoretical and clinical interest, the neurophysiological mechanisms underlying fractal regulation remain largely unknown. The recent discovery that the circadian pacemaker (suprachiasmatic nucleus) plays a crucial role in generating fractal patterns in motor activity and heart rate sheds an entirely new light on both fractal control networks and the function of this master circadian clock, and builds a bridge between the fields of circadian biology and fractal physiology. In this review, we sketch the emerging picture of the developing interdisciplinary field of fractal neurophysiology by examining the circadian system's role in fractal regulation.

From self-organized to extended criticality

We address the issue of criticality that is attracting the attention of an increasing number of neurophysiologists. Our main purpose is to establish the specific nature of some dynamical processes that although physically different, are usually termed as “critical,” and we focus on those characterized by the cooperative interaction of many units. We notice that the term “criticality” has been adopted to denote both noise-induced phase transitions and Self-Organized Criticality (SOC) with no clear connection with the traditional phase transitions, namely the transformation of a thermodynamic system from one state of matter to another. We notice the recent attractive proposal of extended criticality advocated by Bailly and Longo, which is realized through a wide set of critical points rather than emerging as a singularity from a unique value of the control parameter. We study a set of cooperatively firing neurons and we show that for an extended set of interaction couplings the system exhibits a form of temporal complexity similar to that emerging at criticality from ordinary phase transitions. This extended criticality regime is characterized by three main properties: (i) In the ideal limiting case of infinitely large time period, temporal complexity corresponds to Mittag-Leffler complexity; (ii) For large values of the interaction coupling the periodic nature of the process becomes predominant while maintaining to some extent, in the intermediate time asymptotic region, the signature of complexity; (iii) Focusing our attention on firing neuron avalanches, we find two of the popular SOC properties, namely the power indexes 2 and 1.5 respectively for time length and for the intensity of the avalanches. We derive the main conclusion that SOC emerges from extended criticality, thereby explaining the experimental observation of Plenz and Beggs: avalanches occur in time with surprisingly regularity, in apparent conflict with the temporal complexity of physical critical points.

Cooperation-induced topological complexity: a promising road to fault tolerance and hebbian learning

According to an increasing number of researchers intelligence emerges from criticality as a consequence of locality breakdown and long-range correlation, well known properties of phase transition processes. We study a model of interacting units, as an idealization of real cooperative systems such as the brain or a flock of birds, for the purpose of discussing the emergence of long-range correlation from the coupling of any unit with its nearest neighbors. We focus on the critical condition that has been recently shown to maximize information transport and we study the topological structure of the network of dynamically linked nodes. Although the topology of this network depends on the arbitrary choice of correlation threshold, namely the correlation intensity selected to establish a link between two nodes; the numerical calculations of this paper afford some important indications on the dynamically induced topology. The first important property is the emergence of a perception length as large as the flock size, thanks to some nodes with a large number of links, thus playing the leadership role. All the units are equivalent and leadership moves in time from one to another set of nodes, thereby insuring fault tolerance. Then we focus on the correlation threshold generating a scale-free topology with power index ν ≈ 1 and we find that if this topological structure is selected to establish consensus through the linked nodes, the control parameter necessary to generate criticality is close to the critical value corresponding to the all-to-all coupling condition. We find that criticality in this case generates also a third state, corresponding to a total lack of consensus. However, we make a numerical analysis of the dynamically induced network, and we find that it consists of two almost independent structures, each of which is equivalent to a network in the all-to-all coupling condition. This observation confirms that cooperation makes the system evolve toward favoring consensus topological structures. We argue that these results are compatible with both Hebbian learning and fault tolerance.

Criticality in large-scale brain FMRI dynamics unveiled by a novel point process analysis

Functional magnetic resonance imaging (fMRI) techniques have contributed significantly to our understanding of brain function. Current methods are based on the analysis of gradual and continuous changes in the brain blood oxygenated level dependent (BOLD) signal. Departing from that approach, recent work has shown that equivalent results can be obtained by inspecting only the relatively large amplitude BOLD signal peaks, suggesting that relevant information can be condensed in discrete events. This idea is further explored here to demonstrate how brain dynamics at resting state can be captured just by the timing and location of such events, i.e., in terms of a spatiotemporal point process. The method allows, for the first time, to define a theoretical framework in terms of an order and control parameter derived from fMRI data, where the dynamical regime can be interpreted as one corresponding to a system close to the critical point of a second order phase transition. The analysis demonstrates that the resting brain spends most of the time near the critical point of such transition and exhibits avalanches of activity ruled by the same dynamical and statistical properties described previously for neuronal events at smaller scales. Given the demonstrated functional relevance of the resting state brain dynamics, its representation as a discrete process might facilitate large-scale analysis of brain function both in health and disease.

Letting the brain speak for itself

Metaphors of Computation and Information tended to detract attention from the intrinsic modes of neural system functions, uncontaminated by the observer's role in collection, and interpretation of experimental data. Recognizing the self-referential mode of function, and the propensity for self-organization to critical states requires a fundamentally new orientation, based on Complex System Dynamics as non-ergodic, non-stationary processes with inverse-power-law statistical distributions. Accordingly, local cooperative processes, intrinsic to neural structures, and of fractal nature, call for applying Fractional Calculus and models of Random Walks with long-term memory in Theoretical Neuroscience studies.