Being critical of criticality in the brain

Scale-Change Symmetry in the Rules Governing Neural Systems

Similar universal phenomena can emerge in different complex systems when those systems share a common symmetry in their governing laws. In physical systems operating near a critical phase transition, the governing physical laws obey a fractal symmetry; they are the same whether considered at fine or coarse scales. This scale-change symmetry is responsible for universal critical phenomena found across diverse systems. Experiments suggest that the cerebral cortex can also operate near a critical phase transition. Thus we hypothesize that the laws governing cortical dynamics may obey scale-change symmetry. Here we develop a practical approach to test this hypothesis. We confirm, using two different computational models, that neural dynamical laws exhibit scale-change symmetry near a dynamical phase transition. Moreover, we show that as a mouse awakens from anesthesia, scale-change symmetry emerges. Scale-change symmetry of the rules governing cortical dynamics may explain observations of similar critical phenomena across diverse neural systems.

Time-delay model of perceptual decision making in cortical networks

It is known that cortical networks operate on the edge of instability, in which oscillations can appear. However, the influence of this dynamic regime on performance in decision making, is not well understood. In this work, we propose a population model of decision making based on a winner-take-all mechanism. Using this model, we demonstrate that local slow inhibition within the competing neuronal populations can lead to Hopf bifurcation. At the edge of instability, the system exhibits ambiguity in the decision making, which can account for the perceptual switches observed in human experiments. We further validate this model with fMRI datasets from an experiment on semantic priming in perception of ambivalent (male versus female) faces. We demonstrate that the model can correctly predict the drop in the variance of the BOLD within the Superior Parietal Area and Inferior Parietal Area while watching ambiguous visual stimuli.

Critical synchronization and 1/f noise in inhibitory/excitatory rich-club neural networks

In recent years, diverse studies have reported that different brain regions, which are internally densely connected, are also highly connected to each other. This configuration seems to play a key role in integrating and interchanging information between brain areas. Also, changes in the rich-club connectivity and the shift from inhibitory to excitatory behavior of hub neurons have been associated with several diseases. However, there is not a clear understanding about the role of the proportion of inhibitory/excitatory hub neurons, the dynamic consequences of rich-club disconnection, and hub inhibitory/excitatory shifts. Here, we study the synchronization and temporal correlations in the neural Izhikevich model, which comprises excitatory and inhibitory neurons located in a scale-free hierarchical network with rich-club connectivity. We evaluated the temporal autocorrelations and global synchronization dynamics displayed by the system in terms of rich-club connectivity and hub inhibitory/excitatory population. We evaluated the synchrony between pairs of sets of neurons by means of the global lability synchronization, based on the rate of change in the total number of synchronized signals. The results show that for a wide range of excitatory/inhibitory hub ratios the network displays 1/f dynamics with critical synchronization that is concordant with numerous health brain registers, while a network configuration with a vast majority of excitatory hubs mostly exhibits short-term autocorrelations with numerous large avalanches. Furthermore, rich-club connectivity promotes the increase of the global lability of synchrony and the temporal persistence of the system.

Theta-gamma cascades and running speed

Oscillations in the hippocampal local field potential at theta and gamma frequencies are prominent during awake behavior and have demonstrated several behavioral correlates. Both oscillations have been observed to increase in amplitude and frequency as a function of running speed. Previous investigations, however, have examined the relationship between speed and each of these oscillation bands separately. Based on energy cascade models where "…perturbations of slow frequencies cause a cascade of energy dissipation at all frequency scales" (Buzsaki G. Rhythms of the Brain, 2006), we hypothesized that cross-frequency interactions between theta and gamma should increase as a function of speed. We examined these relationships across multiple layers of the CA1 subregion, which correspond to synaptic zones receiving different afferents. Across layers, we found a reliable correlation between the power of theta and the power of gamma, indicative of an amplitude-amplitude relationship. Moreover, there was an increase in the coherence between the power of gamma and the phase of theta, demonstrating increased phase-amplitude coupling with speed. Finally, at higher velocities, phase entrainment between theta and gamma increases. These results have important implications and provide new insights regarding how theta and gamma are integrated for neuronal circuit dynamics, with coupling strength determined by the excitatory drive within the hippocampus. Specifically, rather than arguing that different frequencies can be attributed to different psychological processes, we contend that cognitive processes occur across multiple frequency bands simultaneously with organization occurring as a function of the amount of energy iteratively propagated through the brain. NEW & NOTEWORTHY Often, the theta and gamma oscillations in the hippocampus have been believed to be a consequence of two marginally overlapping phenomena. This perspective, however, runs counter to an alternative hypothesis in which a slow-frequency, high-amplitude oscillation provides energy that cascades into higher frequency, lower amplitude oscillations. We found that as running speed increases, all measures of cross-frequency theta-gamma coupling intensify, providing evidence in favor of the energy cascade hypothesis.

The multidimensional brain

Brain activity takes place in three spatial-plus time dimensions. This rather obvious claim has been recently questioned by papers that, taking into account the big data outburst and novel available computational tools, are starting to unveil a more intricate state of affairs. Indeed, various brain activities and their correlated mental functions can be assessed in terms of trajectories embedded in phase spaces of dimensions higher than the canonical ones. In this review, I show how further dimensions may not just represent a convenient methodological tool that allows a better mathematical treatment of otherwise elusive cortical activities, but may also reflect genuine functional or anatomical relationships among real nervous functions. I then describe how to extract hidden multidimensional information from real or artificial neurodata series, and make clear how our mind dilutes, rather than concentrates as currently believed, inputs coming from the environment. Finally, I argue that the principle "the higher the dimension, the greater the information" may explain the occurrence of mental activities and elucidate the mechanisms of human diseases associated with dimensionality reduction.

Wave Turbulence and Energy Cascade in the Hippocampus

Mesoscale cortical activity can be defined as the organization of activity of large neuron populations into collective action, forming time-dependent patterns such as traveling waves. Although collective action may play an important role in the cross-scale integration of brain activity and in the emergence of cognitive behavior, a comprehensive formulation of the laws governing its dynamics is still lacking. Because collective action processes are macroscopic with respect to neuronal activity, these processes cannot be described directly with methods and models developed for the microscale (individual neurons).To identify the characteristic features of mesoscopic dynamics, and to lay the foundations for a theoretical description of mesoscopic activity in the hippocampus, we conduct a comprehensive examination of observational data of hippocampal local field potential (LFP) recordings. We use the strong correlation between rat running-speed and the LFP power to parameterize the energy input into the hippocampus, and show that both the power and non-linearity of collective action (e.g., theta and gamma rhythms) increase with increased speed. Our results show that collective-action dynamics are stochastic (the precise state of a single neuron is irrelevant), weakly non-linear, and weakly dissipative. These are the principles of the theory of weak turbulence. Therefore, we propose weak turbulence a theoretical framework for the description of mesoscopic activity in the hippocampus. The weak turbulence framework provides a complete description of the cross-scale energy exchange (the energy cascade). It uncovers the mechanism governing major features of LFP spectra and bispectra, such as the physical meaning of the exponent α of power-law LFP spectra (e.g., f ^-α, where f is the frequency), the strengthening of theta-gamma coupling with energy input into the hippocampus, as well as specific phase lags associated with their interaction. Remarkably, the weak turbulence framework is consistent with the theory of self organized criticality, which provides a simple explanation for the existence of the power-law background spectrum. Together with self-organized criticality, weak turbulence could provide a unifying approach to modeling the dynamics of mesoscopic activity.

Separating normosmic and anosmic patients based on entropy evaluation of olfactory event-related potentials

OBJECTIVE

Methods based on electroencephalography (EEG) are used to evaluate brain responses to odors which is challenging due to the relatively low signal-to-noise ratio. This is especially difficult in patients with olfactory loss. In the present study, we aim to establish a method to separate functionally anosmic and normosmic individuals by means of recordings of olfactory event-related potentials (OERP) using an automated tool. Therefore, Shannon entropy was adopted to examine the complexity of the averaged electrophysiological responses.

METHODS

A total of 102 participants received 60 rose-like odorous stimuli at an inter-stimulus interval of 10 s. Olfactory-related brain activity was investigated within three time-windows of equal length; pre-, during-, and post-stimulus.

RESULTS

Based on entropy analysis, patients were correctly diagnosed for anosmia with a 75% success rate.

CONCLUSION

This novel approach can be expected to help clinicians to identify patients with anosmia or patients with early symptoms of neurodegenerative disorders.

SIGNIFICANCE

There is no automated diagnostic tool for anosmic and normosmic patients using OERP. However, detectability of OERP in patients with functional anosmia has been reported to be in the range of 50%.

Oscillatory Patterns of Phase Cone Formations near to Epileptic Spikes Derived from 256-Channel Scalp EEG Data

Our objective was to determine if there are any distinguishable phase cone clustering patterns present near to epileptic spikes. These phase cones arise from episodic phase shifts due to the coordinated activity of cortical neurons at or near to state transitions and can be extracted from the high-density scalp EEG recordings. The phase cone clustering activities in the low gamma band (30-50 Hz) and in the ripple band (80-150 Hz) were extracted from the analytic phase after taking Hilbert transform of the 256-channel high density (dEEG) data of adult patients. We used three subjects in this study. Spatiotemporal contour plots of the unwrapped analytic phase with 1.0 ms intervals were constructed using a montage layout of 256 electrode positions. Stable phase cone patterns were selected based on the criteria that the sign of the spatial gradient did not change for at least three consecutive time samples and the frame velocity was within the range of propagation velocities of cortical axons. These plots exhibited dynamical formation of phase cones which were higher in the seizure area as compared with the nearby surrounding brain areas. Spatiotemporal oscillatory patterns were also visible during ±5 sec period from the location of the spike. These results suggest that the phase cone activity might be useful for noninvasive localization of epileptic sites and also for examining the cortical neurodynamics near to epileptic spikes.

Operating in a Reverberating Regime Enables Rapid Tuning of Network States to Task Requirements

Neural circuits are able to perform computations under very diverse conditions and requirements. The required computations impose clear constraints on their fine-tuning: a rapid and maximally informative response to stimuli in general requires decorrelated baseline neural activity. Such network dynamics is known as asynchronous-irregular. In contrast, spatio-temporal integration of information requires maintenance and transfer of stimulus information over extended time periods. This can be realized at criticality, a phase transition where correlations, sensitivity and integration time diverge. Being able to flexibly switch, or even combine the above properties in a task-dependent manner would present a clear functional advantage. We propose that cortex operates in a “reverberating regime” because it is particularly favorable for ready adaptation of computational properties to context and task. This reverberating regime enables cortical networks to interpolate between the asynchronous-irregular and the critical state by small changes in effective synaptic strength or excitation-inhibition ratio. These changes directly adapt computational properties, including sensitivity, amplification, integration time and correlation length within the local network. We review recent converging evidence that cortex in vivo operates in the reverberating regime, and that various cortical areas have adapted their integration times to processing requirements. In addition, we propose that neuromodulation enables a fine-tuning of the network, so that local circuits can either decorrelate or integrate, and quench or maintain their input depending on task. We argue that this task-dependent tuning, which we call “dynamic adaptive computation,” presents a central organization principle of cortical networks and discuss first experimental evidence.

Role of Network Science in the Study of Anesthetic State Transitions

The heterogeneity of molecular mechanisms, target neural circuits, and neurophysiologic effects of general anesthetics makes it difficult to develop a reliable and drug-invariant index of general anesthesia. No single brain region or mechanism has been identified as the neural correlate of consciousness, suggesting that consciousness might emerge through complex interactions of spatially and temporally distributed brain functions. The goal of this review article is to introduce the basic concepts of networks and explain why the application of network science to general anesthesia could be a pathway to discover a fundamental mechanism of anesthetic-induced unconsciousness. This article reviews data suggesting that reduced network efficiency, constrained network repertoires, and changes in cortical dynamics create inhospitable conditions for information processing and transfer, which lead to unconsciousness. This review proposes that network science is not just a useful tool but a necessary theoretical framework and method to uncover common principles of anesthetic-induced unconsciousness.

Measuring neuronal avalanches in disordered systems with absorbing states

Power-law-shaped avalanche-size distributions are widely used to probe for critical behavior in many different systems, particularly in neural networks. The definition of avalanche is ambiguous. Usually, theoretical avalanches are defined as the activity between a stimulus and the relaxation to an inactive absorbing state. On the other hand, experimental neuronal avalanches are defined by the activity between consecutive silent states. We claim that the latter definition may be extended to some theoretical models to characterize their power-law avalanches and critical behavior. We study a system in which the separation of driving and relaxation time scales emerges from its structure. We apply both definitions of avalanche to our model. Both yield power-law-distributed avalanches that scale with system size in the critical point as expected. Nevertheless, we find restricted power-law-distributed avalanches outside of the critical region within the experimental procedure, which is not expected by the standard theoretical definition. We remark that these results are dependent on the model details.

Beta-Rhythm Oscillations and Synchronization Transition in Network Models of Izhikevich Neurons: Effect of Topology and Synaptic Type

Despite their significant functional roles, beta-band oscillations are least understood. Synchronization in neuronal networks have attracted much attention in recent years with the main focus on transition type. Whether one obtains explosive transition or a continuous transition is an important feature of the neuronal network which can depend on network structure as well as synaptic types. In this study we consider the effect of synaptic interaction (electrical and chemical) as well as structural connectivity on synchronization transition in network models of Izhikevich neurons which spike regularly with beta rhythms. We find a wide range of behavior including continuous transition, explosive transition, as well as lack of global order. The stronger electrical synapses are more conducive to synchronization and can even lead to explosive synchronization. The key network element which determines the order of transition is found to be the clustering coefficient and not the small world effect, or the existence of hubs in a network. These results are in contrast to previous results which use phase oscillator models such as the Kuramoto model. Furthermore, we show that the patterns of synchronization changes when one goes to the gamma band. We attribute such a change to the change in the refractory period of Izhikevich neurons which changes significantly with frequency.

Scale invariance in natural and artificial collective systems: a review

Self-organized collective coordinated behaviour is an impressive phenomenon, observed in a variety of natural and artificial systems, in which coherent global structures or dynamics emerge from local interactions between individual parts. If the degree of collective integration of a system does not depend on size, its level of robustness and adaptivity is typically increased and we refer to it as scale-invariant. In this review, we first identify three main types of self-organized scale-invariant systems: scale-invariant spatial structures, scale-invariant topologies and scale-invariant dynamics. We then provide examples of scale invariance from different domains in science, describe their origins and main features and discuss potential challenges and approaches for designing and engineering artificial systems with scale-invariant properties.

1/f electrophysiological spectra in resting and drug-induced states can be explained by the dynamics of multiple oscillatory relaxation processes

Neurophysiological recordings are dominated by arhythmical activity whose spectra can be characterised by power-law functions, and on this basis are often referred to as reflecting scale-free brain dynamics (1/f^β). Relatively little is known regarding the neural generators and temporal dynamics of this arhythmical behaviour compared to rhythmical behaviour. Here we used Irregularly Resampled AutoSpectral Analysis (IRASA) to quantify β, in both the high (5-100 Hz, β_hf) and low frequency bands (0.1-2.5 Hz, β_lf) in MEG/EEG/ECoG recordings and to separate arhythmical from rhythmical modes of activity, such as, alpha rhythms. In MEG/EEG/ECoG data, we demonstrate that oscillatory alpha power dynamically correlates over time with β_hf and similarly, participants with higher rhythmical alpha power have higher β_hf. In a series of pharmaco-MEG investigations using the GABA reuptake inhibitor tiagabine, the glutamatergic AMPA receptor antagonist perampanel, the NMDA receptor antagonist ketamine and the mixed partial serotonergic agonist LSD, a variety of effects on both β_hf and β_lf were observed. Additionally, strong modulations of β_hf were seen in monkey ECoG data during general anaesthesia using propofol and ketamine. We develop and test a unifying model which can explain, the 1/f nature of electrophysiological spectra, their dynamic interaction with oscillatory rhythms as well as the sensitivity of 1/f activity to drug interventions by considering electrophysiological spectra as being generated by a collection of stochastically perturbed damped oscillators having a distribution of relaxation rates.

Criticality in the Brain: Evidence and Implications for Neuromorphic Computing

We have discovered an unexpected correlation between the operational temperature of the brain and cognitive abilities across a wide variety of animal species. This correlation is extracted from available data in the literature of the temperature range Δ T at which an animal's brain can operate and its encephalization quotient EQ, which can be used as a proxy for cognitive ability. In particular, we found a power-law dependence between Δ T and EQ. These data support the theory that the brain behaves as a critical system where temperature is one of the critical parameters, tuning the performance of the neural network.

Griffiths phases in infinite-dimensional, non-hierarchical modular networks

Griffiths phases (GPs), generated by the heterogeneities on modular networks, have recently been suggested to provide a mechanism, rid of fine parameter tuning, to explain the critical behavior of complex systems. One conjectured requirement for systems with modular structures was that the network of modules must be hierarchically organized and possess finite dimension. We investigate the dynamical behavior of an activity spreading model, evolving on heterogeneous random networks with highly modular structure and organized non-hierarchically. We observe that loosely coupled modules act as effective rare-regions, slowing down the extinction of activation. As a consequence, we find extended control parameter regions with continuously changing dynamical exponents for single network realizations, preserved after finite size analyses, as in a real GP. The avalanche size distributions of spreading events exhibit robust power-law tails. Our findings relax the requirement of hierarchical organization of the modular structure, which can help to rationalize the criticality of modular systems in the framework of GPs.

A Basic Phase Diagram of Neuronal Dynamics

The extreme complexity of the brain has attracted the attention of neuroscientists and other researchers for a long time. More recently, the neuromorphic hardware has matured to provide a new powerful tool to study neuronal dynamics. Here, we study neuronal dynamics using different settings on a neuromorphic chip built with flexible parameters of neuron models. Our unique setting in the network of leaky integrate-and-fire (LIF) neurons is to introduce a weak noise environment. We observed three different types of collective neuronal activities, or phases, separated by sharp boundaries, or phase transitions. From this, we construct a rudimentary phase diagram of neuronal dynamics and demonstrate that a noise-induced chaotic phase (N-phase), which is dominated by neuronal avalanche activity (intermittent aperiodic neuron firing), emerges in the presence of noise and its width grows with the noise intensity. The dynamics can be manipulated in this N-phase. Our results and comparison with clinical data is consistent with the literature and our previous work showing that healthy brain must reside in the N-phase. We argue that the brain phase diagram with further refinement may be used for the diagnosis and treatment of mental disease and also suggest that the dynamics may be manipulated to serve as a means of new information processing (e.g., for optimization). Neuromorphic chips, similar to the one we used but with a variety of neuron models, may be used to further enhance the understanding of human brain function and accelerate the development of neuroscience research.

Hysteresis, neural avalanches, and critical behavior near a first-order transition of a spiking neural network

Many experimental results, both in vivo and in vitro, support the idea that the brain cortex operates near a critical point and at the same time works as a reservoir of precise spatiotemporal patterns. However, the mechanism at the basis of these observations is still not clear. In this paper we introduce a model which combines both these features, showing that scale-free avalanches are the signature of a system posed near the spinodal line of a first-order transition, with many spatiotemporal patterns stored as dynamical metastable attractors. Specifically, we studied a network of leaky integrate-and-fire neurons whose connections are the result of the learning of multiple spatiotemporal dynamical patterns, each with a randomly chosen ordering of the neurons. We found that the network shows a first-order transition between a low-spiking-rate disordered state (down), and a high-rate state characterized by the emergence of collective activity and the replay of one of the stored patterns (up). The transition is characterized by hysteresis, or alternation of up and down states, depending on the lifetime of the metastable states. In both cases, critical features and neural avalanches are observed. Notably, critical phenomena occur at the edge of a discontinuous phase transition, as recently observed in a network of glow lamps.

Efficient Coding and Energy Efficiency Are Promoted by Balanced Excitatory and Inhibitory Synaptic Currents in Neuronal Network

Selective pressure may drive neural systems to process as much information as possible with the lowest energy cost. Recent experiment evidence revealed that the ratio between synaptic excitation and inhibition (E/I) in local cortex is generally maintained at a certain value which may influence the efficiency of energy consumption and information transmission of neural networks. To understand this issue deeply, we constructed a typical recurrent Hodgkin-Huxley network model and studied the general principles that governs the relationship among the E/I synaptic current ratio, the energy cost and total amount of information transmission. We observed in such a network that there exists an optimal E/I synaptic current ratio in the network by which the information transmission achieves the maximum with relatively low energy cost. The coding energy efficiency which is defined as the mutual information divided by the energy cost, achieved the maximum with the balanced synaptic current. Although background noise degrades information transmission and imposes an additional energy cost, we find an optimal noise intensity that yields the largest information transmission and energy efficiency at this optimal E/I synaptic transmission ratio. The maximization of energy efficiency also requires a certain part of energy cost associated with spontaneous spiking and synaptic activities. We further proved this finding with analytical solution based on the response function of bistable neurons, and demonstrated that optimal net synaptic currents are capable of maximizing both the mutual information and energy efficiency. These results revealed that the development of E/I synaptic current balance could lead a cortical network to operate at a highly efficient information transmission rate at a relatively low energy cost. The generality of neuronal models and the recurrent network configuration used here suggest that the existence of an optimal E/I cell ratio for highly efficient energy costs and information maximization is a potential principle for cortical circuit networks.

Adaptive scales of integration and response latencies in a critically-balanced model of the primary visual cortex

The primary visual cortex (V1) integrates information over scales in visual space, which have been shown to vary, in an input-dependent manner, as a function of contrast and other visual parameters. Which algorithms the brain uses to achieve this feat are largely unknown and an open problem in visual neuroscience. We demonstrate that a simple dynamical mechanism can account for this contrast-dependent scale of integration in visuotopic space as well as connect this property to two other stimulus-dependent features of V1: extents of lateral integration on the cortical surface and response latencies.

Scale-freeness or partial synchronization in neural mass phase oscillator networks: Pick one of two?

Modeling and interpreting (partial) synchronous neural activity can be a challenge. We illustrate this by deriving the phase dynamics of two seminal neural mass models: the Wilson-Cowan firing rate model and the voltage-based Freeman model. We established that the phase dynamics of these models differed qualitatively due to an attractive coupling in the first and a repulsive coupling in the latter. Using empirical structural connectivity matrices, we determined that the two dynamics cover the functional connectivity observed in resting state activity. We further searched for two pivotal dynamical features that have been reported in many experimental studies: (1) a partial phase synchrony with a possibility of a transition towards either a desynchronized or a (fully) synchronized state; (2) long-term autocorrelations indicative of a scale-free temporal dynamics of phase synchronization. Only the Freeman phase model exhibited scale-free behavior. Its repulsive coupling, however, let the individual phases disperse and did not allow for a transition into a synchronized state. The Wilson-Cowan phase model, by contrast, could switch into a (partially) synchronized state, but it did not generate long-term correlations although being located close to the onset of synchronization, i.e. in its critical regime. That is, the phase-reduced models can display one of the two dynamical features, but not both.

Beyond the evoked/intrinsic neural process dichotomy

Contemporary functional neuroimaging research has increasingly focused on characterization of intrinsic or "spontaneous" brain activity. Analysis of intrinsic activity is often contrasted with analysis of task-evoked activity that has traditionally been the focus of cognitive neuroscience. But does this evoked/intrinsic dichotomy adequately characterize human brain function? Based on empirical data demonstrating a close functional interdependence between intrinsic and task-evoked activity, we argue that the dichotomy between intrinsic and task-evoked activity as unobserved contributions to brain activity is artificial. We present an alternative picture of brain function in which the brain's spatiotemporal dynamics do not consist of separable intrinsic and task-evoked components, but reflect the enaction of a system of mutual constraints to move the brain into and out of task-appropriate functional configurations. According to this alternative picture, cognitive neuroscientists are tasked with describing both the temporal trajectory of brain activity patterns across time, and the modulation of this trajectory by task states, without separating this process into intrinsic and task-evoked components. We argue that this alternative picture of brain function is best captured in a novel explanatory framework called enabling constraint. Overall, these insights call for a reconceptualization of functional brain activity, and should drive future methodological and empirical efforts.

Optimization by Self-Organized Criticality

Self-organized criticality (SOC) is a phenomenon observed in certain complex systems of multiple interacting components, e.g., neural networks, forest fires, and power grids, that produce power-law distributed avalanche sizes. Here, we report the surprising result that the avalanches from an SOC process can be used to solve non-convex optimization problems. To generate avalanches, we use the Abelian sandpile model on a graph that mirrors the graph of the optimization problem. For optimization, we map the avalanche areas onto search patterns for optimization, while the SOC process receives no feedback from the optimization itself. The resulting method can be applied without parameter tuning to a wide range of optimization problems, as demonstrated on three problems: finding the ground-state of an Ising spin glass, graph coloring, and image segmentation. We find that SOC search is more efficient compared to other random search methods, including simulated annealing, and unlike annealing, it is parameter free, thereby eliminating the time-consuming requirement to tune an annealing temperature schedule.

Coexistence of Stochastic Oscillations and Self-Organized Criticality in a Neuronal Network: Sandpile Model Application

Self-organized criticality (SOC) and stochastic oscillations (SOs) are two theoretically contradictory phenomena that are suggested to coexist in the brain. Recently it has been shown that an accumulation-release process like sandpile dynamics can generate SOC and SOs simultaneously. We considered the effect of the network structure on this coexistence and showed that the sandpile dynamics on a small-world network can produce two power law regimes along with two groups of SOs-two peaks in the power spectrum of the generated signal simultaneously. We also showed that external stimuli in the sandpile dynamics do not affect the coexistence of SOC and SOs but increase the frequency of SOs, which is consistent with our knowledge of the brain.

Inferring multi-scale neural mechanisms with brain network modelling

The neurophysiological processes underlying non-invasive brain activity measurements are incompletely understood. Here, we developed a connectome-based brain network model that integrates individual structural and functional data with neural population dynamics to support multi-scale neurophysiological inference. Simulated populations were linked by structural connectivity and, as a novelty, driven by electroencephalography (EEG) source activity. Simulations not only predicted subjects' individual resting-state functional magnetic resonance imaging (fMRI) time series and spatial network topologies over 20 minutes of activity, but more importantly, they also revealed precise neurophysiological mechanisms that underlie and link six empirical observations from different scales and modalities: (1) resting-state fMRI oscillations, (2) functional connectivity networks, (3) excitation-inhibition balance, (4, 5) inverse relationships between α-rhythms, spike-firing and fMRI on short and long time scales, and (6) fMRI power-law scaling. These findings underscore the potential of this new modelling framework for general inference and integration of neurophysiological knowledge to complement empirical studies.

Neuroscientists can use various techniques to measure activity within the brain without opening up the skull. One of the most common is electroencephalography, or EEG for short. A net of electrodes is attached to the scalp and reveals the patterns of electrical activity occurring in brain tissue. But while EEG is good at revealing electrical activity across the surface of the scalp, it is less effective at linking the observed activity to specific locations in the brain.

Another widely used technique is functional magnetic resonance imaging, or fMRI. A patient, or healthy volunteer, lies inside a scanner containing a large magnet. The scanner tracks changes in the level of oxygen at different regions of the brain to provide a measure of how the activity of these regions changes over time. In contrast to EEG, fMRI is good at pinpointing the location of brain activity, but it is an indirect measure of brain activity as it depends on blood flow and several other factors. In terms of understanding how the brain works, EEG and fMRI thus provide different pieces of the puzzle. But there is no easy way to fit these pieces together.

Other areas of science have used computer models to merge different sources of data to obtain new insights into complex processes. Schirner et al. now adopt this approach to reveal the workings of the brain that underly signals like EEG and fMRI.

After recording structural MRI data from healthy volunteers, Schirner et al. built a computer model of each person’s brain. They then ran simulations with each individual model stimulating it with the person’s EEG to predict the fMRI activity of the same individual. Comparing these predictions with real fMRI data collected at the same time as the EEG confirmed that the predictions were accurate. Importantly, the brain models also displayed many features of neural activity that previously could only be measured by implanting electrodes into the brain.

This new approach provides a way of combining experimental data with theories about how the nervous system works. The resulting models can help generate and test ideas about the mechanisms underlying brain activity. Building models of different brains based on data from individual people could also help reveal the biological basis of differences between individuals. This could in turn provide insights into why some individuals are more vulnerable to certain brain diseases and open up new ways to treat these diseases.

Comparative approximations of criticality in a neural and quantum regime

Under a variety of conditions, stochastic and non-linear systems with many degrees of freedom tend to evolve towards complexity and criticality. Over the last decades, a steady proliferation of models re: far-from-equilibrium thermodynamics of metastable, many-valued systems arose, serving as attributes of a 'critical' attractor landscape. Building off recent data citing trademark aspects of criticality in the brain-including: power-laws, scale-free (1/f) behavior (scale invariance, or scale independence), critical slowing, and avalanches-it has been conjectured that operating at criticality entails functional advantages such as: optimized neural computation and information processing; boosted memory; large dynamical ranges; long-range communication; and an increased ability to react to highly diverse stimuli. In short, critical dynamics provide a necessary condition for neurobiologically significant elements of brain dynamics. Theoretical predictions have been verified in specific models such as Boolean networks, liquid state machines, and neural networks. These findings inspired the neural criticality hypothesis, proposing that the brain operates in a critical state because the associated optimal computational capabilities provide an evolutionarily advantage. This paper develops in three parts: after developing the critical landscape, we will then shift gears to rediscover another inroad to criticality via stochastic quantum field theory and dissipative dynamics. The existence of these two approaches deserves some consideration, given both neural and quantum criticality hypotheses propose specific mechanisms that leverage the same phenomena. This suggests that understanding the quantum approach could help to shed light on brain-based modeling. In the third part, we will turn to Whitehead's actual entities and modes of perception in order to demonstrate a concomitant logic underwriting both models. In the discussion, I briefly motivate a reading of criticality and its properties as responsive to the characterization of tenets from Eastern wisdom traditions.

Relationships between short and fast brain timescales

Brain electric activity exhibits two important features: oscillations with different timescales, characterized by diverse functional and psychological outcomes, and a temporal power law distribution. In order to further investigate the relationships between low- and high- frequency spikes in the brain, we used a variant of the Borsuk-Ulam theorem which states that, when we assess the nervous activity as embedded in a sphere equipped with a fractal dimension, we achieve two antipodal points with similar features (the slow and fast, scale-free oscillations). We demonstrate that slow and fast nervous oscillations mirror each other over time via a sinusoid relationship and provide, through the Bloch theorem from solid-state physics, the possible equation which links the two timescale activities. We show that, based on topological findings, nervous activities occurring in micro-levels are projected to single activities at meso- and macro-levels. This means that brain functions assessed at the higher scale of the whole brain necessarily display a counterpart in the lower ones, and vice versa. Our topological approach makes it possible to assess brain functions both based on entropy, and in the general terms of particle trajectories taking place on donut-like manifolds. Condensed brain activities might give rise to ideas and concepts by combination of different functional and anatomical levels. Furthermore, cognitive phenomena, as well as social activity can be described by the laws of quantum mechanics; memories and decisions exhibit holographic organization. In physics, the term duality refers to a case where two seemingly different systems turn out to be equivalent. This topological duality holds for all the types of spatio-temporal brain activities, independent of their inter- and intra-level relationships, strength, magnitude and boundaries, allowing us to connect the physiological manifestations of consciousness to the electric activities of the brain.

Neuronal avalanches: Where temporal complexity and criticality meet

The model of the current paper is an extension of a previous publication, wherein we have used the leaky integrate-and-fire model on a regular lattice with periodic boundary conditions, and introduced the temporal complexity as a genuine signature of criticality. In that work, the power-law distribution of neural avalanches was a manifestation of supercriticality rather than criticality. Here, however, we show that the continuous solution of the model and replacing the stochastic noise with a Gaussian zero-mean noise leads to the coincidence of power-law display of temporal complexity, and spatiotemporal patterns of neural avalanches at the critical point. We conclude that the source of inconsistency may be a numerical artifact originated by the discrete description of the model which may imply a slow numerical convergence of the avalanche distribution compared to temporal complexity.

Revisiting the Quantum Brain Hypothesis: Toward Quantum (Neuro)biology?

The nervous system is a non-linear dynamical complex system with many feedback loops. A conventional wisdom is that in the brain the quantum fluctuations are self-averaging and thus functionally negligible. However, this intuition might be misleading in the case of non-linear complex systems. Because of an extreme sensitivity to initial conditions, in complex systems the microscopic fluctuations may be amplified and thereby affect the system's behavior. In this way quantum dynamics might influence neuronal computations. Accumulating evidence in non-neuronal systems indicates that biological evolution is able to exploit quantum stochasticity. The recent rise of quantum biology as an emerging field at the border between quantum physics and the life sciences suggests that quantum events could play a non-trivial role also in neuronal cells. Direct experimental evidence for this is still missing but future research should address the possibility that quantum events contribute to an extremely high complexity, variability and computational power of neuronal dynamics.

Temporal profiles of avalanches on networks

An avalanche or cascade occurs when one event causes one or more subsequent events, which in turn may cause further events in a chain reaction. Avalanching dynamics are studied in many disciplines, with a recent focus on average avalanche shapes, i.e., the temporal profiles of avalanches of fixed duration. At the critical point of the dynamics, the rescaled average avalanche shapes for different durations collapse onto a single universal curve. We apply Markov branching process theory to derive an equation governing the average avalanche shape for cascade dynamics on networks. Analysis of the equation at criticality demonstrates that nonsymmetric average avalanche shapes (as observed in some experiments) occur for certain combinations of dynamics and network topology. We give examples using numerical simulations of models for information spreading, neural dynamics, and behavior adoption and we propose simple experimental tests to quantify whether cascading systems are in the critical state.

Signatures of criticality arise from random subsampling in simple population models

The rise of large-scale recordings of neuronal activity has fueled the hope to gain new insights into the collective activity of neural ensembles. How can one link the statistics of neural population activity to underlying principles and theories? One attempt to interpret such data builds upon analogies to the behaviour of collective systems in statistical physics. Divergence of the specific heat-a measure of population statistics derived from thermodynamics-has been used to suggest that neural populations are optimized to operate at a "critical point". However, these findings have been challenged by theoretical studies which have shown that common inputs can lead to diverging specific heat. Here, we connect "signatures of criticality", and in particular the divergence of specific heat, back to statistics of neural population activity commonly studied in neural coding: firing rates and pairwise correlations. We show that the specific heat diverges whenever the average correlation strength does not depend on population size. This is necessarily true when data with correlations is randomly subsampled during the analysis process, irrespective of the detailed structure or origin of correlations. We also show how the characteristic shape of specific heat capacity curves depends on firing rates and correlations, using both analytically tractable models and numerical simulations of a canonical feed-forward population model. To analyze these simulations, we develop efficient methods for characterizing large-scale neural population activity with maximum entropy models. We find that, consistent with experimental findings, increases in firing rates and correlation directly lead to more pronounced signatures. Thus, previous reports of thermodynamical criticality in neural populations based on the analysis of specific heat can be explained by average firing rates and correlations, and are not indicative of an optimized coding strategy. We conclude that a reliable interpretation of statistical tests for theories of neural coding is possible only in reference to relevant ground-truth models.

Dynamical Signatures of Structural Connectivity Damage to a Model of the Brain Posed at Criticality

Synchronization of brain activity fluctuations is believed to represent communication between spatially distant neural processes. These interareal functional interactions develop in the background of a complex network of axonal connections linking cortical and subcortical neurons, termed the human "structural connectome." Theoretical considerations and experimental evidence support the view that the human brain can be modeled as a system operating at a critical point between ordered (subcritical) and disordered (supercritical) phases. Here, we explore the hypothesis that pathologies resulting from brain injury of different etiologies are related to this model of a critical brain. For this purpose, we investigate how damage to the integrity of the structural connectome impacts on the signatures of critical dynamics. Adopting a hybrid modeling approach combining an empirical weighted network of human structural connections with a conceptual model of critical dynamics, we show that lesions located at highly transited connections progressively displace the model toward the subcritical regime. The topological properties of the nodes and links are of less importance when considered independently of their weight in the network. We observe that damage to midline hubs such as the middle and posterior cingulate cortex is most crucial for the disruption of criticality in the model. However, a similar effect can be achieved by targeting less transited nodes and links whose connection weights add up to an equivalent amount. This implies that brain pathology does not necessarily arise due to insult targeted at well-connected areas and that intersubject variability could obscure lesions located at nonhub regions. Finally, we discuss the predictions of our model in the context of clinical studies of traumatic brain injury and neurodegenerative disorders.

Deviations from Critical Dynamics in Interictal Epileptiform Activity

The framework of criticality provides a unifying perspective on neuronal dynamics from in vitro cortical cultures to functioning human brains. Recent findings suggest that a healthy cortex displays critical dynamics, giving rise to scale-free spatiotemporal cascades of activity, termed neuronal avalanches. Pharmacological manipulations of the excitation-inhibition balance (EIB) in cortical cultures were previously shown to result in deviations from criticality and from the power law scaling of avalanche size distribution. To examine the sensitivity of neuronal avalanche metrics to altered EIB in humans, we focused on epilepsy, a neurological disorder characterized by hyperexcitable networks. Using magnetoencephalography, we quantitatively assessed deviations from criticality in the brain dynamics of patients with epilepsy during interictal (between-seizures) activity. Compared with healthy control subjects, epilepsy patients tended to exhibit a higher neural gain and larger avalanches, particularly during interictal epileptiform activity. Moreover, deviations from scale-free behavior were exclusively connected to brief intervals at epileptiform discharges, strengthening the association between deviations from criticality and the instantaneous changes in EIB. The avalanches collected during interictal epileptiform activity had not only a stereotypical size range but also involved particular spatial patterns of activations, as expected for periods of epileptic network dominance. Overall, the neuronal avalanche metrics provide a quantitative novel description of interictal brain activity of patients with epilepsy.

SIGNIFICANCE STATEMENT

Healthy brain dynamics requires a delicate balance between excitatory and inhibitory processes. Several brain disorders, such as epilepsy, are associated with altered excitation-inhibition balance, but assessing this balance using noninvasive tools is still challenging. In this study, we apply the framework of critical brain dynamics to data from epilepsy patients, which were recorded between seizures. We show that metrics of criticality provide a sensitive tool for noninvasive assessment of changes in the balance. Specifically, brain activity of epilepsy patients deviates from healthy critical brain dynamics, particularly during abnormal epileptiform activity. The study offers a novel quantitative perspective on epilepsy and its relation to healthy brain dynamics.

Criticality predicts maximum irregularity in recurrent networks of excitatory nodes

A rigorous understanding of brain dynamics and function requires a conceptual bridge between multiple levels of organization, including neural spiking and network-level population activity. Mounting evidence suggests that neural networks of cerebral cortex operate at a critical regime, which is defined as a transition point between two phases of short lasting and chaotic activity. However, despite the fact that criticality brings about certain functional advantages for information processing, its supporting evidence is still far from conclusive, as it has been mostly based on power law scaling of size and durations of cascades of activity. Moreover, to what degree such hypothesis could explain some fundamental features of neural activity is still largely unknown. One of the most prevalent features of cortical activity in vivo is known to be spike irregularity of spike trains, which is measured in terms of the coefficient of variation (CV) larger than one. Here, using a minimal computational model of excitatory nodes, we show that irregular spiking (CV > 1) naturally emerges in a recurrent network operating at criticality. More importantly, we show that even at the presence of other sources of spike irregularity, being at criticality maximizes the mean coefficient of variation of neurons, thereby maximizing their spike irregularity. Furthermore, we also show that such a maximized irregularity results in maximum correlation between neuronal firing rates and their corresponding spike irregularity (measured in terms of CV). On the one hand, using a model in the universality class of directed percolation, we propose new hallmarks of criticality at single-unit level, which could be applicable to any network of excitable nodes. On the other hand, given the controversy of the neural criticality hypothesis, we discuss the limitation of this approach to neural systems and to what degree they support the criticality hypothesis in real neural networks. Finally, we discuss the limitations of applying our results to real networks and to what degree they support the criticality hypothesis.

Refractory period in network models of excitable nodes: self-sustaining stable dynamics, extended scaling region and oscillatory behavior

Networks of excitable nodes have recently attracted much attention particularly in regards to neuronal dynamics, where criticality has been argued to be a fundamental property. Refractory behavior, which limits the excitability of neurons is thought to be an important dynamical property. We therefore consider a simple model of excitable nodes which is known to exhibit a transition to instability at a critical point (λ = 1), and introduce refractory period into its dynamics. We use mean-field analytical calculations as well as numerical simulations to calculate the activity dependent branching ratio that is useful to characterize the behavior of critical systems. We also define avalanches and calculate probability distribution of their size and duration. We find that in the presence of refractory period the dynamics stabilizes while various parameter regimes become accessible. A sub-critical regime with λ < 1.0, a standard critical behavior with exponents close to critical branching process for λ = 1, a regime with 1 < λ < 2 that exhibits an interesting scaling behavior, and an oscillating regime with λ > 2.0. We have therefore shown that refractory behavior leads to a wide range of scaling as well as periodic behavior which are relevant to real neuronal dynamics.

Optimizing Neural Information Capacity through Discretization

Discretization in neural circuits occurs on many levels, from the generation of action potentials and dendritic integration, to neuropeptide signaling and processing of signals from multiple neurons, to behavioral decisions. It is clear that discretization, when implemented properly, can convey many benefits. However, the optimal solutions depend on both the level of noise and how it impacts a particular computation. This Perspective discusses how current physiological data could potentially be integrated into one theoretical framework based on maximizing information. Key experiments for testing that framework are discussed.

Criticality in the brain: A synthesis of neurobiology, models and cognition

Cognitive function requires the coordination of neural activity across many scales, from neurons and circuits to large-scale networks. As such, it is unlikely that an explanatory framework focused upon any single scale will yield a comprehensive theory of brain activity and cognitive function. Modelling and analysis methods for neuroscience should aim to accommodate multiscale phenomena. Emerging research now suggests that multi-scale processes in the brain arise from so-called critical phenomena that occur very broadly in the natural world. Criticality arises in complex systems perched between order and disorder, and is marked by fluctuations that do not have any privileged spatial or temporal scale. We review the core nature of criticality, the evidence supporting its role in neural systems and its explanatory potential in brain health and disease.

A Unified Theory of Neuro-MRI Data Shows Scale-Free Nature of Connectivity Modes

A primary goal of many neuroimaging studies that use magnetic resonance imaging (MRI) is to deduce the structure-function relationships in the human brain using data from the three major neuro-MRI modalities: high-resolution anatomical, diffusion tensor imaging, and functional MRI. To date, the general procedure for analyzing these data is to combine the results derived independently from each of these modalities. In this article, we develop a new theoretical and computational approach for combining these different MRI modalities into a powerful and versatile framework that combines our recently developed methods for morphological shape analysis and segmentation, simultaneous local diffusion estimation and global tractography, and nonlinear and nongaussian spatial-temporal activation pattern classification and ranking, as well as our fast and accurate approach for nonlinear registration between modalities. This joint analysis method is capable of extracting new levels of information that is not achievable from any of those single modalities alone. A theoretical probabilistic framework based on a reformulation of prior information and available interdependencies between modalities through a joint coupling matrix and an efficient computational implementation allows construction of quantitative functional, structural, and effective brain connectivity modes and parcellation. This new method provides an overall increase of resolution, accuracy, level of detail, and information content and has the potential to be instrumental in the clinical adaptation of neuro-MRI modalities, which, when jointly analyzed, provide a more comprehensive view of a subject's structure-function relations, while the current standard, wherein single-modality methods are analyzed separately, leaves a critical gap in an integrated view of a subject's neuorphysiological state. As one example of this increased sensitivity, we demonstrate that the jointly estimated structural and functional dependencies of mode power follow the same power law decay with the same exponent.

Modular co-organization of functional connectivity and scale-free dynamics in the human brain

Scale-free neuronal dynamics and interareal correlations are emergent characteristics of spontaneous brain activity. How such dynamics and the anatomical patterns of neuronal connectivity are mutually related in brain networks has, however, remained unclear. We addressed this relationship by quantifying the network colocalization of scale-free neuronal activity-both neuronal avalanches and long-range temporal correlations (LRTCs)-and functional connectivity (FC) by means of intracranial and noninvasive human resting-state electrophysiological recordings. We found frequency-specific colocalization of scale-free dynamics and FC so that the interareal couplings of LRTCs and the propagation of neuronal avalanches were most pronounced in the predominant pathways of FC. Several control analyses and the frequency specificity of network colocalization showed that the results were not trivial by-products of either brain dynamics or our analysis approach. Crucially, scale-free neuronal dynamics and connectivity also had colocalized modular structures at multiple levels of network organization, suggesting that modules of FC would be endowed with partially independent dynamic states. These findings thus suggest that FC and scale-free dynamics-hence, putatively, neuronal criticality as well-coemerge in a hierarchically modular structure in which the modules are characterized by dense connectivity, avalanche propagation, and shared dynamic states.

Criticality meets learning: Criticality signatures in a self-organizing recurrent neural network

Many experiments have suggested that the brain operates close to a critical state, based on signatures of criticality such as power-law distributed neuronal avalanches. In neural network models, criticality is a dynamical state that maximizes information processing capacities, e.g. sensitivity to input, dynamical range and storage capacity, which makes it a favorable candidate state for brain function. Although models that self-organize towards a critical state have been proposed, the relation between criticality signatures and learning is still unclear. Here, we investigate signatures of criticality in a self-organizing recurrent neural network (SORN). Investigating criticality in the SORN is of particular interest because it has not been developed to show criticality. Instead, the SORN has been shown to exhibit spatio-temporal pattern learning through a combination of neural plasticity mechanisms and it reproduces a number of biological findings on neural variability and the statistics and fluctuations of synaptic efficacies. We show that, after a transient, the SORN spontaneously self-organizes into a dynamical state that shows criticality signatures comparable to those found in experiments. The plasticity mechanisms are necessary to attain that dynamical state, but not to maintain it. Furthermore, onset of external input transiently changes the slope of the avalanche distributions – matching recent experimental findings. Interestingly, the membrane noise level necessary for the occurrence of the criticality signatures reduces the model’s performance in simple learning tasks. Overall, our work shows that the biologically inspired plasticity and homeostasis mechanisms responsible for the SORN’s spatio-temporal learning abilities can give rise to criticality signatures in its activity when driven by random input, but these break down under the structured input of short repeating sequences.

A distribution model of functional connectome based on criticality and energy constraints

The analysis of the network structure of the functional connectivity data constructed from fMRI images provides basic information about functions and features of the brain activity. We focus on the two features which are considered as relevant to the brain activity, the criticality and the constraint regarding energy consumptions. Within a wide variety of complex systems, the critical state occurs associated with a phase transition between distinct phases, random one and order one. Although the hypothesis that human brain activity is also in a state of criticality is supported by some experimental results, it still remains controversial. One issue is that experimental distributions exhibit deviations from the power law predicted by the criticality. Based on the assumption that constraints on brain from the biological costs cause these deviations, we derive a distribution model. The evaluation using the information criteria indicates an advantage of this model in fitting to experimental data compared to other representative distribution models, the truncated power law and the power law. Our findings also suggest that the mechanism underlying this model is closely related to the cost effective behavior in human brain with maximizing the network efficiency for the given network cost.

Neocortical activity is stimulus- and scale-invariant

Mounting evidence supports the hypothesis that the cortex operates near a critical state, defined as the transition point between order (large-scale activity) and disorder (small-scale activity). This criticality is manifested by power law distribution of the size and duration of spontaneous cascades of activity, which are referred as neuronal avalanches. The existence of such neuronal avalanches has been confirmed by several studies both in vitro and in vivo, among different species and across multiple spatial scales. However, despite the prevalence of scale free activity, still very little is known concerning whether and how the scale-free nature of cortical activity is altered during external stimulation. To address this question, we performed in vivo two-photon population calcium imaging of layer 2/3 neurons in primary visual cortex of behaving mice during visual stimulation and conducted statistical analyses on the inferred spike trains. Our investigation for each mouse and condition revealed power law distributed neuronal avalanches, and irregular spiking individual neurons. Importantly, both the avalanche and the spike train properties remained largely unchanged for different stimuli, while the cross-correlation structure varied with stimuli. Our results establish that microcircuits in the visual cortex operate near the critical regime, while rearranging functional connectivity in response to varying sensory inputs.

Injured Brains and Adaptive Networks: The Benefits and Costs of Hyperconnectivity

A common finding in human functional brain-imaging studies is that damage to neural systems paradoxically results in enhanced functional connectivity between network regions, a phenomenon commonly referred to as 'hyperconnectivity'. Here, we describe the various ways that hyperconnectivity operates to benefit a neural network following injury while simultaneously negotiating the trade-off between metabolic cost and communication efficiency. Hyperconnectivity may be optimally expressed by increasing connections through the most central and metabolically efficient regions (i.e., hubs). While adaptive in the short term, we propose that chronic hyperconnectivity may leave network hubs vulnerable to secondary pathological processes over the life span due to chronically elevated metabolic stress. We conclude by offering novel, testable hypotheses for advancing our understanding of the role of hyperconnectivity in systems-level brain plasticity in neurological disorders.

Spontaneous cortical activity is transiently poised close to criticality

Brain activity displays a large repertoire of dynamics across the sleep-wake cycle and even during anesthesia. It was suggested that criticality could serve as a unifying principle underlying the diversity of dynamics. This view has been supported by the observation of spontaneous bursts of cortical activity with scale-invariant sizes and durations, known as neuronal avalanches, in recordings of mesoscopic cortical signals. However, the existence of neuronal avalanches in spiking activity has been equivocal with studies reporting both its presence and absence. Here, we show that signs of criticality in spiking activity can change between synchronized and desynchronized cortical states. We analyzed the spontaneous activity in the primary visual cortex of the anesthetized cat and the awake monkey, and found that neuronal avalanches and thermodynamic indicators of criticality strongly depend on collective synchrony among neurons, LFP fluctuations, and behavioral state. We found that synchronized states are associated to criticality, large dynamical repertoire and prolonged epochs of eye closure, while desynchronized states are associated to sub-criticality, reduced dynamical repertoire, and eyes open conditions. Our results show that criticality in cortical dynamics is not stationary, but fluctuates during anesthesia and between different vigilance states.

Scale Invariance in Lateral Head Scans During Spatial Exploration

Universality connects various natural phenomena through physical principles governing their dynamics, and has provided broadly accepted answers to many complex questions, including information processing in neuronal systems. However, its significance in behavioral systems is still elusive. Lateral head scanning (LHS) behavior in rodents might contribute to spatial navigation by actively managing (optimizing) the available sensory information. Our findings of scale invariant distributions in LHS lifetimes, interevent intervals and event magnitudes, provide evidence for the first time that the optimization takes place at a critical point in LHS dynamics. We propose that the LHS behavior is responsible for preprocessing of the spatial information content, critical for subsequent foolproof encoding by the respective downstream neural networks.

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

Critical states are sometimes identified experimentally through power-law statistics or universal scaling functions. We show here that such features naturally emerge from networks in self-sustained irregular regimes away from criticality. In these regimes, statistical physics theory of large interacting systems predict a regime where the nodes have independent and identically distributed dynamics. We thus investigated the statistics of a system in which units are replaced by independent stochastic surrogates and found the same power-law statistics, indicating that these are not sufficient to establish criticality. We rather suggest that these are universal features of large-scale networks when considered macroscopically. These results put caution on the interpretation of scaling laws found in nature.

Repertoire of mesoscopic cortical activity is not reduced during anesthesia

Consciousness has been linked to the repertoire of brain states at various spatiotemporal scales. Anesthesia is thought to modify consciousness by altering information integration in cortical and thalamocortical circuits. At a mesoscopic scale, neuronal populations in the cortex form synchronized ensembles whose characteristics are presumably state-dependent but this has not been rigorously tested. In this study, spontaneous neuronal activity was recorded with 64-contact microelectrode arrays in primary visual cortex of chronically instrumented, unrestrained rats under stepwise decreasing levels of desflurane anesthesia (8%, 6%, 4%, and 2% inhaled concentrations) and wakefulness (0% concentration). Negative phases of the local field potentials formed compact, spatially contiguous activity patterns (CAPs) that were not due to chance. The number of CAPs was 120% higher in wakefulness and deep anesthesia associated with burst-suppression than at intermediate levels of consciousness. The frequency distribution of CAP sizes followed a power-law with slope -1.5 in relatively deep anesthesia (8-6%) but deviated from that at the lighter levels. Temporal variance and entropy of CAP sizes were lowest in wakefulness (76% and 24% lower at 0% than at 8% desflurane, respectively) but changed little during recovery of consciousness. CAPs categorized by K-means clustering were conserved at all anesthesia levels and wakefulness, although their proportion changed in a state-dependent manner. These observations yield new knowledge about the dynamic landscape of ongoing population activity in sensory cortex at graded levels of anesthesia. The repertoire of population activity and self-organized criticality at the mesoscopic scale do not appear to contribute to anesthetic suppression of consciousness, which may instead depend on large-scale effects, more subtle dynamic properties, or changes outside of primary sensory cortex.