Waves, bumps, and patterns in neural field theories

A coupled neural field model for the standard consolidation theory

The standard consolidation theory states that short-term memories located in the hippocampus enable the consolidation of long-term memories in the neocortex. In other words, the neocortex slowly learns long-term memories with a transient support of the hippocampus that quickly learns unstable memories. However, it is not clear yet what could be the neurobiological mechanisms underlying these differences in learning rates and memory time-scales. Here, we propose a novel modeling approach of the standard consolidation theory, that focuses on its potential neurobiological mechanisms. In addition to synaptic plasticity and spike frequency adaptation, our model incorporates adult neurogenesis in the dentate gyrus as well as the difference in size between the neocortex and the hippocampus, that we associate with distance-dependent synaptic plasticity. We also take into account the interconnected spatial structure of the involved brain areas, by incorporating the above neurobiological mechanisms in a coupled neural field framework, where each area is represented by a separate neural field with intra- and inter-area connections. To our knowledge, this is the first attempt to apply neural fields to this process. Using numerical simulations and mathematical analysis, we explore the short-term and long-term dynamics of the model upon alternance of phases of hippocampal replay and retrieval cue of an external input. This external input is encodable as a memory pattern in the form of a multiple bump attractor pattern in the individual neural fields. In the model, hippocampal memory patterns become encoded first, before neocortical ones, because of the smaller distances between the bumps of the hippocampal memory patterns. As a result, retrieval of the input pattern in the neocortex at short time-scales necessitates the additional input delivered by the memory pattern of the hippocampus. Neocortical memory patterns progressively consolidate at longer times, up to a point where their retrieval does not need the support of the hippocampus anymore. At longer times, perturbation of the hippocampal neural fields by neurogenesis erases the hippocampus pattern, leading to a final state where the memory pattern is exclusively evoked in the neocortex. Therefore, the dynamics of our model successfully reproduces the main features of the standard consolidation theory. This suggests that neurogenesis in the hippocampus and distance-dependent synaptic plasticity coupled to synaptic depression and spike frequency adaptation, are indeed critical neurobiological processes in memory consolidation.

Criticality and partial synchronization analysis in Wilson-Cowan and Jansen-Rit neural mass models

Synchronization is a phenomenon observed in neuronal networks involved in diverse brain activities. Neural mass models such as Wilson-Cowan (WC) and Jansen-Rit (JR) manifest synchronized states. Despite extensive research on these models over the past several decades, their potential of manifesting second-order phase transitions (SOPT) and criticality has not been sufficiently acknowledged. In this study, two networks of coupled WC and JR nodes with small-world topologies were constructed and Kuramoto order parameter (KOP) was used to quantify the amount of synchronization. In addition, we investigated the presence of SOPT using the synchronization coefficient of variation. Both networks reached high synchrony by changing the coupling weight between their nodes. Moreover, they exhibited abrupt changes in the synchronization at certain values of the control parameter not necessarily related to a phase transition. While SOPT was observed only in JR model, neither WC nor JR model showed power-law behavior. Our study further investigated the global synchronization phenomenon that is known to exist in pathological brain states, such as seizure. JR model showed global synchronization, while WC model seemed to be more suitable in producing partially synchronized patterns.

Dynamics of neural fields with exponential temporal kernel

We consider the standard neural field equation with an exponential temporal kernel. We analyze the time-independent (static) and time-dependent (dynamic) bifurcations of the equilibrium solution and the emerging spatiotemporal wave patterns. We show that an exponential temporal kernel does not allow static bifurcations such as saddle-node, pitchfork, and in particular, static Turing bifurcations. However, the exponential temporal kernel possesses the important property that it takes into account the finite memory of past activities of neurons, which Green's function does not. Through a dynamic bifurcation analysis, we give explicit bifurcation conditions. Hopf bifurcations lead to temporally non-constant, but spatially constant solutions, but Turing-Hopf bifurcations generate spatially and temporally non-constant solutions, in particular, traveling waves. Bifurcation parameters are the coefficient of the exponential temporal kernel, the transmission speed of neural signals, the time delay rate of synapses, and the ratio of excitatory to inhibitory synaptic weights.

Representing stimulus motion with waves in adaptive neural fields

Traveling waves of neural activity emerge in cortical networks both spontaneously and in response to stimuli. The spatiotemporal structure of waves can indicate the information they encode and the physiological processes that sustain them. Here, we investigate the stimulus-response relationships of traveling waves emerging in adaptive neural fields as a model of visual motion processing. Neural field equations model the activity of cortical tissue as a continuum excitable medium, and adaptive processes provide negative feedback, generating localized activity patterns. Synaptic connectivity in our model is described by an integral kernel that weakens dynamically due to activity-dependent synaptic depression, leading to marginally stable traveling fronts (with attenuated backs) or pulses of a fixed speed. Our analysis quantifies how weak stimuli shift the relative position of these waves over time, characterized by a wave response function we obtain perturbatively. Persistent and continuously visible stimuli model moving visual objects. Intermittent flashes that hop across visual space can produce the experience of smooth apparent visual motion. Entrainment of waves to both kinds of moving stimuli are well characterized by our theory and numerical simulations, providing a mechanistic description of the perception of visual motion.

Neural dynamic foundations of a theory of higher cognition: the case of grounding nested phrases

Because cognitive competences emerge in evolution and development from the sensory-motor domain, we seek a neural process account for higher cognition in which all representations are necessarily grounded in perception and action. The challenge is to understand how hallmarks of higher cognition, productivity, systematicity, and compositionality, may emerge from such a bottom-up approach. To address this challenge, we present key ideas from Dynamic Field Theory which postulates that neural populations are organized by recurrent connectivity to create stable localist representations. Dynamic instabilities enable the autonomous generation of sequences of mental states. The capacity to apply neural circuitry across broad sets of inputs that emulates the function call postulated in symbolic computation emerges through coordinate transforms implemented in neural gain fields. We show how binding localist neural representations through a shared index dimension enables conceptual structure, in which the interdependence among components of a representation is flexibly expressed. We demonstrate these principles in a neural dynamic architecture that represents and perceptually grounds nested relational and action phrases. Sequences of neural processing steps are generated autonomously to attentionally select the referenced objects and events in a manner that is sensitive to their interdependencies. This solves the problem of 2 and the massive binding problem in expressions such as "the small tree that is to the left of the lake which is to the left of the large tree". We extend earlier work by incorporating new types of grammatical constructions and a larger vocabulary. We discuss the DFT framework relative to other neural process accounts of higher cognition and assess the scope and challenges of such neural theories.

Modelling cortical network dynamics

We have investigated the theoretical constraints of the interactions between coupled cortical columns. Each cortical column consists of a set of neural populations where each population is modelled as a neural mass. The existence of semi-stable states within a cortical column is dependent on the type of interaction between the neuronal populations, i.e., the form of the synaptic kernels. Current-to-current coupling has been shown, in contrast to potential-to-current coupling, to create semi-stable states within a cortical column. The interaction between semi-stable states of the cortical columns is studied where we derive the dynamics for the collected activity. For small excitations the dynamics follow the Kuramoto model; however, in contrast to previous work we derive coupled equations between phase and amplitude dynamics with the possibility of defining connectivity as a stationary and dynamic variable. The turbulent flow of phase dynamics which occurs in networks of Kuramoto oscillators would indicate turbulent changes in dynamic connectivity for coupled cortical columns which is something that has been recorded in epileptic seizures. We used the results we derived to estimate a seizure propagation model which allowed for inversions using the Laplace assumption (Dynamic Causal Modelling). The seizure propagation model was trialed on simulated data, and future work will investigate the estimation of the connectivity matrix from empirical data. This model can be used to predict changes in seizure evolution after virtual changes in the connectivity network, something that could be of clinical use when applied to epilepsy surgical cases.

Self-Organization of Nonlinearly Coupled Neural Fluctuations Into Synergistic Population Codes

Neural activity in the brain exhibits correlated fluctuations that may strongly influence the properties of neural population coding. However, how such correlated neural fluctuations may arise from the intrinsic neural circuit dynamics and subsequently affect the computational properties of neural population activity remains poorly understood. The main difficulty lies in resolving the nonlinear coupling between correlated fluctuations with the overall dynamics of the system. In this study, we investigate the emergence of synergistic neural population codes from the intrinsic dynamics of correlated neural fluctuations in a neural circuit model capturing realistic nonlinear noise coupling of spiking neurons. We show that a rich repertoire of spatial correlation patterns naturally emerges in a bump attractor network and further reveals the dynamical regime under which the interplay between differential and noise correlations leads to synergistic codes. Moreover, we find that negative correlations may induce stable bound states between two bumps, a phenomenon previously unobserved in firing rate models. These noise-induced effects of bump attractors lead to a number of computational advantages including enhanced working memory capacity and efficient spatiotemporal multiplexing and can account for a range of cognitive and behavioral phenomena related to working memory. This study offers a dynamical approach to investigating realistic correlated neural fluctuations and insights to their roles in cortical computations.

The inhibitory control of traveling waves in cortical networks

Propagating waves of activity can be evoked and can occur spontaneously in vivo and in vitro in cerebral cortex. These waves are thought to be instrumental in the propagation of information across cortical regions and as a means to modulate the sensitivity of neurons to subsequent stimuli. In normal tissue, the waves are sparse and tightly controlled by inhibition and other negative feedback processes. However, alterations of this balance between excitation and inhibition can lead to pathological behavior such as seizure-type dynamics (with low inhibition) or failure to propagate (with high inhibition). We develop a spiking one-dimensional network of neurons to explore the reliability and control of evoked waves and compare this to a cortical slice preparation where the excitability can be pharmacologically manipulated. We show that the waves enhance sensitivity of the cortical network to stimuli in specific spatial and temporal ways. To gain further insight into the mechanisms of propagation and transitions to pathological behavior, we derive a mean-field model for the synaptic activity. We analyze the mean-field model and a piece-wise constant approximation of it and study the stability of the propagating waves as spatial and temporal properties of the inhibition are altered. We show that that the transition to seizure-like activity is gradual but that the loss of propagation is abrupt and can occur via either the loss of existence of the wave or through a loss of stability leading to complex patterns of propagation.

In vivo ephaptic coupling allows memory network formation

It is increasingly clear that memories are distributed across multiple brain areas. Such "engram complexes" are important features of memory formation and consolidation. Here, we test the hypothesis that engram complexes are formed in part by bioelectric fields that sculpt and guide the neural activity and tie together the areas that participate in engram complexes. Like the conductor of an orchestra, the fields influence each musician or neuron and orchestrate the output, the symphony. Our results use the theory of synergetics, machine learning, and data from a spatial delayed saccade task and provide evidence for in vivo ephaptic coupling in memory representations.

Multi-modal and multi-model interrogation of large-scale functional brain networks

Existing whole-brain models are generally tailored to the modelling of a particular data modality (e.g., fMRI or MEG/EEG). We propose that despite the differing aspects of neural activity each modality captures, they originate from shared network dynamics. Building on the universal principles of self-organising delay-coupled nonlinear systems, we aim to link distinct features of brain activity - captured across modalities - to the dynamics unfolding on a macroscopic structural connectome. To jointly predict connectivity, spatiotemporal and transient features of distinct signal modalities, we consider two large-scale models - the Stuart Landau and Wilson and Cowan models - which generate short-lived 40 Hz oscillations with varying levels of realism. To this end, we measure features of functional connectivity and metastable oscillatory modes (MOMs) in fMRI and MEG signals - and compare them against simulated data. We show that both models can represent MEG functional connectivity (FC), functional connectivity dynamics (FCD) and generate MOMs to a comparable degree. This is achieved by adjusting the global coupling and mean conduction time delay and, in the WC model, through the inclusion of balance between excitation and inhibition. For both models, the omission of delays dramatically decreased the performance. For fMRI, the SL model performed worse for FCD and MOMs, highlighting the importance of balanced dynamics for the emergence of spatiotemporal and transient patterns of ultra-slow dynamics. Notably, optimal working points varied across modalities and no model was able to achieve a correlation with empirical FC higher than 0.4 across modalities for the same set of parameters. Nonetheless, both displayed the emergence of FC patterns that extended beyond the constraints of the anatomical structure. Finally, we show that both models can generate MOMs with empirical-like properties such as size (number of brain regions engaging in a mode) and duration (continuous time interval during which a mode appears). Our results demonstrate the emergence of static and dynamic properties of neural activity at different timescales from networks of delay-coupled oscillators at 40 Hz. Given the higher dependence of simulated FC on the underlying structural connectivity, we suggest that mesoscale heterogeneities in neural circuitry may be critical for the emergence of parallel cross-modal functional networks and should be accounted for in future modelling endeavours.

Coherent noise enables probabilistic sequence replay in spiking neuronal networks

Animals rely on different decision strategies when faced with ambiguous or uncertain cues. Depending on the context, decisions may be biased towards events that were most frequently experienced in the past, or be more explorative. A particular type of decision making central to cognition is sequential memory recall in response to ambiguous cues. A previously developed spiking neuronal network implementation of sequence prediction and recall learns complex, high-order sequences in an unsupervised manner by local, biologically inspired plasticity rules. In response to an ambiguous cue, the model deterministically recalls the sequence shown most frequently during training. Here, we present an extension of the model enabling a range of different decision strategies. In this model, explorative behavior is generated by supplying neurons with noise. As the model relies on population encoding, uncorrelated noise averages out, and the recall dynamics remain effectively deterministic. In the presence of locally correlated noise, the averaging effect is avoided without impairing the model performance, and without the need for large noise amplitudes. We investigate two forms of correlated noise occurring in nature: shared synaptic background inputs, and random locking of the stimulus to spatiotemporal oscillations in the network activity. Depending on the noise characteristics, the network adopts various recall strategies. This study thereby provides potential mechanisms explaining how the statistics of learned sequences affect decision making, and how decision strategies can be adjusted after learning.

Simulations approaching data: cortical slow waves in inferred models of the whole hemisphere of mouse

The development of novel techniques to record wide-field brain activity enables estimation of data-driven models from thousands of recording channels and hence across large regions of cortex. These in turn improve our understanding of the modulation of brain states and the richness of traveling waves dynamics. Here, we infer data-driven models from high-resolution in-vivo recordings of mouse brain obtained from wide-field calcium imaging. We then assimilate experimental and simulated data through the characterization of the spatio-temporal features of cortical waves in experimental recordings. Inference is built in two steps: an inner loop that optimizes a mean-field model by likelihood maximization, and an outer loop that optimizes a periodic neuro-modulation via direct comparison of observables that characterize cortical slow waves. The model reproduces most of the features of the non-stationary and non-linear dynamics present in the high-resolution in-vivo recordings of the mouse brain. The proposed approach offers new methods of characterizing and understanding cortical waves for experimental and computational neuroscientists.

Global dynamics of neural mass models

Neural mass models are used to simulate cortical dynamics and to explain the electrical and magnetic fields measured using electro- and magnetoencephalography. Simulations evince a complex phase-space structure for these kinds of models; including stationary points and limit cycles and the possibility for bifurcations and transitions among different modes of activity. This complexity allows neural mass models to describe the itinerant features of brain dynamics. However, expressive, nonlinear neural mass models are often difficult to fit to empirical data without additional simplifying assumptions: e.g., that the system can be modelled as linear perturbations around a fixed point. In this study we offer a mathematical analysis of neural mass models, specifically the canonical microcircuit model, providing analytical solutions describing slow changes in the type of cortical activity, i.e. dynamical itinerancy. We derive a perturbation analysis up to second order of the phase flow, together with adiabatic approximations. This allows us to describe amplitude modulations in a relatively simple mathematical format providing analytic proof-of-principle for the existence of semi-stable states of cortical dynamics at the scale of a cortical column. This work allows for model inversion of neural mass models, not only around fixed points, but over regions of phase space that encompass transitions among semi or multi-stable states of oscillatory activity. Crucially, these theoretical results speak to model inversion in the context of multiple semi-stable brain states, such as the transition between interictal, pre-ictal and ictal activity in epilepsy.

Mixed-mode chimera states in pendula networks

We report the emergence of peculiar chimera states in networks of identical pendula with global phase-lagged coupling. The states reported include both rotating and quiescent modes, i.e., with non-zero and zero average frequencies. This kind of mixed-mode chimeras may be interpreted as images of bump states known in neuroscience in the context of modeling the working memory. We illustrate this striking phenomenon for a network of N = 100 coupled pendula, followed by a detailed description of the minimal non-trivial case of N = 3. Parameter regions for five characteristic types of the system behavior are identified, which consist of two mixed-mode chimeras with one and two rotating pendula, classical weak chimera with all three pendula rotating, synchronous rotation, and quiescent state. The network dynamics is multistable: up to four of the states can coexist in the system phase state as demonstrated through the basins of attraction. The analysis suggests that the robust mixed-mode chimera states can generically describe the complex dynamics of diverse pendula-like systems widespread in nature.

Generative Models of Brain Dynamics

This review article gives a high-level overview of the approaches across different scales of organization and levels of abstraction. The studies covered in this paper include fundamental models in computational neuroscience, nonlinear dynamics, data-driven methods, as well as emergent practices. While not all of these models span the intersection of neuroscience, AI, and system dynamics, all of them do or can work in tandem as generative models, which, as we argue, provide superior properties for the analysis of neuroscientific data. We discuss the limitations and unique dynamical traits of brain data and the complementary need for hypothesis- and data-driven modeling. By way of conclusion, we present several hybrid generative models from recent literature in scientific machine learning, which can be efficiently deployed to yield interpretable models of neural dynamics.

NNMT: Mean-Field Based Analysis Tools for Neuronal Network Models

Mean-field theory of neuronal networks has led to numerous advances in our analytical and intuitive understanding of their dynamics during the past decades. In order to make mean-field based analysis tools more accessible, we implemented an extensible, easy-to-use open-source Python toolbox that collects a variety of mean-field methods for the leaky integrate-and-fire neuron model. The Neuronal Network Mean-field Toolbox (NNMT) in its current state allows for estimating properties of large neuronal networks, such as firing rates, power spectra, and dynamical stability in mean-field and linear response approximation, without running simulations. In this article, we describe how the toolbox is implemented, show how it is used to reproduce results of previous studies, and discuss different use-cases, such as parameter space explorations, or mapping different network models. Although the initial version of the toolbox focuses on methods for leaky integrate-and-fire neurons, its structure is designed to be open and extensible. It aims to provide a platform for collecting analytical methods for neuronal network model analysis, such that the neuroscientific community can take maximal advantage of them.

Gell-Mann-Low Criticality in Neural Networks

Criticality is deeply related to optimal computational capacity. The lack of a renormalized theory of critical brain dynamics, however, so far limits insights into this form of biological information processing to mean-field results. These methods neglect a key feature of critical systems: the interaction between degrees of freedom across all length scales, required for complex nonlinear computation. We present a renormalized theory of a prototypical neural field theory, the stochastic Wilson-Cowan equation. We compute the flow of couplings, which parametrize interactions on increasing length scales. Despite similarities with the Kardar-Parisi-Zhang model, the theory is of a Gell-Mann-Low type, the archetypal form of a renormalizable quantum field theory. Here, nonlinear couplings vanish, flowing towards the Gaussian fixed point, but logarithmically slowly, thus remaining effective on most scales. We show this critical structure of interactions to implement a desirable trade-off between linearity, optimal for information storage, and nonlinearity, required for computation.

Bumps, chimera states, and Turing patterns in systems of coupled active rotators

Self-organized coherence-incoherence patterns, called chimera states, have first been reported in systems of Kuramoto oscillators. For coupled excitable units, similar patterns where coherent units are at rest are called bump states. Here, we study bumps in an array of active rotators coupled by nonlocal attraction and global repulsion. We demonstrate how they can emerge in a supercritical scenario from completely coherent Turing patterns: a single incoherent unit appears in a homoclinic bifurcation, undergoing subsequent transitions to quasiperiodic and chaotic behavior, which eventually transforms into extensive chaos with many incoherent units. We present different types of transitions and explain the formation of coherence-incoherence patterns according to the classical paradigm of short-range activation and long-range inhibition.

A data-informed mean-field approach to mapping of cortical parameter landscapes

Constraining the many biological parameters that govern cortical dynamics is computationally and conceptually difficult because of the curse of dimensionality. This paper addresses these challenges by proposing (1) a novel data-informed mean-field (MF) approach to efficiently map the parameter space of network models; and (2) an organizing principle for studying parameter space that enables the extraction biologically meaningful relations from this high-dimensional data. We illustrate these ideas using a large-scale network model of the Macaque primary visual cortex. Of the 10-20 model parameters, we identify 7 that are especially poorly constrained, and use the MF algorithm in (1) to discover the firing rate contours in this 7D parameter cube. Defining a "biologically plausible" region to consist of parameters that exhibit spontaneous Excitatory and Inhibitory firing rates compatible with experimental values, we find that this region is a slightly thickened codimension-1 submanifold. An implication of this finding is that while plausible regimes depend sensitively on parameters, they are also robust and flexible provided one compensates appropriately when parameters are varied. Our organizing principle for conceptualizing parameter dependence is to focus on certain 2D parameter planes that govern lateral inhibition: Intersecting these planes with the biologically plausible region leads to very simple geometric structures which, when suitably scaled, have a universal character independent of where the intersections are taken. In addition to elucidating the geometry of the plausible region, this invariance suggests useful approximate scaling relations. Our study offers, for the first time, a complete characterization of the set of all biologically plausible parameters for a detailed cortical model, which has been out of reach due to the high dimensionality of parameter space.

Interpolating between bumps and chimeras

A "bump" refers to a group of active neurons surrounded by quiescent ones while a "chimera" refers to a pattern in a network in which some oscillators are synchronized while the remainder are asynchronous. Both types of patterns have been studied intensively but are sometimes conflated due to their similar appearance and existence in similar types of networks. Here, we numerically study a hybrid system that linearly interpolates between a network of theta neurons that supports a bump at one extreme and a network of phase oscillators that supports a chimera at the other extreme. Using the Ott/Antonsen ansatz, we derive the equation describing the hybrid network in the limit of an infinite number of oscillators and perform bifurcation analysis on this equation. We find that neither the bump nor chimera persists over the whole range of parameters, and the hybrid system shows a variety of other states such as spatiotemporal chaos, traveling waves, and modulated traveling waves.

Modulation of the dynamical state in cortical network models

Cortical neural responses can be modulated by various factors, such as stimulus inputs and the behavior state of the animal. Understanding the circuit mechanisms underlying modulations of network dynamics is important to understand the flexibility of circuit computations. Identifying the dynamical state of a network is an important first step to predict network responses to external stimulus and top-down modulatory inputs. Models in stable or unstable dynamical regimes require different analytic tools to estimate the network responses to inputs and the structure of neural variability. In this article, I review recent cortical models of state-dependent responses and their predictions about the underlying modulatory mechanisms.

The impact of a closed-loop thalamocortical model on the spatiotemporal dynamics of cortical and thalamic traveling waves

Propagation of activity in spatially structured neuronal networks has been observed in awake, anesthetized, and sleeping brains. How these wave patterns emerge and organize across brain structures, and how network connectivity affects spatiotemporal neural activity remains unclear. Here, we develop a computational model of a two-dimensional thalamocortical network, which gives rise to emergent traveling waves similar to those observed experimentally. We illustrate how spontaneous and evoked oscillatory activity in space and time emerge using a closed-loop thalamocortical architecture, sustaining smooth waves in the cortex and staggered waves in the thalamus. We further show that intracortical and thalamocortical network connectivity, cortical excitation/inhibition balance, and thalamocortical or corticothalamic delay can independently or jointly change the spatiotemporal patterns (radial, planar and rotating waves) and characteristics (speed, direction, and frequency) of cortical and thalamic traveling waves. Computer simulations predict that increased thalamic inhibition induces slower cortical frequencies and that enhanced cortical excitation increases traveling wave speed and frequency. Overall, our results provide insight into the genesis and sustainability of thalamocortical spatiotemporal patterns, showing how simple synaptic alterations cause varied spontaneous and evoked wave patterns. Our model and simulations highlight the need for spatially spread neural recordings to uncover critical circuit mechanisms for brain functions.

A model integrating multiple processes of synchronization and coherence for information instantiation within a cortical area

What is the form of dynamic, e.g., sensory, information in the mammalian cortex? Information in the cortex is modeled as a coherence map of a mixed chimera state of synchronous, phasic, and disordered minicolumns. The theoretical model is built on neurophysiological evidence. Complex spatiotemporal information is instantiated through a system of interacting biological processes that generate a synchronized cortical area, a coherent aperture. Minicolumn elements are grouped in macrocolumns in an array analogous to a phased-array radar, modeled as an aperture, a "hole through which radiant energy flows." Coherence maps in a cortical area transform inputs from multiple sources into outputs to multiple targets, while reducing complexity and entropy. Coherent apertures can assume extremely large numbers of different information states as coherence maps, which can be communicated among apertures with corresponding very large bandwidths. The coherent aperture model incorporates considerable reported research, integrating five conceptually and mathematically independent processes: 1) a damped Kuramoto network model, 2) a pumped area field potential, 3) the gating of nearly coincident spikes, 4) the coherence of activity across cortical lamina, and 5) complex information formed through functions in macrocolumns. Biological processes and their interactions are described in equations and a functional circuit such that the mathematical pieces can be assembled the same way the neurophysiological ones are. The model can be conceptually convolved over the specifics of local cortical areas within and across species. A coherent aperture becomes a node in a graph of cortical areas with a corresponding distribution of information.

Parcels and particles: Markov blankets in the brain

At the inception of human brain mapping, two principles of functional anatomy underwrote most conceptions-and analyses-of distributed brain responses: namely, functional segregation and integration. There are currently two main approaches to characterizing functional integration. The first is a mechanistic modeling of connectomics in terms of directed effective connectivity that mediates neuronal message passing and dynamics on neuronal circuits. The second phenomenological approach usually characterizes undirected functional connectivity (i.e., measurable correlations), in terms of intrinsic brain networks, self-organized criticality, dynamical instability, and so on. This paper describes a treatment of effective connectivity that speaks to the emergence of intrinsic brain networks and critical dynamics. It is predicated on the notion of Markov blankets that play a fundamental role in the self-organization of far from equilibrium systems. Using the apparatus of the renormalization group, we show that much of the phenomenology found in network neuroscience is an emergent property of a particular partition of neuronal states, over progressively coarser scales. As such, it offers a way of linking dynamics on directed graphs to the phenomenology of intrinsic brain networks.

Traveling waves in non-local pulse-coupled networks

Traveling phase waves are commonly observed in recordings of the cerebral cortex and are believed to organize behavior across different areas of the brain. We use this as motivation to analyze a one-dimensional network of phase oscillators that are nonlocally coupled via the phase response curve (PRC) and the Dirac delta function. Existence of waves is proven and the dispersion relation is computed. Using the theory of distributions enables us to write and solve an associated stability problem. First and second order perturbation theory is applied to get analytic insight and we show that long waves are stable while short waves are unstable. We apply the results to PRCs that come from mitral neurons. We extend the results to smooth pulse-like coupling by reducing the nonlocal equation to a local one and solving the associated boundary value problem.

Graph neural fields: A framework for spatiotemporal dynamical models on the human connectome

Tools from the field of graph signal processing, in particular the graph Laplacian operator, have recently been successfully applied to the investigation of structure-function relationships in the human brain. The eigenvectors of the human connectome graph Laplacian, dubbed "connectome harmonics", have been shown to relate to the functionally relevant resting-state networks. Whole-brain modelling of brain activity combines structural connectivity with local dynamical models to provide insight into the large-scale functional organization of the human brain. In this study, we employ the graph Laplacian and its properties to define and implement a large class of neural activity models directly on the human connectome. These models, consisting of systems of stochastic integrodifferential equations on graphs, are dubbed graph neural fields, in analogy with the well-established continuous neural fields. We obtain analytic predictions for harmonic and temporal power spectra, as well as functional connectivity and coherence matrices, of graph neural fields, with a technique dubbed CHAOSS (shorthand for Connectome-Harmonic Analysis Of Spatiotemporal Spectra). Combining graph neural fields with appropriate observation models allows for estimating model parameters from experimental data as obtained from electroencephalography (EEG), magnetoencephalography (MEG), or functional magnetic resonance imaging (fMRI). As an example application, we study a stochastic Wilson-Cowan graph neural field model on a high-resolution connectome graph constructed from diffusion tensor imaging (DTI) and structural MRI data. We show that the model equilibrium fluctuations can reproduce the empirically observed harmonic power spectrum of resting-state fMRI data, and predict its functional connectivity, with a high level of detail. Graph neural fields natively allow the inclusion of important features of cortical anatomy and fast computations of observable quantities for comparison with multimodal empirical data. They thus appear particularly suitable for modelling whole-brain activity at mesoscopic scales, and opening new potential avenues for connectome-graph-based investigations of structure-function relationships.

Neural field models with transmission delays and diffusion

A neural field models the large scale behaviour of large groups of neurons. We extend previous results for these models by including a diffusion term into the neural field, which models direct, electrical connections. We extend known and prove new sun-star calculus results for delay equations to be able to include diffusion and explicitly characterise the essential spectrum. For a certain class of connectivity functions in the neural field model, we are able to compute its spectral properties and the first Lyapunov coefficient of a Hopf bifurcation. By examining a numerical example, we find that the addition of diffusion suppresses non-synchronised steady-states while favouring synchronised oscillatory modes.

Stability analysis of a neural field self-organizing map

We provide theoretical conditions guaranteeing that a self-organizing map efficiently develops representations of the input space. The study relies on a neural field model of spatiotemporal activity in area 3b of the primary somatosensory cortex. We rely on Lyapunov's theory for neural fields to derive theoretical conditions for stability. We verify the theoretical conditions by numerical experiments. The analysis highlights the key role played by the balance between excitation and inhibition of lateral synaptic coupling and the strength of synaptic gains in the formation and maintenance of self-organizing maps.

Synchronization, Stochasticity, and Phase Waves in Neuronal Networks With Spatially-Structured Connectivity

Oscillations in the beta/low gamma range (10–45 Hz) are recorded in diverse neural structures. They have successfully been modeled as sparsely synchronized oscillations arising from reciprocal interactions between randomly connected excitatory (E) pyramidal cells and local interneurons (I). The synchronization of spatially distant oscillatory spiking E–I modules has been well-studied in the rate model framework but less so for modules of spiking neurons. Here, we first show that previously proposed modifications of rate models provide a quantitative description of spiking E–I modules of Exponential Integrate-and-Fire (EIF) neurons. This allows us to analyze the dynamical regimes of sparsely synchronized oscillatory E–I modules connected by long-range excitatory interactions, for two modules, as well as for a chain of such modules. For modules with a large number of neurons (> 10⁵), we obtain results similar to previously obtained ones based on the classic deterministic Wilson-Cowan rate model, with the added bonus that the results quantitatively describe simulations of spiking EIF neurons. However, for modules with a moderate (~ 10⁴) number of neurons, stochastic variations in the spike emission of neurons are important and need to be taken into account. On the one hand, they modify the oscillations in a way that tends to promote synchronization between different modules. On the other hand, independent fluctuations on different modules tend to disrupt synchronization. The correlations between distant oscillatory modules can be described by stochastic equations for the oscillator phases that have been intensely studied in other contexts. On shorter distances, we develop a description that also takes into account amplitude modes and that quantitatively accounts for our simulation data. Stochastic dephasing of neighboring modules produces transient phase gradients and the transient appearance of phase waves. We propose that these stochastically-induced phase waves provide an explanative framework for the observations of traveling waves in the cortex during beta oscillations.

EEG microstate periodicity explained by rotating phase patterns of resting-state alpha oscillations

Spatio-temporal patterns in electroencephalography (EEG) can be described by microstate analysis, a discrete approximation of the continuous electric field patterns produced by the cerebral cortex. Resting-state EEG microstates are largely determined by alpha frequencies (8-12 Hz) and we recently demonstrated that microstates occur periodically with twice the alpha frequency. To understand the origin of microstate periodicity, we analyzed the analytic amplitude and the analytic phase of resting-state alpha oscillations independently. In continuous EEG data we found rotating phase patterns organized around a small number of phase singularities which varied in number and location. The spatial rotation of phase patterns occurred with the underlying alpha frequency. Phase rotors coincided with periodic microstate motifs involving the four canonical microstate maps. The analytic amplitude showed no oscillatory behaviour and was almost static across time intervals of 1-2 alpha cycles, resulting in the global pattern of a standing wave. In n=23 healthy adults, time-lagged mutual information analysis of microstate sequences derived from amplitude and phase signals of awake eyes-closed EEG records showed that only the phase component contributed to the periodicity of microstate sequences. Phase sequences showed mutual information peaks at multiples of 50 ms and the group average had a main peak at 100 ms (10 Hz), whereas amplitude sequences had a slow and monotonous information decay. This result was confirmed by an independent approach combining temporal principal component analysis (tPCA) and autocorrelation analysis. We reproduced our observations in a generic model of EEG oscillations composed of coupled non-linear oscillators (Stuart-Landau model). Phase-amplitude dynamics similar to experimental EEG occurred when the oscillators underwent a supercritical Hopf bifurcation, a common feature of many computational models of the alpha rhythm. These findings explain our previous description of periodic microstate recurrence and its relation to the time scale of alpha oscillations. Moreover, our results corroborate the predictions of computational models and connect experimentally observed EEG patterns to properties of critical oscillator networks.

Direction-selective motion discrimination by traveling waves in visual cortex

The majority of neurons in primary visual cortex respond selectively to bars of light that have a specific orientation and move in a specific direction. The spatial and temporal responses of such neurons are non-separable. How neurons accomplish that computational feat without resort to explicit time delays is unknown. We propose a novel neural mechanism whereby visual cortex computes non-separable responses by generating endogenous traveling waves of neural activity that resonate with the space-time signature of the visual stimulus. The spatiotemporal characteristics of the response are defined by the local topology of excitatory and inhibitory lateral connections in the cortex. We simulated the interaction between endogenous traveling waves and the visual stimulus using spatially distributed populations of excitatory and inhibitory neurons with Wilson-Cowan dynamics and inhibitory-surround coupling. Our model reliably detected visual gratings that moved with a given speed and direction provided that we incorporated neural competition to suppress false motion signals in the opposite direction. The findings suggest that endogenous traveling waves in visual cortex can impart direction-selectivity on neural responses without resort to explicit time delays. They also suggest a functional role for motion opponency in eliminating false motion signals.

It is well established that the so-called ‘simple cells’ of the primary visual cortex respond preferentially to oriented bars of light that move across the visual field with a particular speed and direction. The spatiotemporal responses of such neurons are said to be non-separable because they cannot be constructed from independent spatial and temporal neural mechanisms. Contemporary theories of how neurons compute non-separable responses typically rely on finely tuned transmission delays between signals from disparate regions of the visual field. However the existence of such delays is controversial. We propose an alternative neural mechanism for computing non-separable responses that does not require transmission delays. It instead relies on the predisposition of the cortical tissue to spontaneously generate spatiotemporal waves of neural activity that travel with a particular speed and direction. We propose that the endogenous wave activity resonates with the visual stimulus to elicit direction-selective neural responses to visual motion. We demonstrate the principle in computer models and show that competition between opposing neurons robustly enhances their ability to discriminate between visual gratings that move in opposite directions.

Nonlocal Raman response in Kerr resonators: Moving temporal localized structures and bifurcation structure

A ring resonator made of a silica-based optical fiber is a paradigmatic system for the generation of dissipative localized structures or dissipative solitons. We analyze the effect of the non-instantaneous nonlinear response of the fused silica or the Raman response on the formation of localized structures. After reducing the generalized Lugiato-Lefever to a simple and generic bistable model with a nonlocal Raman effect, we investigate analytically the formation of moving temporal localized structures. This reduction is valid close to the nascent bistability regime, where the system undergoes a second-order critical point marking the onset of a hysteresis loop. The interaction between fronts allows for the stabilization of temporal localized structures. Without the Raman effect, moving temporal localized structures do not exist, as shown in M. G. Clerc, S. Coulibaly, and M. Tlidi, Phys. Rev. Res. 2, 013024 (2020). The detailed derivation of the speed and the width associated with these structures is presented. We characterize numerically in detail the bifurcation structure and stability associated with the moving temporal localized states. The numerical results of the governing equations are in close agreement with analytical predictions.

Before and beyond the Wilson-Cowan equations

The Wilson-Cowan equations represent a landmark in the history of computational neuroscience. Along with the insights Wilson and Cowan offered for neuroscience, they crystallized an approach to modeling neural dynamics and brain function. Although their iconic equations are used in various guises today, the ideas that led to their formulation and the relationship to other approaches are not well known. Here, we give a little context to some of the biological and theoretical concepts that lead to the Wilson-Cowan equations and discuss how to extend beyond them.

Moving bumps in theta neuron networks

We consider large networks of theta neurons on a ring, synaptically coupled with an asymmetric kernel. Such networks support stable "bumps" of activity, which move along the ring if the coupling kernel is asymmetric. We investigate the effects of the kernel asymmetry on the existence, stability, and speed of these moving bumps using continuum equations formally describing infinite networks. Depending on the level of heterogeneity within the network, we find complex sequences of bifurcations as the amount of asymmetry is varied, in strong contrast to the behavior of a classical neural field model.

A model for focal seizure onset, propagation, evolution, and progression

We developed a neural network model that can account for major elements common to human focal seizures. These include the tonic-clonic transition, slow advance of clinical semiology and corresponding seizure territory expansion, widespread EEG synchronization, and slowing of the ictal rhythm as the seizure approaches termination. These were reproduced by incorporating usage-dependent exhaustion of inhibition in an adaptive neural network that receives global feedback inhibition in addition to local recurrent projections. Our model proposes mechanisms that may underline common EEG seizure onset patterns and status epilepticus, and postulates a role for synaptic plasticity in the emergence of epileptic foci. Complex patterns of seizure activity and bi-stable seizure end-points arise when stochastic noise is included. With the rapid advancement of clinical and experimental tools, we believe that this model can provide a roadmap and potentially an in silico testbed for future explorations of seizure mechanisms and clinical therapies.

Bioelectrical domain walls in homogeneous tissues

Electrical signaling in biology is typically associated with action potentials, transient spikes in membrane voltage that return to baseline. Hodgkin-Huxley and related conductance-based models of electrophysiology belong to a more general class of reaction-diffusion equations which could, in principle, support spontaneous emergence of patterns of membrane voltage which are stable in time but structured in space. Here we show theoretically and experimentally that homogeneous or nearly homogeneous tissues can undergo spontaneous spatial symmetry breaking through a purely electrophysiological mechanism, leading to formation of domains with different resting potentials separated by stable bioelectrical domain walls. Transitions from one resting potential to another can occur through long-range migration of these domain walls. We map bioelectrical domain wall motion using all-optical electrophysiology in an engineered cell line and in human induced pluripotent stem cell (iPSC)-derived myoblasts. Bioelectrical domain wall migration may occur during embryonic development and during physiological signaling processes in polarized tissues. These results demonstrate that nominally homogeneous tissues can undergo spontaneous bioelectrical symmetry breaking.

Bumps and oscillons in networks of spiking neurons

We study localized patterns in an exact mean-field description of a spatially extended network of quadratic integrate-and-fire neurons. We investigate conditions for the existence and stability of localized solutions, so-called bumps, and give an analytic estimate for the parameter range, where these solutions exist in parameter space, when one or more microscopic network parameters are varied. We develop Galerkin methods for the model equations, which enable numerical bifurcation analysis of stationary and time-periodic spatially extended solutions. We study the emergence of patterns composed of multiple bumps, which are arranged in a snake-and-ladder bifurcation structure if a homogeneous or heterogeneous synaptic kernel is suitably chosen. Furthermore, we examine time-periodic, spatially localized solutions (oscillons) in the presence of external forcing, and in autonomous, recurrently coupled excitatory and inhibitory networks. In both cases, we observe period-doubling cascades leading to chaotic oscillations.

Cortical stimulation in aphasia following ischemic stroke: toward model-guided electrical neuromodulation

The aim of this paper is to integrate different bodies of research including brain traveling waves, brain neuromodulation, neural field modeling and post-stroke language disorders in order to explore the opportunity of implementing model-guided, cortical neuromodulation for the treatment of post-stroke aphasia. Worldwide according to WHO, strokes are the second leading cause of death and the third leading cause of disability. In ischemic stroke, there is not enough blood supply to provide enough oxygen and nutrients to parts of the brain, while in hemorrhagic stroke, there is bleeding within the enclosed cranial cavity. The present paper focuses on ischemic stroke. We first review accumulating observations of traveling waves occurring spontaneously or triggered by external stimuli in healthy subjects as well as in patients with brain disorders. We examine the putative functions of these waves and focus on post-stroke aphasia observed when brain language networks become fragmented and/or partly silent, thus perturbing the progression of traveling waves across perilesional areas. Secondly, we focus on a simplified model based on the current literature in the field and describe cortical traveling wave dynamics and their modulation. This model uses a biophysically realistic integro-differential equation describing spatially distributed and synaptically coupled neural networks producing traveling wave solutions. The model is used to calculate wave parameters (speed, amplitude and/or frequency) and to guide the reconstruction of the perturbed wave. A stimulation term is included in the model to restore wave propagation to a reasonably good level. Thirdly, we examine various issues related to the implementation model-guided neuromodulation in the treatment of post-stroke aphasia given that closed-loop invasive brain stimulation studies have recently produced encouraging results. Finally, we suggest that modulating traveling waves by acting selectively and dynamically across space and time to facilitate wave propagation is a promising therapeutic strategy especially at a time when a new generation of closed-loop cortical stimulation systems is about to arrive on the market.

The hyperbolic model for edge and texture detection in the primary visual cortex

The modeling of neural fields in the visual cortex involves geometrical structures which describe in mathematical formalism the functional architecture of this cortical area. The case of contour detection and orientation tuning has been extensively studied and has become a paradigm for the mathematical analysis of image processing by the brain. Ten years ago an attempt was made to extend these models by replacing orientation (an angle) with a second-order tensor built from the gradient of the image intensity, and it was named the structure tensor. This assumption does not follow from biological observations (experimental evidence is still lacking) but from the idea that the effectiveness of texture processing with the structure tensor in computer vision may well be exploited by the brain itself. The drawback is that in this case the geometry is not Euclidean but hyperbolic instead, which complicates the analysis substantially. The purpose of this review is to present the methodology that was developed in a series of papers to investigate this quite unusual problem, specifically from the point of view of tuning and pattern formation. These methods, which rely on bifurcation theory with symmetry in the hyperbolic context, might be of interest for the modeling of other features such as color vision or other brain functions.

Spherical-harmonics mode decomposition of neural field equations

Large-scale neural networks can be described in the spatial continuous limit by neural field equations. For large-scale brain networks, the connectivity is typically translationally variant and imposes a large computational burden upon simulations. To reduce this burden, we take a semiquantitative approach and study the dynamics of neural fields described by a delayed integrodifferential equation. We decompose the connectivity into spatially variant and invariant contributions, which typically comprise the short- and long-range fiber systems, respectively. The neural fields are mapped on the two-dimensional spherical surface, which is choice consistent with routine mappings of cortical surfaces. Then, we perform mathematically a mode decomposition of the neural field equation into spherical harmonic basis functions. A spatial truncation of the leading orders at low wave number is consistent with the spatially coherent pattern formation of large-scale patterns observed in simulations and empirical brain imaging data and leads to a low-dimensional representation of the dynamics of the neural fields, bearing promise for an acceleration of the numerical simulations by orders of magnitude.

Investigating the role of gap junctions in seizure wave propagation

The effect of gap junctions as well as the biological mechanisms behind seizure wave propagation is not completely understood. In this work, we use a simple neural field model to study the possible influence of gap junctions specifically on cortical wave propagation that has been observed in vivo preceding seizure termination. We consider a voltage-based neural field model consisting of an excitatory and an inhibitory population as well as both chemical and gap junction-like synapses. We are able to approximate important properties of cortical wave propagation previously observed in vivo before seizure termination. This model adds support to existing evidence from models and clinical data suggesting a key role of gap junctions in seizure wave propagation. In particular, we found that in this model gap junction-like connectivity determines the propagation of one-bump or two-bump traveling wave solutions with features consistent with the clinical data. For sufficiently increased gap junction connectivity, wave solutions cease to exist. Moreover, gap junction connectivity needs to be sufficiently low or moderate to permit the existence of linearly stable solutions of interest.

Emergence of spatially periodic diffusive waves in small-world neuronal networks

It has been observed in experiment that the anatomical structure of neuronal networks in the brain possesses the feature of small-world networks. Yet how the small-world structure affects network dynamics remains to be fully clarified. Here we study the dynamics of a class of small-world networks consisting of pulse-coupled integrate-and-fire (I&F) neurons. Under stochastic Poisson drive, we find that the activity of the entire network resembles diffusive waves. To understand its underlying mechanism, we analyze the simplified regular-lattice network consisting of firing-rate-based neurons as an approximation to the original I&F small-world network. We demonstrate both analytically and numerically that, with strongly coupled connections, in the absence of noise, the activity of the firing-rate-based regular-lattice network spatially forms a static grating pattern that corresponds to the spatial distribution of the firing rate observed in the I&F small-world neuronal network. We further show that the spatial grating pattern with different phases comprise the continuous attractor of both the I&F small-world and firing-rate-based regular-lattice network dynamics. In the presence of input noise, the activity of both networks is perturbed along the continuous attractor, which gives rise to the diffusive waves. Our numerical simulations and theoretical analysis may potentially provide insights into the understanding of the generation of wave patterns observed in cortical networks.

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.

Computing by modulating spontaneous cortical activity patterns as a mechanism of active visual processing

Cortical populations produce complex spatiotemporal activity spontaneously without sensory inputs. However, the fundamental computational roles of such spontaneous activity remain unclear. Here, we propose a new neural computation mechanism for understanding how spontaneous activity is actively involved in cortical processing: Computing by Modulating Spontaneous Activity (CMSA). Using biophysically plausible circuit models, we demonstrate that spontaneous activity patterns with dynamical properties, as found in empirical observations, are modulated or redistributed by external stimuli to give rise to neural responses. We find that this CMSA mechanism of generating neural responses provides profound computational advantages, such as actively speeding up cortical processing. We further reveal that the CMSA mechanism provides a unifying explanation for many experimental findings at both the single-neuron and circuit levels, and that CMSA in response to natural stimuli such as face images is the underlying neurophysiological mechanism of perceptual "bubbles" as found in psychophysical studies.

Scaling of a Large-Scale Simulation of Synchronous Slow-Wave and Asynchronous Awake-Like Activity of a Cortical Model With Long-Range Interconnections

Cortical synapse organization supports a range of dynamic states on multiple spatial and temporal scales, from synchronous slow wave activity (SWA), characteristic of deep sleep or anesthesia, to fluctuating, asynchronous activity during wakefulness (AW). Such dynamic diversity poses a challenge for producing efficient large-scale simulations that embody realistic metaphors of short- and long-range synaptic connectivity. In fact, during SWA and AW different spatial extents of the cortical tissue are active in a given timespan and at different firing rates, which implies a wide variety of loads of local computation and communication. A balanced evaluation of simulation performance and robustness should therefore include tests of a variety of cortical dynamic states. Here, we demonstrate performance scaling of our proprietary Distributed and Plastic Spiking Neural Networks (DPSNN) simulation engine in both SWA and AW for bidimensional grids of neural populations, which reflects the modular organization of the cortex. We explored networks up to 192 × 192 modules, each composed of 1,250 integrate-and-fire neurons with spike-frequency adaptation, and exponentially decaying inter-modular synaptic connectivity with varying spatial decay constant. For the largest networks the total number of synapses was over 70 billion. The execution platform included up to 64 dual-socket nodes, each socket mounting 8 Intel Xeon Haswell processor cores @ 2.40 GHz clock rate. Network initialization time, memory usage, and execution time showed good scaling performances from 1 to 1,024 processes, implemented using the standard Message Passing Interface (MPI) protocol. We achieved simulation speeds of between 2.3 × 109 and 4.1 × 109 synaptic events per second for both cortical states in the explored range of inter-modular interconnections.

A gradient flow formulation for the stochastic Amari neural field model

We study stochastic Amari-type neural field equations, which are mean-field models for neural activity in the cortex. We prove that under certain assumptions on the coupling kernel, the neural field model can be viewed as a gradient flow in a nonlocal Hilbert space. This makes all gradient flow methods available for the analysis, which could previously not be used, as it was not known, whether a rigorous gradient flow formulation exists. We show that the equation is well-posed in the nonlocal Hilbert space in the sense that solutions starting in this space also remain in it for all times and space-time regularity results hold for the case of spatially correlated noise. Uniqueness of invariant measures, ergodic properties for the associated Feller semigroups, and several examples of kernels are also discussed.

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.

Dynamic Computation in Visual Thalamocortical Networks

Contemporary neurodynamical frameworks, such as coordination dynamics and winnerless competition, posit that the brain approximates symbolic computation by transitioning between metastable attractive states. This article integrates these accounts with electrophysiological data suggesting that coherent, nested oscillations facilitate information representation and transmission in thalamocortical networks. We review the relationship between criticality, metastability, and representational capacity, outline existing methods for detecting metastable oscillatory patterns in neural time series data, and evaluate plausible spatiotemporal coding schemes based on phase alignment. We then survey the circuitry and the mechanisms underlying the generation of coordinated alpha and gamma rhythms in the primate visual system, with particular emphasis on the pulvinar and its role in biasing visual attention and awareness. To conclude the review, we begin to integrate this perspective with longstanding theories of consciousness and cognition.

Control theory in biology and medicine : Introduction to the special issue

From September-December 2017, the Mathematical Biosciences Institute at Ohio State University hosted a series of workshops on control theory in biology and medicine, including workshops on control and modulation of neuronal and motor systems, control of cellular and molecular systems, control of disease / personalized medicine across heterogeneous populations, and sensorimotor control of animals and robots. This special issue presents tutorials and research articles by several of the participants in the MBI workshops.