The human alpha (8 -12 Hz) rhythm is one of the most prominent, robust, and widely studied attributes of ongoing cortical activity. Contrary to the prevalent notion that it simply "waxes and wanes," spontaneous alpha activity bursts erratically between two distinct modes of activity. We now establish a mechanism for this multistable phenomenon in resting-state cortical recordings by characterizing the complex dynamics of a biophysical model of macroscopic corticothalamic activity. This is achieved by studying the predicted activity of cortical and thalamic neuronal populations in this model as a function of its dynamic stability and the role of nonspecific synaptic noise. We hence find that fluctuating noisy inputs into thalamic neurons elicit spontaneous bursts between low-and high-amplitude alpha oscillations when the system is near a particular type of dynamical instability, namely a subcritical Hopf bifurcation. When the postsynaptic potentials associated with these noisy inputs are modulated by cortical feedback, the SD of power within each of these modes scale in proportion to their mean, showing remarkable concordance with empirical data. Our state-dependent corticothalamic model hence exhibits multistability and scale-invariant fluctuations-key features of resting-state cortical activity and indeed, of human perception, cognition, and behavior-thus providing a unified account of these apparently divergent phenomena.

The brain is widely assumed to be a paradigmatic example of a complex, self-organizing system. As such, it should exhibit the classic hallmarks of nonlinearity, multistability, and "nondiffusivity" (large coherent fluctuations). Surprisingly, at least at the very large scale of neocortical dynamics, there is little empirical evidence to support this, and hence most computational and methodological frameworks for healthy brain activity have proceeded very reasonably from a purely linear and diffusive perspective. By studying the temporal fluctuations of power in human resting-state electroencephalograms, we show that, although these simple properties may hold true at some temporal scales, there is strong evidence for bistability and nondiffusivity in key brain rhythms. Bistability is manifest as nonclassic bursting between high-and low-amplitude modes in the alpha rhythm. Nondiffusivity is expressed through the irregular appearance of high amplitude "extremal" events in beta rhythm power fluctuations. The statistical robustness of these observations was confirmed through comparison with Gaussian-rendered phase-randomized surrogate data. Although there is a good conceptual framework for understanding bistability in cortical dynamics, the implications of the extremal events challenge existing frameworks for understanding large-scale brain systems.

Multistability and scale-invariant fluctuations occur in a wide variety of biological organisms from bacteria to humans as well as financial, chemical and complex physical systems. Multistability refers to noise driven switches between multiple weakly stable states. Scale-invariant fluctuations arise when there is an approximately constant ratio between the mean and standard deviation of a system's fluctuations. Both are an important property of human perception, movement, decision making and computation and they occur together in the human alpha rhythm, imparting it with complex dynamical behavior. Here, we elucidate their fundamental dynamical mechanisms in a canonical model of nonlinear bifurcations under stochastic fluctuations. We find that the co-occurrence of multistability and scale-invariant fluctuations mandates two important dynamical properties: Multistability arises in the presence of a subcritical Hopf bifurcation, which generates co-existing attractors, whilst the introduction of multiplicative (state-dependent) noise ensures that as the system jumps between these attractors, fluctuations remain in constant proportion to their mean and their temporal statistics become long-tailed. The simple algebraic construction of this model affords a systematic analysis of the contribution of stochastic and nonlinear processes to cortical rhythms, complementing a recently proposed biophysical model. Similar dynamics also occur in a kinetic model of gene regulation, suggesting universality across a broad class of biological phenomena.

Variability of evoked single-trial responses despite constant input or task is a feature of large-scale brain signals recorded by fMRI. Initial evidence signified relevance of fMRI signal variability for perception and behavior. Yet the underlying intrinsic neuronal sources have not been previously substantiated. Here, we address this issue using simultaneous EEG-fMRI and real-time classification of ongoing alpharhythm states triggering visual stimulation in human subjects. We investigated whether spontaneous neuronal oscillations-as reflected in the posterior alpha rhythm-account for variability of evoked fMRI responses. Based on previous work, we specifically hypothesized linear superposition of fMRI activity related to fluctuations of ongoing alpha rhythm and a visually evoked fMRI response. We observed that spontaneous alpha-rhythm power fluctuations largely explain evoked fMRI response variance in extrastriate, thalamic, and cerebellar areas. For extrastriate areas, we confirmed the linear superposition hypothesis. We hence linked evoked fMRI response variability to an intrinsic rhythm's power fluctuations. These findings contribute to our conceptual understanding of how brain rhythms can account for trial-by-trial variability in stimulus processing.