Neuronal activity in the brain has been traditionally separated into ongoing activity, which lacks a particular sensory stimulus or movement, and evoked activity, which is the response to a well-defined stimulus. Yet, studies have consistently shown that ongoing activity predicts stimulus response (Arieli et al., 1996; Tsodyks et al., 1999; Luczak et al., 2009) and behavioral outcome (Supèr et al., 2003; Womelsdorf et al., 2006) suggesting that both forms of activity are closely related. During the last decade, ongoing activity has been found to organize as neuronal avalanches (Beggs and Plenz, 2003) — scale-invariant activity cascades whose size, duration, waveform and inter-cascade intervals are governed by power laws (for review see Chialvo (2010) and Plenz (2012)). From microscale to macroscale organization, neuronal avalanches describe spontaneous firing of local pyramidal neuron groups in awake rodents in vivo (Bellay et al., 2015), the ongoing local field potential (LFP) and subthreshold population activity in rodents and nonhuman primates (Gireesh and Plenz, 2008; Petermann et al., 2009; Scott et al., 2014), as well as resting activity in humans observed using functional magnetic resonance imaging (fMRI) (Fraiman and Chialvo, 2012; Tagliazucchi et al., 2012; Haimovici et al., 2013), magnetoencephalography (MEG) (Palva et al., 2013; Shriki et al., 2013) and electrocorticography (ECoG) (Priesemann et al., 2014). Neuronal avalanches signify a cortical network at a critical or near-critical state (Chialvo, 2010; Plenz, 2012), which experiments (Shew et al., 2009; Gautam et al., 2015) and simulations (Beggs and Plenz, 2003; Bertschinger and Natschläger, 2004; Haldeman and Beggs, 2005; Kinouchi and Copelli, 2006; Rämö et al., 2007; Nykter et al., 2008) have shown to hold numerous advantages in information processing, such as maximum mnemonic repertoire size (Haldeman and Beggs, 2005; Shew et al., 2011; Fagerholm et al., 2016), information diversity (Nykter et al., 2008) and dynamic range (Kinouchi and Copelli, 2006; Shew et al., 2009; Gautam et al., 2015), optimized computational capabilities (Bertschinger and Natschläger, 2004), information transmission (Beggs and Plenz, 2003; Rämö et al., 2007; Fagerholm et al., 2016), sensory discrimination (Tomen et al., 2014) and learning (de Arcangelis and Herrmann, 2010). Yet, despite these potential advantages, evidence for neuronal avalanches outside the realm of ongoing activity, specifically, during behaviorally relevant motor and cognitive tasks, has been controversial.

Recent analysis of neuronal avalanches during sustained dynamic movie stimulation in the visual cortex of turtle (Shew et al., 2015) and visuo-motor tasks in human fMRI (Fagerholm et al., 2015) and EEG (Arviv et al., 2015) suggest transient dynamics that differ from scale-invariant avalanches. Similarly, neural models implied transient deviations from avalanche dynamics during strong external inputs (Millman et al., 2010; Taylor et al., 2013; Hartley et al., 2014; Stepp et al., 2015; Williams-García et al., 2014). Here, we hypothesize that transient changes in activity during evoked responses could simply reflect a change in the rate of avalanches rather than a transition to a different dynamical regime, for example a change in synchronization. Such a rate change likely requires modifications to standard avalanche analysis, as the standard was originally introduced for stationary rates of activity (Beggs and Plenz, 2003). In fact, by taking changes in activity rate into account, here we show in the behaving nonhuman primate that cortical dynamics during sensory and cognitive processing exhibit clear power law organization in line with the maintenance of neuronal avalanches. In line with these findings, our phase-synchronization analysis demonstrates that synchrony levels are maintained for ongoing and task-related activity. We expand on our experimental results using computational models and explore the sensitivity of this advanced approach to distinguish avalanche dynamics from states with changed synchrony, for example supercritical dynamics. Our results extend previous work on the dynamics of resting state activity and suggest that the optimization principles identified for critical dynamics during resting state activity extend to and sculpt cortical information processing during behavior.

In order to demonstrate avalanche maintenance during evoked activity that is independent of specific brain region or behavioral paradigm, we chose two types of tasks and two cortical areas in two adult macaque monkeys. High-density microelectrode arrays (MEA; 10 × 10; corner electrodes missing; inter electrode distance Δd = 400 µm) were chronically implanted in left prefrontal cortex (PFC) of monkey A and left premotor cortex (PM) of monkey B. We recorded extracellular unit activity (0.3–3 kHz) to confirm involvement of recording sites in behavioral tasks as well as the LFP (1–100 Hz), which is commonly used for avalanche analysis (Beggs and Plenz, 2003; Petermann et al., 2009; Scott et al., 2014).

Avalanche maintenance during a cue-triggered cognitive task

In our first behavioral condition, we asked whether scale-invariant dynamics can be found when strong transient activity changes due to cognitive processing are present. To this end, we analyzed avalanche dynamics in the dorsal-lateral prefrontal cortex (dlPFC; left hemisphere) for a visual-motor mapping task (monkey A; Figures 1–3). Monkey A applied a learned rule (>95% success rate) to two individually presented visual cues, which instructed the retrieval of food from a left feeder (left trials) or right feeder (right trials).

Figure 1

Download asset Open asset

Figure 2

Download asset Open asset

Figure 3

Download asset Open asset

It is known that the dlPFC is involved in such visual-motor mapping task (Asaad et al., 1998; Puig and Miller, 2012). Consistently, in all of 29 putative single-units recorded in the dlPFC (1 day, 1 hr recording session), unit activity was variable with transient increases or decreases in firing rate around cue on- and off-set that depended on left or right trials (Figure 1A; three examples shown). No significant change in Fano Factor (Churchland et al., 2010) was found during task execution (Figure 1B; see Materials and methods). On the other hand, the average LFP changed consistently for each electrode exhibiting two waveforms that differed between left and right trials (Figure 1C). We extracted peak times of the negative LFP (nLFP) at a threshold thr = –2 SD for further analysis (see Materials and methods), and observed a rapid ~8-fold increase in nLFP rate ~0.5 s after cue presentation, when the monkey was preparing to reach for the left feeder (Figure 1D, single electrode; Figure 1E, all electrodes). A much weaker, more delayed increase was observed for right trials. Based on the delayed increase in nLFP rate for left trials, we divided trials into a baseline period before cue-onset (BASE; –0.4–0 s), and an early (EARLY; 0–0.4 s) and late epoch (LATE; 0.4–0.8 s) after cue-onset.

To study neuronal avalanches, two 1 hr recording sessions from 2 consecutive days resulting in 322 trials (day 1: 139, day 2: 183) were combined. In short, for each trial successive nLFPs on the array were concatenated into contiguous clusters at temporal resolution Δt (Beggs and Plenz, 2003). A cluster started with the transition from an empty time bin of width Δt to a time bin with at least one nLFP on the array. The concatenation continued until a time bin with no nLFP was encountered (see also Figure 2A). The temporal resolution Δt was defined by the inverse of the average nLFP rate on the array and is not a free parameter (Beggs and Plenz, 2003). In line with standard avalanche analysis, we kept Δt fixed (‘fixed binning’) and determined its value from the full recording to be Δt = 4.75 ms for monkey A. We found that the probability density function of cluster size s, that is the number of nLFPs in a cluster, followed a power law, for both left and right trials, with slope α close to –1.5 and a cut-off for s > 96, the maximal number of MEA electrodes (Klaus et al., 2011; Yu et al., 2014) (Figure 1F, solid black and green lines for left and right trials, respectively; exponential vs. power law: ${\displaystyle p<{10}^{-5}}$; see Materials and methods). Importantly, these hallmarks of avalanche dynamics were destroyed by trial shuffling, which removes spatial trial-to-trial correlations while maintaining average rate changes (Figure 1F, broken lines; exponential vs. power law: ${\displaystyle 1-p<{10}^{-5}}$). This control demonstrates that the power law in cluster sizes, indicative of avalanche dynamics, does not simply arise from transient activity changes during behaviorally relevant periods.

By combining clusters from different behavioral epochs into one single distribution, however, transient deviations from scale-invariance could cancel each other, giving the appearance of maintained avalanche dynamics during behavior. Indeed, when separating clusters according to trial epochs BASE, EARLY and LATE, systematic differences in the corresponding avalanche raster were revealed. Specifically, the increase in nLFP rate during LATE for left trials yielded significantly larger avalanches during the period prior to reach (Figure 2B, red dots; Figure 2C, orange and purple segments indicate significantly high and low values, respectively; $p\le 0.05$; see Materials and methods). Accordingly, a preponderance of large avalanches was visible in the corresponding size distribution during LATE for left trials (Figure 3A and B, left), whereas size distributions remained close to a power law during BASE and EARLY and for right trials (Figure 3A and B; ${\displaystyle p<{10}^{-5}}$, all cases). In the case of right trials, where the transient changes in activity rate are much less pronounced, size distributions are more similar among all epochs, with a small tendency for larger avalanches during the LATE epoch (Figure 3A and B, right; ${\displaystyle p<{10}^{-5}}$, all cases).

These deviations from scale-invariance compared to BASE, on the other hand, were abated when the temporal resolution Δt to define clusters was obtained from the average nLFP rate for each trial i and epoch separately (${\displaystyle \mathrm{\Delta }{t}_i^{BASE,EARLY,LATE}}$; Figure 2D–F). We define this approach as ‘adaptive binning’. It decreases Δt during periods of high nLFP rate, while increasing Δt for low activity periods (Figure 2D). Indeed, adaptive binning obtained scale-invariant power laws in avalanche sizes that were similar for each epoch (Figure 3C,D; $p=0.06\left(0.08\right)$ and $p=0.11\left(0.1\right)$ for EARLY (LATE) compared to BASE, for left and right trials respectively; see Materials and methods), with a change in exponent consistent with previous results (Beggs and Plenz, 2003; Petermann et al., 2009; Priesemann et al., 2014) and a cut-off for s > 96 (Klaus et al., 2011; Yu et al., 2014). Comparing BASE and LATE avalanches for each trial, the bias toward larger ones in the latter epoch seen when employing fixed binning is significantly decreased after adaptive binning for left trials (Figure 3E, left), while during right trials, in which scale-free avalanches were obtained even for fixed binning, the bias does not change after adaptive binning. As can be seen in the corresponding raster plots and trial averages (Figure 2E,F), avalanche size decreased in LATE concomitant with an increase in avalanche rate (more prominently for left trials). We conclude that the apparent over-representations of large avalanches during behavioral epochs were due to a mismatched temporal resolution Δt by not accounting for systematic changes in nLFP rate during different epochs.

We further explored if nLFP rate within an epoch correlates with behavioral outcome. Indeed, reaction time, defined as the time it took the monkey to reach the reward after the visual cue was switched off, did not correlate with nLFP rate (left trials: R = 0.036, p=0.73; right trials: R = −0.039, p=0.71; p-values indicate the likelihood of obtaining a correlation at least as high by chance), avalanche rate (left trials: R = 0.041, p=0.68 for fixed Δt and R = −0.022, p=0.79 for adaptive Δt; right trials: R = −0.033, p=0.74 for fixed Δt and R = 0.012, p=0.87 for adaptive Δt) or avalanche size (left trials: R = 0.085, p=0.52 for fixed Δt and R = 0.062, p=0.66 for adaptive Δt; right trials: R = 0.058, p=0.63 for fixed Δt and R = −0.025, p=0.72 for adaptive Δt), in line with our finding of a scale-free regime regardless of rates for each trial. Given the monkey’s high performance (>95% successful trials), rates could not be analyzed for correct/incorrect trial outcome.

Avalanche maintenance during self-initiated movements

In monkey B, we studied activity changes during self-initiated movements. The monkey, without further cue or prompting, initiated the touching of a pad with her right hand, after which a reward was provided. Similar to monkey A, two 1 hr sessions were recorded in 2 consecutive days. Since results were highly consistent between different days, the data were pooled together, resulting in 663 trials analyzed (day 1: 319; day 2: 344). As the recording site was within the arm representation area within the PM, recorded activity clearly reflected task execution (Figure 4). About 200 ms before touching the pad, 18/43 (~42%) putative single-units increased in firing rate, whereas 5/43 units (~12%) decreased their firing rate around the time of touching (Figure 4A,B). Fano Factor of unit activity decreased in line with previous reports (Figure 4C) (Churchland et al., 2010). Similarly, the LFP transiently turned negative ~200 ms before touching for all electrodes (Figure 4D), together with a ~4-fold increase in the nLFP rate and a slight decrease thereafter (Figure 4E), which readily separated activity into a baseline epoch (BASE; –0.8 to –0.4 s relative to the touch), followed by a pre-touch (PRE; –0.4 to 0 s) and post-touch epoch (POST; 0 to 0.4 s) (Figure 4F).

Figure 4

Download asset Open asset

When all data across epochs was analyzed using fixed binning, that is Δt based on the full recording (Δt = 10.3 ms), avalanche sizes followed a power law with –1.5 exponent and cut-off at ~96 electrodes (Figure 4G; exponential vs. power law: ${\displaystyle p<{10}^{-5}}$), which was destroyed by trial shuffling (exponential vs. power law: ${\displaystyle 1-p<{10}^{-5}}$), again demonstrating that scale-invariant LFP clusters were based on the trial-by-trial correlations between cortical sites and did not simply arise from transient changes in nLFP rate.

Evaluation of avalanche raster plots based on epochs revealed systematic differences in the corresponding avalanche raster (Figure 5A,B). Specifically, the increase in nLFP rate during PRE yielded significantly more avalanches and of larger size during that period, whereas the converse was true during POST when nLFP rate decreased (Figure 5A; colored dots; Figure 5B; orange and purple segments respectively; $p\le 0.05$; see Materials and methods). The corresponding size distributions exhibited deviations from scale-invariance obtained for BASE (Figure 6A and B, left) with an excess of large avalanches for PRE and a deficit of large avalanches for POST, despite all distributions remaining close to power laws (exponential vs. power law: ${\displaystyle p<{10}^{-5}}$, all three epochs).

Figure 5

Download asset Open asset

Figure 6

Download asset Open asset

As in the case for the cognitive task, these systematic deviations in size distributions were abated when calculating Δt from the average nLFP rate for each trial i and epoch (${\displaystyle \mathrm{\Delta }{t}_i^{BASE,PRE,POST}}$; Figure 6A and B, right). Through adaptive binning we once again obtained scale-invariant power laws in avalanche size that were similar for each epoch ($p=0.177$, PRE vs. BASE; $p=0.05$, POST vs. BASE), with the bias towards larger avalanches during PRE being significantly reduced (Figure 6C). As can be seen in the corresponding raster plots and trial averages (Figure 5C,D), avalanche size decreased in PRE concomitant with an increase in avalanche rate confirming our findings in monkey A. During POST, under representation of large avalanches due to mismatched temporal resolution Δt can also be partially compensated for by adaptive binning.

Impact of ‘adaptive binning’ in capturing rate change vs. change in dynamics

We next illustrate the expected deviation from scale-invariance when rates change and its remediation using ‘adaptive binning’ in a simple schematic. In the hypothetical example shown in Figure 7A, we consider a baseline period with nLFP rate of 3 Hz per electrode, followed by a period with the rate increased to 9 Hz (Figure 7A, left). Fixed binning, which sets Δt according to the average rate of 6 Hz, underestimates the temporal resolution for baseline activity, whereas it overestimates Δt for the high activity period. The resulting power law size distributions will be steeper for the baseline period and shallower for the high activity period (Beggs and Plenz, 2003). The combined distribution will show an upward bend at the cross-over of the individual distributions (Figure 7A; middle, solid line) which seems to suggest a deviation from avalanche dynamics with an overabundance of large clusters similar to what would be expected for supercritical dynamics. Adaptive binning, that is calculating Δt for each period separately approximates the exponents for both power laws, resulting in a final distribution that matches a power law with an intermediate slope (Figure 7A, right) (Beggs and Plenz, 2003; Petermann et al., 2009; Priesemann et al., 2014). We note that simply using shorter time bins for all data would retain different slopes for different regimes and thus is not an alternative to adaptive binning.

Figure 7

Download asset Open asset

Although we illustrated this principle using two sequential periods with different rates, the same holds for repeat trials with different rates. In Figure 7B, we separated trials of monkey B into low- and high-rate responses based on PRE epochs (Figure 7B, left). Size distributions obtained from high-rate trials at fixed binning exhibit an overabundance of large avalanches (Figure 7B, middle, arrow). With adaptive binning, deviations from the baseline distribution decreased and all distributions approached power laws with a –1.2 exponent (Figure 7B, right). This demonstrates that trials in monkey B exhibit similar dynamics as those observed for monkey A or, more generally, that deviations from scale-invariance observed in the two experimental regimes can be explained by nLFP rate changes alone.

‘Adaptive’ thresholding as an alternative approach to ‘adaptive’ binning

Adaptive binning reduces Δt for periods of high rates, thereby reducing the number of nLFPs per time bin. Alternatively, one can increase the threshold for nLFP detection during these periods, likewise reducing the number of nLFPs encountered per time bin (see also [Petermann et al., 2009]). The reverse argument holds when periods of low rates are encountered. We tracked the mean and standard deviation of the LFP in successive windows of 100 ms width and used a moving threshold thr = mean – 2SD to identify nLFPs within each window (Figure 8A,C). Thus, activity modulations induced by the task were mainly reflected in nLFP amplitude, but not rate. This approach, which trades temporal resolution for local sensitivity, collapsed distributions from different periods (grey and orange areas in Figure 8B,C) in monkey A into a similar power law, indicative of sustained avalanche dynamics during task performance (Figure 8D,E).

Figure 8

Download asset Open asset

Our findings suggest that cortical avalanches are maintained in nonhuman primates during information processing, thus adding avalanche dynamics to the many links found between ongoing and evoked activities in the brain (Arieli et al., 1996; Tsodyks et al., 1999; Womelsdorf et al., 2006). We suggest that optimal information processing capabilities of avalanche dynamics (Beggs and Plenz, 2003; Bertschinger and Natschläger, 2004; Haldeman and Beggs, 2005; Kinouchi and Copelli, 2006; Rämö et al., 2007; Nykter et al., 2008; Shew et al., 2009; Tomen et al., 2014; Gautam et al., 2015) might guide sensory, motor and cognitive processing in the brain.

During ongoing activity, avalanche dynamics have been demonstrated in various brain areas including the primary visual (Hahn et al., 2010), primary motor (Petermann et al., 2009), somatosensory (Gireesh and Plenz, 2008), premotor and prefrontal (Yu et al., 2014) cortices, reflecting a common dynamic feature regardless of specific cortical areas. Therefore, we aimed to extract a similar rule that govern avalanches during evoked activity. To achieve this goal while minimizing the number of nonhuman primates used, we chose to record activity from different cortical sites and different behavioral tasks in two monkeys, while the common condition is that the recorded region was strongly activated by the corresponding task. Monkey A performed a visual-motor mapping task, in which the correct actions have to be applied to specific cues (Wise and Murray, 2000). This task is known to involve the dorsal-lateral prefrontal cortex (Asaad et al., 1998; Puig and Miller, 2012), in line with our finding of differential spiking as well as LFP responses for the left and right cues. Conversely, the target area in our self-initiated motor task for monkey B was the premotor cortex, which is involved in movement planning (Weinrich et al., 1984; Brasted and Wise, 2004) and self-initiated movement (Hoffstaedter et al., 2013). Consistent with previous studies, we found that, at both the single neuron and population level, the premotor cortex exhibited strong activation a few hundred milliseconds before the self-initiated touch. Our highly consistent results obtained for these different cortical areas and behavioral tasks suggest that avalanche dynamics might be a general principle governing cortical dynamics underlying various behaviors.

Our trial-shuffling control demonstrated that the power law in avalanche sizes reflects correlations that are maintained during behavior on a trial-by-trial basis and that this statistical feature did not simply arise from non-stationary rates. Recently, power laws, for example Zipf’s law, have been shown to arise when the range of event probabilities becomes enlarged, for example from underlying common correlations due to common inputs typically captured in latent variables (Schwab et al., 2014; Aitchison et al., 2016). Because trial-shuffling maintains such latent variables as captured in the average evoked response, yet, abolishes scale-invariance, our findings suggest that the scale-invariance encountered in our behavioral data reflects intrinsic correlation structure of trial-to-trial variability and not common correlations.

Originally introduced for resting or spontaneous activity (Beggs and Plenz, 2003), a fixed temporal resolution Δt assumes stationary rates, which is typically violated during sensory epochs or periods of movement when event rates change systematically. For avalanche analysis, the temporal resolution Δt is based on the inter-event interval distribution of the population, for example successive nLFP events on the array. This temporal scale Δt is not a free parameter. It further depends on the spatial scale, that is inter-electrode distance Δd, as well as local minimal event threshold thr. Importantly, all three parameters Δt, Δd and thr have been experimentally linked and shown to affect the slope of the power law in a predictable manner (Beggs and Plenz, 2003; Petermann et al., 2009). In our initial approach, that is adaptive binning, we kept the threshold thr and inter electrode distance Δd constant while changing Δt in accordance with the change in event rate, which in our experiments increased up to a factor of 10. This recovered power law statistics, which was confirmed by our simulation in which rate changes were created by increasing avalanche rate. These findings are in line with reported changes of the power law slope with changes in Δt (Beggs and Plenz, 2003; Petermann et al., 2009; Priesemann et al., 2014).

Our results highlight the importance of the principle of ‘separation of time scales’ (Bak et al., 1988; Vespignani and Zapperi, 1998; Plenz, 2012; Priesemann et al., 2014) for measuring avalanche dynamics. That is, the time scale at which cascading events unfold within individual avalanches, characterized by the avalanche duration, should be significantly shorter than the time scale at which different avalanches emerge, characterized by the inter-avalanche interval. Estimating the temporal resolution Δt from the inverse of the event rate provides a reasonable separation of time scales, because it balances two potential errors – false concatenation of separate avalanches and false separation of a single avalanche (Beggs and Plenz, 2003; Plenz, 2012).

In an alternative approach, we kept spatial resolution Δd and temporal resolution Δt fixed while changing threshold thr. This approach does not affect the power law organization, as demonstrated in nonhuman primates (Petermann et al., 2009) and human fMRI (Tagliazucchi et al., 2012). Particularly when Δt is difficult to change, changing thr might provide an alternative means to reestablish a separation of time scales, since large local events are less common than small local events. Our recovery of a power law for large local events during high activity periods demonstrates the analogous roles time and local event size play in avalanche dynamics. Our separate demonstration that synchrony does not change between BASE, EARLY/PRE and LATE/POST conditions further supports our finding that avalanches are maintained despite changes in event rates. Conversely, if thresholds are calculated individually for different temporal segments corresponding to different conditions of interest and the resulting distributions deviate from power laws along with concurrent changes in synchrony measures, then one has strong indication that the underlying dynamics deviate from criticality (Meisel et al., 2013).

It is important to emphasize that adaptive binning and thresholding are only tools to reveal the signatures of scale-free dynamics after they become blurred by increased drive and absence of separation of time-scales. They are not meant to be a perfect theoretical solution to the problem and have obvious limitations. For instance, there is no clear answer as to how large the window for calculating adaptive binning or thresholding should be. If it is too small, the inherent fluctuations of a critical system will be regarded as rate changes and therefore no power law dynamics can be observed. If it is too large, the rate changes can be averaged out inside them, and therefore the problem of the non-stationarity would remain. In the case of our experiments, we had a clear indication of when we should expect a rate change, which gave us a good parameter to choose the window size (~400 ms). Another problem arrives for extreme rate regimes. For very high drive, in which Δt would be much shorter than the time it takes for spikes to propagate from a neuron to its neighbors, it is unlikely that a power law can be recovered since activity becomes uncorrelated at that temporal resolution. On the other extreme, if there is almost no activity to be observed, a very large bin may require unrealistic long recordings to produce power laws. These two extreme cases may explain some of the residual error we found for left trials during PRE for monkey A (slightly larger avalanches even after adaptive binning) and during POST for monkey B (slightly smaller avalanches even after adaptive binning).

Another important point is that the advantages of criticality within a certain area allow it to optimally process information, regardless of the nature of the input received by that area. In fact, one of the advantages of criticality is maximal dynamic range (Kinouchi and Copelli, 2006; Shew et al., 2009; Gautam et al., 2015), which means that the downstream network which employs criticality will be able to respond properly for a larger range of different inputs. Nevertheless, it is possible downstream neurons employ a sort of adaptive binning mechanism. It has been suggested in simulations that the activity of an upstream network can influence how individual downstream neurons perform spatial and temporal integration at the synaptic level (Bernander et al., 1991; Rapp et al., 1992). The more input a cell receives, the more synchronized these incoming spikes have to be in order to produce a spike on the post-synaptic cell, effectively increasing the temporal resolution by which this downstream cell process its input (Bernander et al., 1991), similarly to how our proposed adaptive binning method works.

Neuronal avalanche dynamics are predominantly located in cortical layers 2/3 (Stewart and Plenz, 2006; Petermann et al., 2009) and exhibit strong non-linear components involving the interaction between cortical sites, that is the interactions of neurons or neuronal groups (Thiagarajan et al., 2010; Yu et al., 2011; Plenz, 2012). Accordingly, maintained avalanche dynamics during self-initiated and cue-related responses suggest non-linear interactions in superficial layers during cortical processing, which numerous studies have shown to be indeed the case. For example, evoked responses in superficial layers of visual cortex in the nonhuman primate have been found to be non-linear, potentially involving the local recurrent network (Williams and Shapley, 2007; Xing et al., 2012). Similarly, functional connectivity based on the interaction from neuronal firing in layer 2/3 of nonhuman primates has been found to contribute to motor coding (Hatsopoulos et al., 1998) and to be as robust or even superior over tuning curves in predicting motor outcome (Stevenson et al., 2012). With respect to their spatiotemporal spread, avalanches reveal a fractal dimension when analyzed individually (Plenz, 2012), with nearest-neighbor relationships in the average (Gireesh and Plenz, 2008; Yu et al., 2014), in line with compact spatiotemporal spreading of average evoked responses reported, for example, for premotor cortex (Rubino et al., 2006).

Our current findings may retroactively explain previous reports that found deviations from avalanche dynamics during stimulus presentation. Arviv et al. (2015) found power law avalanche size distribution when considering the full 1 s window for a visual detection task in human MEG recordings. This is consistent with our results when considering the entire period of task performance (~2 s). They also observed a trend toward supercritical dynamics during transient activity periods, just as we found predominantly large avalanches during high rate periods. Our analysis, though, shows adaptive binning as one potential approach to distinguish true supercritical dynamics from avalanche dynamics. Fagerholm et al. (2015) reported subcritical dynamics during focused attention in human EEG, in line with the observation of reduced correlations during such episodes (Cohen and Maunsell, 2009; Mitchell et al., 2009; Harris and Thiele, 2011). However, changes in activity levels were not reported in that study, and size distributions were evaluated outside their cut-off (Yu et al., 2014). Finally, early stimulus responses in the turtle’s visual system (Shew et al., 2015) have been reported to deviate from avalanches, a finding which might reflect the statistics of the stimulus used rather than intrinsic cortical dynamics.

Numerical simulations of systems under external drive or with non-zero spontaneous neuronal activity can never be truly critical, since the quiescent phase (and with it the transition connecting it to an active phase) disappears under these conditions (Taylor et al., 2013; Hartley et al., 2014; Williams-García et al., 2014). This effect is increasingly pronounced in simulations with small network size and large activity drive. Given the relatively small change in firing rate observed during behavior and the large size of the cortical network, our findings are in line with simulations that show weakly driven critical system to exhibit power law behavior when the temporal resolution is matched to the activity rate (Hartley et al., 2014).

In conclusion, our findings provide strong evidence that the scale-free organization of neuronal avalanches is maintained during evoked responses in premotor and prefrontal cortex, despite systematic fluctuations in firing rates. We therefore suggest that neuronal information processing during tasks might capitalize on various functional benefits shown for critical dynamics. This calls for further investigation into how critical dynamics may explain the execution of specific brain functions. Of equal importance, our results also shed new, methodological light on how to identify true deviations from avalanche dynamics under pathological conditions.

1. 1. Aitchison L
   2. Corradi N
   3. Latham PE

   (2016) Zipf's law arises naturally when there are underlying, unobserved variables

   PLOS Computational Biology 12:e1005110.

   https://doi.org/10.1371/journal.pcbi.1005110

   • PubMed
   • Google Scholar
2. 1. Arieli A
   2. Sterkin A
   3. Grinvald A
   4. Aertsen A

   (1996) Dynamics of ongoing activity: explanation of the large variability in evoked cortical responses

   Science 273:1868–1871.

   https://doi.org/10.1126/science.273.5283.1868

   • PubMed
   • Google Scholar
3. 1. Arviv O
   2. Goldstein A
   3. Shriki O

   (2015) Near-critical dynamics in stimulus-evoked activity of the human brain and its relation to spontaneous resting-state activity

   Journal of Neuroscience 35:13927–13942.

   https://doi.org/10.1523/JNEUROSCI.0477-15.2015

   • PubMed
   • Google Scholar
4. 1. Asaad WF
   2. Rainer G
   3. Miller EK

   (1998) Neural activity in the primate prefrontal cortex during associative learning

   Neuron 21:1399–1407.

   https://doi.org/10.1016/S0896-6273(00)80658-3

   • PubMed
   • Google Scholar
5. 1. Bak P
   2. Tang C
   3. Wiesenfeld K

   (1988) Self-organized criticality

   Physical Review A 38:364–374.

   https://doi.org/10.1103/PhysRevA.38.364

   • PubMed
   • Google Scholar
6. 1. Beggs JM
   2. Plenz D

   (2003)

   Neuronal avalanches in neocortical circuits

   Journal of Neuroscience 23:11167–11177.

   • PubMed
   • Google Scholar
7. 1. Bellay T
   2. Klaus A
   3. Seshadri S
   4. Plenz D

   (2015) Irregular spiking of pyramidal neurons organizes as scale-invariant neuronal avalanches in the awake state

   eLife 4:e07224.

   https://doi.org/10.7554/eLife.07224

   • PubMed
   • Google Scholar
8. 1. Bernander O
   2. Douglas RJ
   3. Martin KA
   4. Koch C

   (1991) Synaptic background activity influences spatiotemporal integration in single pyramidal cells

   PNAS 88:11569–11573.

   https://doi.org/10.1073/pnas.88.24.11569

   • PubMed
   • Google Scholar
9. 1. Bertschinger N
   2. Natschläger T

   (2004) Real-time computation at the edge of chaos in recurrent neural networks

   Neural Computation 16:1413–1436.

   https://doi.org/10.1162/089976604323057443

   • PubMed
   • Google Scholar
10. 1. Brasted PJ
    2. Wise SP

    (2004) Comparison of learning-related neuronal activity in the dorsal premotor cortex and striatum

    European Journal of Neuroscience 19:721–740.

    https://doi.org/10.1111/j.0953-816X.2003.03181.x

    • PubMed
    • Google Scholar
11. 1. Chialvo DR

    (2010) Emergent complex neural dynamics

    Nature Physics 6:744–750.

    https://doi.org/10.1038/nphys1803

    • Google Scholar
12. 1. Churchland MM
    2. Yu BM
    3. Cunningham JP
    4. Sugrue LP
    5. Cohen MR
    6. Corrado GS
    7. Newsome WT
    8. Clark AM
    9. Hosseini P
    10. Scott BB
    11. Bradley DC
    12. Smith MA
    13. Kohn A
    14. Movshon JA
    15. Armstrong KM
    16. Moore T
    17. Chang SW
    18. Snyder LH
    19. Lisberger SG
    20. Priebe NJ
    21. Finn IM
    22. Ferster D
    23. Ryu SI
    24. Santhanam G
    25. Sahani M
    26. Shenoy KV

    (2010) Stimulus onset quenches neural variability: a widespread cortical phenomenon

    Nature Neuroscience 13:369–378.

    https://doi.org/10.1038/nn.2501

    • PubMed
    • Google Scholar
13. 1. Clauset A
    2. Shalizi CR
    3. Newman MEJ

    (2009) Power-law distributions in empirical data

    SIAM Review 51:661–703.

    https://doi.org/10.1137/070710111

    • Google Scholar
14. 1. Cohen MR
    2. Maunsell JH

    (2009) Attention improves performance primarily by reducing interneuronal correlations

    Nature Neuroscience 12:1594–1600.

    https://doi.org/10.1038/nn.2439

    • PubMed
    • Google Scholar
15. 1. de Arcangelis L
    2. Herrmann HJ

    (2010) Learning as a phenomenon occurring in a critical state

    PNAS 107:3977–3981.

    https://doi.org/10.1073/pnas.0912289107

    • PubMed
    • Google Scholar
16. 1. Efron B

    (1979) Bootstrap methods: another look at the jackknife

    The Annals of Statistics 7:1–26.

    https://doi.org/10.1214/aos/1176344552

    • Google Scholar
17. 1. Fagerholm ED
    2. Lorenz R
    3. Scott G
    4. Dinov M
    5. Hellyer PJ
    6. Mirzaei N
    7. Leeson C
    8. Carmichael DW
    9. Sharp DJ
    10. Shew WL
    11. Leech R

    (2015) Cascades and cognitive state: focused attention incurs subcritical dynamics

    Journal of Neuroscience 35:4626–4634.

    https://doi.org/10.1523/JNEUROSCI.3694-14.2015

    • PubMed
    • Google Scholar
18. 1. Fagerholm ED
    2. Scott G
    3. Shew WL
    4. Song C
    5. Leech R
    6. Knöpfel T
    7. Sharp DJ

    (2016) Cortical entropy, mutual information and scale-free dynamics in waking mice

    Cerebral Cortex 26:3945–3952.

    https://doi.org/10.1093/cercor/bhw200

    • PubMed
    • Google Scholar
19. 1. Fraiman D
    2. Chialvo DR

    (2012) What kind of noise is brain noise: anomalous scaling behavior of the resting brain activity fluctuations

    Frontiers in Physiology 3:307.

    https://doi.org/10.3389/fphys.2012.00307

    • PubMed
    • Google Scholar
20. 1. Gautam SH
    2. Hoang TT
    3. McClanahan K
    4. Grady SK
    5. Shew WL

    (2015) Maximizing sensory dynamic range by tuning the cortical state to criticality

    PLoS Computational Biology 11:e1004576.

    https://doi.org/10.1371/journal.pcbi.1004576

    • PubMed
    • Google Scholar
21. 1. Gerstein GL
    2. Perkel DH

    (1969) Simultaneously recorded trains of action potentials: analysis and functional interpretation

    Science 164:828–830.

    https://doi.org/10.1126/science.164.3881.828

    • PubMed
    • Google Scholar
22. 1. Gerstein GL
    2. Perkel DH

    (1972) Mutual temporal relationships among neuronal spike trains. statistical techniques for display and analysis

    Biophysical journal 12:453–473.

    https://doi.org/10.1016/S0006-3495(72)86097-1

    • PubMed
    • Google Scholar
23. 1. Gireesh ED
    2. Plenz D

    (2008) Neuronal avalanches organize as nested theta- and beta/gamma-oscillations during development of cortical layer 2/3

    PNAS 105:7576–7581.

    https://doi.org/10.1073/pnas.0800537105

    • PubMed
    • Google Scholar
24. 1. Hahn G
    2. Petermann T
    3. Havenith MN
    4. Yu S
    5. Singer W
    6. Plenz D
    7. Nikolic D

    (2010) Neuronal avalanches in spontaneous activity in vivo

    Journal of Neurophysiology 104:3312–3322.

    https://doi.org/10.1152/jn.00953.2009

    • PubMed
    • Google Scholar
25. 1. Haimovici A
    2. Tagliazucchi E
    3. Balenzuela P
    4. Chialvo DR

    (2013) Brain organization into resting state networks emerges at criticality on a model of the human connectome

    Physical Review Letters 110:178101.

    https://doi.org/10.1103/PhysRevLett.110.178101

    • PubMed
    • Google Scholar
26. 1. Haldeman C
    2. Beggs JM

    (2005) Critical branching captures activity in living neural networks and maximizes the number of metastable States

    Physical Review Letters 94:058101.

    https://doi.org/10.1103/PhysRevLett.94.058101

    • PubMed
    • Google Scholar
27. 1. Harris KD
    2. Thiele A

    (2011) Cortical state and attention

    Nature Reviews Neuroscience 12:509–523.

    https://doi.org/10.1038/nrn3084

    • PubMed
    • Google Scholar
28. 1. Hartley C
    2. Taylor TJ
    3. Kiss IZ
    4. Farmer SF
    5. Berthouze L

    (2014) Identification of criticality in neuronal avalanches: ii. a theoretical and empirical investigation of the driven case

    The Journal of Mathematical Neuroscience 4:9.

    https://doi.org/10.1186/2190-8567-4-9

    • PubMed
    • Google Scholar
29. 1. Hatsopoulos NG
    2. Ojakangas CL
    3. Paninski L
    4. Donoghue JP

    (1998) Information about movement direction obtained from synchronous activity of motor cortical neurons

    PNAS 95:15706–15711.

    https://doi.org/10.1073/pnas.95.26.15706

    • PubMed
    • Google Scholar
30. 1. Hoffstaedter F
    2. Grefkes C
    3. Zilles K
    4. Eickhoff SB

    (2013) The "what" and "when" of self-initiated movements

    Cerebral cortex 23:520–530.

    https://doi.org/10.1093/cercor/bhr391

    • PubMed
    • Google Scholar
31. 1. Kinouchi O
    2. Copelli M

    (2006) Optimal dynamical range of excitable networks at criticality

    Nature Physics 2:348–351.

    https://doi.org/10.1038/nphys289

    • Google Scholar
32. 1. Kitzbichler MG
    2. Smith ML
    3. Christensen SR
    4. Bullmore E

    (2009) Broadband criticality of human brain network synchronization

    PLoS Computational Biology 5:e1000314.

    https://doi.org/10.1371/journal.pcbi.1000314

    • Google Scholar
33. 1. Klaus A
    2. Yu S
    3. Plenz D

    (2011) Statistical analyses support power law distributions found in neuronal avalanches

    PLoS One 6:e19779.

    https://doi.org/10.1371/journal.pone.0019779

    • PubMed
    • Google Scholar
34. 1. Larremore DB
    2. Shew WL
    3. Ott E
    4. Restrepo JG

    (2011) Effects of network topology, transmission delays, and refractoriness on the response of coupled excitable systems to a stochastic stimulus

    Chaos: An Interdisciplinary Journal of Nonlinear Science 21:025117.

    https://doi.org/10.1063/1.3600760

    • PubMed
    • Google Scholar
35. 1. Luczak A
    2. Barthó P
    3. Harris KD

    (2009) Spontaneous events outline the realm of possible sensory responses in neocortical populations

    Neuron 62:413–425.

    https://doi.org/10.1016/j.neuron.2009.03.014

    • PubMed
    • Google Scholar
36. 1. Meisel C
    2. Olbrich E
    3. Shriki O
    4. Achermann P

    (2013) Fading signatures of critical brain dynamics during sustained wakefulness in humans

    Journal of Neuroscience 33:17363–17372.

    https://doi.org/10.1523/JNEUROSCI.1516-13.2013

    • PubMed
    • Google Scholar
37. 1. Meisel C
    2. Plenz D
    3. Schulze-Bonhage A
    4. Reichmann H

    (2016) Quantifying antiepileptic drug effects using intrinsic excitability measures

    Epilepsia 57:e210–e215.

    https://doi.org/10.1111/epi.13517

    • PubMed
    • Google Scholar
38. 1. Millman D
    2. Mihalas S
    3. Kirkwood A
    4. Niebur E

    (2010) Self-organized criticality occurs in non-conservative neuronal networks during Up states

    Nature Physics 6:801–805.

    https://doi.org/10.1038/nphys1757

    • PubMed
    • Google Scholar
39. 1. Mitchell JF
    2. Sundberg KA
    3. Reynolds JH

    (2009) Spatial attention decorrelates intrinsic activity fluctuations in macaque area V4

    Neuron 63:879–888.

    https://doi.org/10.1016/j.neuron.2009.09.013

    • PubMed
    • Google Scholar
40. 1. Mitz AR
    2. Boring SA
    3. Wise SP
    4. Lebedev MA

    (2001) A novel food-delivery device for neurophysiological and neuropsychological studies in monkeys

    Journal of Neuroscience Methods 109:129–135.

    https://doi.org/10.1016/S0165-0270(01)00406-X

    • PubMed
    • Google Scholar
41. 1. Nykter M
    2. Price ND
    3. Larjo A
    4. Aho T
    5. Kauffman SA
    6. Yli-Harja O
    7. Shmulevich I

    (2008) Critical networks exhibit maximal information diversity in structure-dynamics relationships

    Physical Review Letters 100:058702.

    https://doi.org/10.1103/PhysRevLett.100.058702

    • PubMed
    • Google Scholar
42. 1. Palva JM
    2. Zhigalov A
    3. Hirvonen J
    4. Korhonen O
    5. Linkenkaer-Hansen K
    6. Palva S

    (2013) Neuronal long-range temporal correlations and avalanche dynamics are correlated with behavioral scaling laws

    PNAS 110:3585–3590.

    https://doi.org/10.1073/pnas.1216855110

    • PubMed
    • Google Scholar
43. 1. Petermann T
    2. Thiagarajan TC
    3. Lebedev MA
    4. Nicolelis MA
    5. Chialvo DR
    6. Plenz D

    (2009) Spontaneous cortical activity in awake monkeys composed of neuronal avalanches

    PNAS 106:15921–15926.

    https://doi.org/10.1073/pnas.0904089106

    • PubMed
    • Google Scholar
44. 1. Plenz D

    (2012) Neuronal avalanches and coherence potentials

    The European Physical Journal Special Topics 205:259–301.

    https://doi.org/10.1140/epjst/e2012-01575-5

    • Google Scholar
45. 1. Priesemann V
    2. Wibral M
    3. Valderrama M
    4. Pröpper R
    5. Le Van Quyen M
    6. Geisel T
    7. Triesch J
    8. Nikolić D
    9. Munk MH

    (2014) Spike avalanches in vivo suggest a driven, slightly subcritical brain state

    Frontiers in Systems Neuroscience 8:108.

    https://doi.org/10.3389/fnsys.2014.00108

    • PubMed
    • Google Scholar
46. 1. Puig MV
    2. Miller EK

    (2012) The role of prefrontal dopamine D1 receptors in the neural mechanisms of associative learning

    Neuron 74:874–886.

    https://doi.org/10.1016/j.neuron.2012.04.018

    • PubMed
    • Google Scholar
47. 1. Rapp M
    2. Yarom Y
    3. Segev I

    (1992) The impact of parallel fiber background activity on the cable properties of cerebellar purkinje cells

    Neural Computation 4:518–533.

    https://doi.org/10.1162/neco.1992.4.4.518

    • Google Scholar
48. 1. Rämö P
    2. Kauffman S
    3. Kesseli J
    4. Yli-Harja O

    (2007) Measures for information propagation in boolean networks

    Physica D: Nonlinear Phenomena 227:100–104.

    https://doi.org/10.1016/j.physd.2006.12.005

    • Google Scholar
49. 1. Ribeiro TL
    2. Ribeiro S
    3. Belchior H
    4. Caixeta F
    5. Copelli M

    (2014) Undersampled critical branching processes on small-world and random networks fail to reproduce the statistics of spike avalanches

    PLoS One 9:e94992.

    https://doi.org/10.1371/journal.pone.0094992

    • PubMed
    • Google Scholar
50. 1. Rubino D
    2. Robbins KA
    3. Hatsopoulos NG

    (2006) Propagating waves mediate information transfer in the motor cortex

    Nature Neuroscience 9:1549–1557.

    https://doi.org/10.1038/nn1802

    • PubMed
    • Google Scholar
51. 1. Schwab DJ
    2. Nemenman I
    3. Mehta P

    (2014) Zipf's law and criticality in multivariate data without fine-tuning

    Physical Review Letters 113:068102.

    https://doi.org/10.1103/PhysRevLett.113.068102

    • PubMed
    • Google Scholar
52. 1. Scott G
    2. Fagerholm ED
    3. Mutoh H
    4. Leech R
    5. Sharp DJ
    6. Shew WL
    7. Knöpfel T

    (2014) Voltage imaging of waking mouse cortex reveals emergence of critical neuronal dynamics

    Journal of Neuroscience 34:16611–16620.

    https://doi.org/10.1523/JNEUROSCI.3474-14.2014

    • PubMed
    • Google Scholar
53. 1. Shew WL
    2. Clawson WP
    3. Pobst J
    4. Karimipanah Y
    5. Wright NC
    6. Wessel R

    (2015) Adaptation to sensory input tunes visual cortex to criticality

    Nature Physics 11:659–663.

    https://doi.org/10.1038/nphys3370

    • Google Scholar
54. 1. Shew WL
    2. Yang H
    3. Petermann T
    4. Roy R
    5. Plenz D

    (2009) Neuronal avalanches imply maximum dynamic range in cortical networks at criticality

    Journal of Neuroscience 29:15595–15600.

    https://doi.org/10.1523/JNEUROSCI.3864-09.2009

    • PubMed
    • Google Scholar
55. 1. Shew WL
    2. Yang H
    3. Yu S
    4. Roy R
    5. Plenz D

    (2011) Information capacity and transmission are maximized in balanced cortical networks with neuronal avalanches

    Journal of Neuroscience 31:55–63.

    https://doi.org/10.1523/JNEUROSCI.4637-10.2011

    • PubMed
    • Google Scholar
56. 1. Shriki O
    2. Alstott J
    3. Carver F
    4. Holroyd T
    5. Henson RN
    6. Smith ML
    7. Coppola R
    8. Bullmore E
    9. Plenz D

    (2013) Neuronal avalanches in the resting MEG of the human brain

    Journal of Neuroscience 33:7079–7090.

    https://doi.org/10.1523/JNEUROSCI.4286-12.2013

    • PubMed
    • Google Scholar
57. 1. Stepp N
    2. Plenz D
    3. Srinivasa N

    (2015) Synaptic plasticity enables adaptive self-tuning critical networks

    PLoS Computational Biology 11:e1004043.

    https://doi.org/10.1371/journal.pcbi.1004043

    • PubMed
    • Google Scholar
58. 1. Stevenson IH
    2. London BM
    3. Oby ER
    4. Sachs NA
    5. Reimer J
    6. Englitz B
    7. David SV
    8. Shamma SA
    9. Blanche TJ
    10. Mizuseki K
    11. Zandvakili A
    12. Hatsopoulos NG
    13. Miller LE
    14. Kording KP

    (2012) Functional connectivity and tuning curves in populations of simultaneously recorded neurons

    PLoS Computational Biology 8:e1002775.

    https://doi.org/10.1371/journal.pcbi.1002775

    • PubMed
    • Google Scholar
59. 1. Stewart CV
    2. Plenz D

    (2006) Inverted-U profile of dopamine-NMDA-mediated spontaneous avalanche recurrence in superficial layers of rat prefrontal cortex

    Journal of Neuroscience 26:8148–8159.

    https://doi.org/10.1523/JNEUROSCI.0723-06.2006

    • PubMed
    • Google Scholar
60. 1. Supèr H
    2. van der Togt C
    3. Spekreijse H
    4. Lamme VA

    (2003)

    Internal state of monkey primary visual cortex (V1) predicts figure-ground perception

    Journal of Neuroscience 23:3407–3414.

    • PubMed
    • Google Scholar
61. 1. Tagliazucchi E
    2. Balenzuela P
    3. Fraiman D
    4. Chialvo DR

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

    Frontiers in Physiology 3:15.

    https://doi.org/10.3389/fphys.2012.00015

    • PubMed
    • Google Scholar
62. 1. Taylor TJ
    2. Hartley C
    3. Simon PL
    4. Kiss IZ
    5. Berthouze L

    (2013) Identification of criticality in neuronal avalanches: I. a theoretical investigation of the non-driven case

    The Journal of Mathematical Neuroscience 3:5.

    https://doi.org/10.1186/2190-8567-3-5

    • PubMed
    • Google Scholar
63. 1. Thiagarajan TC
    2. Lebedev MA
    3. Nicolelis MA
    4. Plenz D

    (2010) Coherence potentials: loss-less, all-or-none network events in the cortex

    PLoS Biology 8:e1000278.

    https://doi.org/10.1371/journal.pbio.1000278

    • PubMed
    • Google Scholar
64. 1. Tomen N
    2. Rotermund D
    3. Ernst U

    (2014) Marginally subcritical dynamics explain enhanced stimulus discriminability under attention

    Frontiers in Systems Neuroscience 8:151.

    https://doi.org/10.3389/fnsys.2014.00151

    • PubMed
    • Google Scholar
65. 1. Tsodyks M
    2. Kenet T
    3. Grinvald A
    4. Arieli A

    (1999) Linking spontaneous activity of single cortical neurons and the underlying functional architecture

    Science 286:1943–1946.

    https://doi.org/10.1126/science.286.5446.1943

    • PubMed
    • Google Scholar
66. 1. Vespignani A
    2. Zapperi S

    (1998) How self-organized criticality works: A unified mean-field picture

    Physical Review E 57:6345–6362.

    https://doi.org/10.1103/PhysRevE.57.6345

    • Google Scholar
67. 1. Weinrich M
    2. Wise SP
    3. Mauritz KH

    (1984) A neurophysiological study of the premotor cortex in the rhesus monkey

    Brain 107:385–414.

    https://doi.org/10.1093/brain/107.2.385

    • PubMed
    • Google Scholar
68. 1. Williams PE
    2. Shapley RM

    (2007) A dynamic nonlinearity and spatial phase specificity in macaque V1 neurons

    Journal of Neuroscience 27:5706–5718.

    https://doi.org/10.1523/JNEUROSCI.4743-06.2007

    • PubMed
    • Google Scholar
69. 1. Williams-García RV
    2. Moore M
    3. Beggs JM
    4. Ortiz G

    (2014) Quasicritical brain dynamics on a nonequilibrium Widom line

    Physical Review E 90:062714.

    https://doi.org/10.1103/PhysRevE.90.062714

    • PubMed
    • Google Scholar
70. 1. Wise SP
    2. Murray EA

    (2000) Arbitrary associations between antecedents and actions

    Trends in Neurosciences 23:271–276.

    https://doi.org/10.1016/S0166-2236(00)01570-8

    • PubMed
    • Google Scholar
71. 1. Womelsdorf T
    2. Fries P
    3. Mitra PP
    4. Desimone R

    (2006) Gamma-band synchronization in visual cortex predicts speed of change detection

    Nature 439:733–736.

    https://doi.org/10.1038/nature04258

    • PubMed
    • Google Scholar
72. 1. Xing D
    2. Yeh CI
    3. Burns S
    4. Shapley RM

    (2012) Laminar analysis of visually evoked activity in the primary visual cortex

    PNAS 109:13871–13876.

    https://doi.org/10.1073/pnas.1201478109

    • PubMed
    • Google Scholar
73. 1. Yang H
    2. Shew WL
    3. Roy R
    4. Plenz D

    (2012) Maximal variability of phase synchrony in cortical networks with neuronal avalanches

    Journal of Neuroscience 32:1061–1072.

    https://doi.org/10.1523/JNEUROSCI.2771-11.2012

    • PubMed
    • Google Scholar
74. 1. Yu S
    2. Klaus A
    3. Yang H
    4. Plenz D

    (2014) Scale-invariant neuronal avalanche dynamics and the cut-off in size distributions

    PLoS One 9:e99761.

    https://doi.org/10.1371/journal.pone.0099761

    • PubMed
    • Google Scholar
75. 1. Yu S
    2. Yang H
    3. Nakahara H
    4. Santos GS
    5. Nikolić D
    6. Plenz D

    (2011) Higher-order interactions characterized in cortical activity

    Journal of Neuroscience 31:17514–17526.

    https://doi.org/10.1523/JNEUROSCI.3127-11.2011

    • PubMed
    • Google Scholar
76. 1. Zapperi S
    2. Lauritsen KB
    3. Stanley HE

    (1995) Self-organized branching processes: mean-field theory for avalanches

    Physical Review Letters 75:4071–4074.

    https://doi.org/10.1103/PhysRevLett.75.4071

    • PubMed
    • Google Scholar

Author details

1. Shan Yu

   Section on Critical Brain Dynamics, National Institute of Mental Health, Bethesda, United States

   Present address

   Institute of Automation, and Center for Excellence in Brain Science and Intelligence Technology, Chinese Academy of Sciences, Beijing, China

   Contribution

   Conceptualization, Data curation, Software, Formal analysis, Supervision, Validation, Investigation, Visualization, Methodology, Writing—original draft, Writing—review and editing

   Contributed equally with

   Tiago L Ribeiro

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-9008-6658

2. Tiago L Ribeiro

   Section on Critical Brain Dynamics, National Institute of Mental Health, Bethesda, United States

   Contribution

   Software, Formal analysis, Investigation, Visualization, Methodology, Writing—original draft, Writing—review and editing

   Contributed equally with

   Shan Yu

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0003-3195-9284

3. Christian Meisel

   Section on Critical Brain Dynamics, National Institute of Mental Health, Bethesda, United States

   Contribution

   Software, Formal analysis, Visualization, Methodology, Writing—review and editing

   Competing interests

   No competing interests declared

4. Samantha Chou

   Section on Critical Brain Dynamics, National Institute of Mental Health, Bethesda, United States

   Contribution

   Data curation, Writing—review and editing

   Competing interests

   No competing interests declared

5. Andrew Mitz

   Laboratory of Neuropsychology, National Institute of Mental Health, Bethesda, United States

   Contribution

   Resources, Data curation, Methodology

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-8045-0970

6. Richard Saunders

   Laboratory of Neuropsychology, National Institute of Mental Health, Bethesda, United States

   Contribution

   Resources, Data curation, Methodology

   Competing interests

   No competing interests declared

7. Dietmar Plenz

   Section on Critical Brain Dynamics, National Institute of Mental Health, Bethesda, United States

   Contribution

   Conceptualization, Resources, Software, Supervision, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing—original draft, Project administration, Writing—review and editing

   For correspondence

   plenzd@mail.nih.gov

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-0008-3657

• 1,425

  views

• 232

  downloads

• 59

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.