Neural Field Theory of Corticothalamic Prediction With Control Systems Analysis

Unified theory of alpha, mu, and tau rhythms via eigenmodes of brain activity

A compact description of the frequency structure and topography of human alpha-band rhythms is obtained by use of the first four brain activity eigenmodes previously derived from corticothalamic neural field theory. Just two eigenmodes that overlap in frequency are found to reproduce the observed topography of the classical alpha rhythm for subjects with a single, occipitally concentrated alpha peak in their electroencephalograms. Alpha frequency splitting and relative amplitudes of double alpha peaks are explored analytically and numerically within this four-mode framework using eigenfunction expansion and perturbation methods. These effects are found to result primarily from the different eigenvalues and corticothalamic gains corresponding to the eigenmodes. Three modes with two non-overlapping frequencies suffice to reproduce the observed topography for subjects with a double alpha peak, where the appearance of a distinct second alpha peak requires an increase of the corticothalamic gain of higher eigenmodes relative to the first. Conversely, alpha blocking is inferred to be linked to a relatively small attention-dependent reduction of the gain of the relevant eigenmodes, whose effect is enhanced by the near-critical state of the brain and whose sign is consistent with inferences from neural field theory. The topographies and blocking of the mu and tau rhythms within the alpha-band are explained analogously via eigenmodes. Moreover, the observation of three rhythms in the alpha band is due to there being exactly three members of the first family of spatially nonuniform modes. These results thus provide a simple, unified description of alpha band rhythms and enable experimental observations of spectral structure and topography to be linked directly to theory and underlying physiology.

Neural field theory of adaptive effects on auditory evoked responses and mismatch negativity in multifrequency stimulus sequences

Physiologically based neural field theory (NFT) of the corticothalamic system, including adaptation, is used to calculate the responses evoked by trains of auditory stimuli that differ in frequency. In oddball paradigms, fully distinguishable frequencies lead to different standard (common stimulus) and deviant (rare stimulus) responses; the signal obtained by subtracting the standard response from the deviant is termed the mismatch negativity (MMN). In this analysis, deviant responses are found to correspond to unadapted cortex, whereas the part of auditory cortex that processes the standard stimuli adapts over several stimulus presentations until the final standard response form is achieved. No higher-order memory processes are invoked. In multifrequency experiments, the deviant response approaches the standard one as the deviant frequency approaches that of the standard and analytic criteria for this effect to be obtained. It is shown that these criteria can also be used to understand adaptation in random tone sequences. A method of probing MMNs and adaptation in random tone sequences is suggested to makes more use of such data.

Neural Field Theory of Evoked Response Sequences and Mismatch Negativity With Adaptation

Physiologically based neural field theory of the corticothalamic system is used to calculate the responses evoked by trains of auditory stimuli that correspond to different cortical locations via the tonotopic map. The results are shown to account for standard and deviant evoked responses to frequent and rare stimuli, respectively, in the auditory oddball paradigms widely used in human cognitive studies, and the so-called mismatch negativity between them. It also reproduces a wide range of other effects and variants, including the mechanism by which a change in standard responses relative to deviants can develop through adaptation, different responses when two deviants are presented in a row or a standard is presented after two deviants, relaxation of standard responses back to deviant form after a stimulus-free period, and more complex sequences. Some cases are identified in which adaptation does not account for the whole difference between standard and deviant responses. The results thus provide a systematic means to determine how much of the response is due to adaptation in the system comprising the primary auditory cortex and medial geniculate nucleus, and how much requires involvement of higher-level processing.

Determination of Dynamic Brain Connectivity via Spectral Analysis

Spectral analysis based on neural field theory is used to analyze dynamic connectivity via methods based on the physical eigenmodes that are the building blocks of brain dynamics. These approaches integrate over space instead of averaging over time and thereby greatly reduce or remove the temporal averaging effects, windowing artifacts, and noise at fine spatial scales that have bedeviled the analysis of dynamical functional connectivity (FC). The dependences of FC on dynamics at various timescales, and on windowing, are clarified and the results are demonstrated on simple test cases, demonstrating how modes provide directly interpretable insights that can be related to brain structure and function. It is shown that FC is dynamic even when the brain structure and effective connectivity are fixed, and that the observed patterns of FC are dominated by relatively few eigenmodes. Common artifacts introduced by statistical analyses that do not incorporate the physical nature of the brain are discussed and it is shown that these are avoided by spectral analysis using eigenmodes. Unlike most published artificially discretized “resting state networks” and other statistically-derived patterns, eigenmodes overlap, with every mode extending across the whole brain and every region participating in every mode—just like the vibrations that give rise to notes of a musical instrument. Despite this, modes are independent and do not interact in the linear limit. It is argued that for many purposes the intrinsic limitations of covariance-based FC instead favor the alternative of tracking eigenmode coefficients vs. time, which provide a compact representation that is directly related to biophysical brain dynamics.

Modal-Polar Representation of Evoked Response Potentials in Multiple Arousal States

An expansion of the corticothalamic transfer function into eigenmodes and resonant poles is used to derive a simple formula for evoked response potentials (ERPs) in various states of arousal. The transfer function corresponds to the cortical response to an external stimulus, which encodes all the information and properties of the linear system. This approach links experimental observations of resonances and characteristic timescales in brain activity with physically based neural field theory (NFT). The present work greatly simplifies the formula of the analytical ERP, and separates its spatial part (eigenmodes) from the temporal part (poles). Within this framework, calculations involve contour integrations that yield an explicit expression for ERPs. The dominant global mode is considered explicitly in more detail to study how the ERP varies with time in this mode and to illustrate the method. For each arousal state in sleep and wake, the resonances of the system are determined and it is found that five poles are sufficient to study the main dynamics of the system in waking eyes-open and eyes-closed states. Similarly, it is shown that six poles suffice to reproduce ERPs in rapid-eye movement sleep, sleep state 1, and sleep state 2 states, whereas just four poles suffice to reproduce the dynamics in slow wave sleep. Thus, six poles are sufficient to preserve the main global ERP dynamics of the system for all states of arousal. These six poles correspond to the dominant resonances of the system at slow-wave, alpha, and beta frequencies. These results provide the basis for simplified analytic treatment of brain dynamics and link observations more closely to theory.

Brain dynamics and structure-function relationships via spectral factorization and the transfer function

It is shown how the brain's linear transfer function provides a means to systematically analyze brain connectivity and dynamics, and to infer connectivity, eigenmodes, and activity measures such as spectra, evoked responses, coherence, and causality, all of which are widely used in brain monitoring. In particular, the Wilson spectral factorization algorithm is outlined and used to efficiently obtain linear transfer functions from experimental two-point correlation functions. The algorithm is tested on a series of brain-like structures of increasing complexity which include time delays, asymmetry, two-dimensionality, and complex network connectivity. These tests are used to verify the algorithm is suitable for application to brain dynamics, specify sampling requirements for experimental time series, and to verify that its runtime is short enough to obtain accurate results for systems of similar size to current experiments. The results can equally well be applied to inference of the transfer function in complex linear systems other than brains.

Whole-Brain Models to Explore Altered States of Consciousness from the Bottom Up

The scope of human consciousness includes states departing from what most of us experience as ordinary wakefulness. These altered states of consciousness constitute a prime opportunity to study how global changes in brain activity relate to different varieties of subjective experience. We consider the problem of explaining how global signatures of altered consciousness arise from the interplay between large-scale connectivity and local dynamical rules that can be traced to known properties of neural tissue. For this purpose, we advocate a research program aimed at bridging the gap between bottom-up generative models of whole-brain activity and the top-down signatures proposed by theories of consciousness. Throughout this paper, we define altered states of consciousness, discuss relevant signatures of consciousness observed in brain activity, and introduce whole-brain models to explore the biophysics of altered consciousness from the bottom-up. We discuss the potential of our proposal in view of the current state of the art, give specific examples of how this research agenda might play out, and emphasize how a systematic investigation of altered states of consciousness via bottom-up modeling may help us better understand the biophysical, informational, and dynamical underpinnings of consciousness.

Neural Field Theory of Evoked Response Potentials With Attentional Gain Dynamics

A generalized neural field model of large-scale activity in the corticothalamic system is used to predict standard evoked potentials. This model embodies local feedbacks that modulate the gains of neural activity as part of the response to incoming stimuli and thus enables both activity changes and effective connectivity changes to be calculated as parts of a generalized evoked response, and their relative contributions to be determined. The results show that incorporation of gain modulations enables a compact and physically justifiable description of the differences in gain between background-EEG and standard-ERP conditions, with the latter able to be initiated from the background state, rather than requiring distinct parameters as in earlier work. In particular, top-down gains are found to be reduced during an ERP, consistent with recent theoretical suggestions that the role of internal models is diminished in favor of external inputs when the latter change suddenly. The static-gain and modulated-gain system transfer functions are analyzed via control theory in terms of system resonances that were recently shown to implement data filtering whose gain adjustments can be interpreted as attention. These filters are shown to govern early and late features in standard evoked responses and their gain parameters are shown to be dynamically adjusted in a way that implements a form of attention. The results show that dynamically modulated resonant filters responsible for the low-frequency oscillations in an evoked potential response have different parameters than those responsible for low-frequency resting EEG responses, while both responses share similar mid- and high-frequency resonant filters. These results provide a biophysical mechanism by which oscillatory activity in the theta, alpha, and beta frequency ranges of an evoked response are modulated as reflections of attention; notably theta is enhanced and alpha suppressed during the latter parts of the ERP. Furthermore, the model enables the part of the ERP response induced by gain modulations to be estimated and interpreted in terms of attention.

Neural Field Theory of Corticothalamic Attention With Control System Analysis

Neural field theory is used to analyze attention by extending an existing model of the large-scale activity in the corticothalamic system to incorporate local feedbacks that modulate the gains of neural connectivity as part of the response to incoming stimuli. Treatment of both activity changes and connectivity changes as part of a generalized response enables generalized linear transfer functions of the combined response to be derived. These are then analyzed and interpreted via control theory in terms of stimulus-driven changes in system resonances that were recently shown to implement data filtering and prediction of the inputs. Using simple visual stimuli as a test case, it is shown that the gain response can implement attention by evaluating two main features of the stimuli: the magnitude and the rate of change, by increasing the weight placed on the rate of change in response to sudden changes, while reducing the contribution of stimuli value in tandem. These changes of filter parameters are shown to improve the prediction of the upcoming stimuli based on its recent time course. This outcome is analogous to controller-parameter tuning for performance enhancement in engineering control theory.

Neural field theory of evoked response potentials in a spherical brain geometry

Evoked response potentials (ERPs) are calculated in spherical and planar geometries using neural field theory of the corticothalamic system. The ERP is modeled as an impulse response and the resulting modal effects of spherical corticothalamic dynamics are explored, showing that results for spherical and planar geometries converge in the limit of large brain size. Cortical modal effects can lead to a double-peak structure in the ERP time series. It is found that the main difference between infinite planar geometry and spherical geometry is that the ERP peak is sharper and stronger in the spherical geometry. It is also found that the magnitude of the response decreases with increasing spatial width of the stimulus at the cortex. The peak is slightly delayed at large angles from the stimulus point, corresponding to group velocities of 6-10 m s^{-1}. Strong modal effects are found in the spherical geometry, with the lowest few modes sufficing to describe the main features of ERPs, except very near to spatially narrow stimuli.

Physical brain connectomics

Brain connectivity and structure-function relationships are analyzed from a physical perspective in place of common graph-theoretic and statistical approaches that overwhelmingly ignore the brain's physical structure and geometry. Field theory is used to define connectivity tensors in terms of bare and dressed propagators, and discretized representations are implemented that respect the physical nature and dimensionality of the quantities involved, retain the correct continuum limit, and enable diagrammatic analysis. Eigenfunction analysis is used to simultaneously characterize and probe patterns of brain connectivity and activity, in place of statistical or phenomenological patterns. Physically based measures that characterize the connectivity are then developed in coordinate and spectral domains; some of which generalize or rectify graph-theoretic measures to implement correct dimensionality and continuum limits, and some replace graph-theoretic quantities. Traditional graph-based measures are shown to be highly prone to artifacts introduced by discretization and threshold, often because essential physical constraints have not been imposed, dimensionality has not been included, and/or distinctions between scalar, vector, and tensor quantities have not been considered. The results can replace them in ways that converge correctly and measure properties of brain structure, rather than of its discretization, and thus potentially enable physical interpretation of the many phenomenological results in the literature. Geometric effects are shown to dominate in determining many brain properties and care must be taken not to interpret geometric differences as differences in intrinsic neural connectivity. The results demonstrate the need to use systematic physical methods to analyze the brain and the potential of such methods to obtain new insights from data, make new predictions for experimental test, and go beyond phenomenological classification to dynamics and mechanisms.