Derivation of a neural field model from a network of theta neurons

Chimera states in brain networks: Empirical neural vs. modular fractal connectivity

Complex spatiotemporal patterns, called chimera states, consist of coexisting coherent and incoherent domains and can be observed in networks of coupled oscillators. The interplay of synchrony and asynchrony in complex brain networks is an important aspect in studies of both the brain function and disease. We analyse the collective dynamics of FitzHugh-Nagumo neurons in complex networks motivated by its potential application to epileptology and epilepsy surgery. We compare two topologies: an empirical structural neural connectivity derived from diffusion-weighted magnetic resonance imaging and a mathematically constructed network with modular fractal connectivity. We analyse the properties of chimeras and partially synchronized states and obtain regions of their stability in the parameter planes. Furthermore, we qualitatively simulate the dynamics of epileptic seizures and study the influence of the removal of nodes on the network synchronizability, which can be useful for applications to epileptic surgery.

The Dynamics of Networks of Identical Theta Neurons

We consider finite and infinite all-to-all coupled networks of identical theta neurons. Two types of synaptic interactions are investigated: instantaneous and delayed (via first-order synaptic processing). Extensive use is made of the Watanabe/Strogatz (WS) ansatz for reducing the dimension of networks of identical sinusoidally-coupled oscillators. As well as the degeneracy associated with the constants of motion of the WS ansatz, we also find continuous families of solutions for instantaneously coupled neurons, resulting from the reversibility of the reduced model and the form of the synaptic input. We also investigate a number of similar related models. We conclude that the dynamics of networks of all-to-all coupled identical neurons can be surprisingly complicated.

Optimal design of tweezer control for chimera states

Chimera states are complex spatio-temporal patterns which consist of coexisting domains of spatially coherent and incoherent dynamics in systems of coupled oscillators. In small networks, chimera states usually exhibit short lifetimes and erratic drifting of the spatial position of the incoherent domain. A tweezer feedback control scheme can stabilize and fix the position of chimera states. We analyze the action of the tweezer control in small nonlocally coupled networks of Van der Pol and FitzHugh-Nagumo oscillators, and determine the ranges of optimal control parameters. We demonstrate that the tweezer control scheme allows for stabilization of chimera states with different shapes, and can be used as an instrument for controlling the coherent domains size, as well as the maximum average frequency difference of the oscillators.

Firing rate equations require a spike synchrony mechanism to correctly describe fast oscillations in inhibitory networks

Recurrently coupled networks of inhibitory neurons robustly generate oscillations in the gamma band. Nonetheless, the corresponding Wilson-Cowan type firing rate equation for such an inhibitory population does not generate such oscillations without an explicit time delay. We show that this discrepancy is due to a voltage-dependent spike-synchronization mechanism inherent in networks of spiking neurons which is not captured by standard firing rate equations. Here we investigate an exact low-dimensional description for a network of heterogeneous canonical Class 1 inhibitory neurons which includes the sub-threshold dynamics crucial for generating synchronous states. In the limit of slow synaptic kinetics the spike-synchrony mechanism is suppressed and the standard Wilson-Cowan equations are formally recovered as long as external inputs are also slow. However, even in this limit synchronous spiking can be elicited by inputs which fluctuate on a time-scale of the membrane time-constant of the neurons. Our meanfield equations therefore represent an extension of the standard Wilson-Cowan equations in which spike synchrony is also correctly described.

Population models describing the average activity of large neuronal ensembles are a powerful mathematical tool to investigate the principles underlying cooperative function of large neuronal systems. However, these models do not properly describe the phenomenon of spike synchrony in networks of neurons. In particular, they fail to capture the onset of synchronous oscillations in networks of inhibitory neurons. We show that this limitation is due to a voltage-dependent synchronization mechanism which is naturally present in spiking neuron models but not captured by traditional firing rate equations. Here we investigate a novel set of macroscopic equations which incorporate both firing rate and membrane potential dynamics, and that correctly generate fast inhibition-based synchronous oscillations. In the limit of slow-synaptic processing oscillations are suppressed, and the model reduces to an equation formally equivalent to the Wilson-Cowan model.

Synchronization of coupled active rotators by common noise

We study the effect of common noise on coupled active rotators. While such a noise always facilitates synchrony, coupling may be attractive (synchronizing) or repulsive (desynchronizing). We develop an analytical approach based on a transformation to approximate angle-action variables and averaging over fast rotations. For identical rotators, we describe a transition from full to partial synchrony at a critical value of repulsive coupling. For nonidentical rotators, the most nontrivial effect occurs at moderate repulsive coupling, where a juxtaposition of phase locking with frequency repulsion (anti-entrainment) is observed. We show that the frequency repulsion obeys a nontrivial power law.

Oscillators that sync and swarm

Synchronization occurs in many natural and technological systems, from cardiac pacemaker cells to coupled lasers. In the synchronized state, the individual cells or lasers coordinate the timing of their oscillations, but they do not move through space. A complementary form of self-organization occurs among swarming insects, flocking birds, or schooling fish; now the individuals move through space, but without conspicuously altering their internal states. Here we explore systems in which both synchronization and swarming occur together. Specifically, we consider oscillators whose phase dynamics and spatial dynamics are coupled. We call them swarmalators, to highlight their dual character. A case study of a generalized Kuramoto model predicts five collective states as possible long-term modes of organization. These states may be observable in groups of sperm, Japanese tree frogs, colloidal suspensions of magnetic particles, and other biological and physical systems in which self-assembly and synchronization interact.

Collective self-organized behavior can be observed in a variety of systems such as colloids and microswimmers. Here O’Keeffe et al. propose a model of oscillators which move in space and tend to synchronize with neighboring oscillators and outline five types of collective self-organized states.

Synchrony-induced modes of oscillation of a neural field model

We investigate the modes of oscillation of heterogeneous ring networks of quadratic integrate-and-fire (QIF) neurons with nonlocal, space-dependent coupling. Perturbations of the equilibrium state with a particular wave number produce transient standing waves with a specific temporal frequency, analogously to those in a tense string. In the neuronal network, the equilibrium corresponds to a spatially homogeneous, asynchronous state. Perturbations of this state excite the network's oscillatory modes, which reflect the interplay of episodes of synchronous spiking with the excitatory-inhibitory spatial interactions. In the thermodynamic limit, an exact low-dimensional neural field model describing the macroscopic dynamics of the network is derived. This allows us to obtain formulas for the Turing eigenvalues of the spatially homogeneous state and hence to obtain its stability boundary. We find that the frequency of each Turing mode depends on the corresponding Fourier coefficient of the synaptic pattern of connectivity. The decay rate instead is identical for all oscillation modes as a consequence of the heterogeneity-induced desynchronization of the neurons. Finally, we numerically compute the spectrum of spatially inhomogeneous solutions branching from the Turing bifurcation, showing that similar oscillatory modes operate in neural bump states and are maintained away from onset.

Synchronization scenarios in the Winfree model of coupled oscillators

Fifty years ago Arthur Winfree proposed a deeply influential mean-field model for the collective synchronization of large populations of phase oscillators. Here we provide a detailed analysis of the model for some special, analytically tractable cases. Adopting the thermodynamic limit, we derive an ordinary differential equation that exactly describes the temporal evolution of the macroscopic variables in the Ott-Antonsen invariant manifold. The low-dimensional model is then thoroughly investigated for a variety of pulse types and sinusoidal phase response curves (PRCs). Two structurally different synchronization scenarios are found, which are linked via the mutation of a Bogdanov-Takens point. From our results, we infer a general rule of thumb relating pulse shape and PRC offset with each scenario. Finally, we compare the exact synchronization threshold with the prediction of the averaging approximation given by the Kuramoto-Sakaguchi model. At the leading order, the discrepancy appears to behave as an odd function of the PRC offset.

An optimal strategy for epilepsy surgery: Disruption of the rich-club?

Surgery is a therapeutic option for people with epilepsy whose seizures are not controlled by anti-epilepsy drugs. In pre-surgical planning, an array of data modalities, often including intra-cranial EEG, is used in an attempt to map regions of the brain thought to be crucial for the generation of seizures. These regions are then resected with the hope that the individual is rendered seizure free as a consequence. However, post-operative seizure freedom is currently sub-optimal, suggesting that the pre-surgical assessment may be improved by taking advantage of a mechanistic understanding of seizure generation in large brain networks. Herein we use mathematical models to uncover the relative contribution of regions of the brain to seizure generation and consequently which brain regions should be considered for resection. A critical advantage of this modeling approach is that the effect of different surgical strategies can be predicted and quantitatively compared in advance of surgery. Herein we seek to understand seizure generation in networks with different topologies and study how the removal of different nodes in these networks reduces the occurrence of seizures. Since this a computationally demanding problem, a first step for this aim is to facilitate tractability of this approach for large networks. To do this, we demonstrate that predictions arising from a neural mass model are preserved in a lower dimensional, canonical model that is quicker to simulate. We then use this simpler model to study the emergence of seizures in artificial networks with different topologies, and calculate which nodes should be removed to render the network seizure free. We find that for scale-free and rich-club networks there exist specific nodes that are critical for seizure generation and should therefore be removed, whereas for small-world networks the strategy should instead focus on removing sufficient brain tissue. We demonstrate the validity of our approach by analysing intra-cranial EEG recordings from a database comprising 16 patients who have undergone epilepsy surgery, revealing rich-club structures within the obtained functional networks. We show that the postsurgical outcome for these patients was better when a greater proportion of the rich club was removed, in agreement with our theoretical predictions.

Uncovering low dimensional macroscopic chaotic dynamics of large finite size complex systems

In the last decade, it has been shown that a large class of phase oscillator models admit low dimensional descriptions for the macroscopic system dynamics in the limit of an infinite number N of oscillators. The question of whether the macroscopic dynamics of other similar systems also have a low dimensional description in the infinite N limit has, however, remained elusive. In this paper, we show how techniques originally designed to analyze noisy experimental chaotic time series can be used to identify effective low dimensional macroscopic descriptions from simulations with a finite number of elements. We illustrate and verify the effectiveness of our approach by applying it to the dynamics of an ensemble of globally coupled Landau-Stuart oscillators for which we demonstrate low dimensional macroscopic chaotic behavior with an effective 4-dimensional description. By using this description, we show that one can calculate dynamical invariants such as Lyapunov exponents and attractor dimensions. One could also use the reconstruction to generate short-term predictions of the macroscopic dynamics.

Frequency and phase synchronization in large groups: Low dimensional description of synchronized clapping, firefly flashing, and cricket chirping

A common observation is that large groups of oscillatory biological units often have the ability to synchronize. A paradigmatic model of such behavior is provided by the Kuramoto model, which achieves synchronization through coupling of the phase dynamics of individual oscillators, while each oscillator maintains a different constant inherent natural frequency. Here we consider the biologically likely possibility that the oscillatory units may be capable of enhancing their synchronization ability by adaptive frequency dynamics. We propose a simple augmentation of the Kuramoto model which does this. We also show that, by the use of a previously developed technique [Ott and Antonsen, Chaos 18, 037113 (2008)], it is possible to reduce the resulting dynamics to a lower dimensional system for the macroscopic evolution of the oscillator ensemble. By employing this reduction, we investigate the dynamics of our system, finding a characteristic hysteretic behavior and enhancement of the quality of the achieved synchronization.

Modeling the network dynamics of pulse-coupled neurons

We derive a mean-field approximation for the macroscopic dynamics of large networks of pulse-coupled theta neurons in order to study the effects of different network degree distributions and degree correlations (assortativity). Using the ansatz of Ott and Antonsen [Chaos 18, 037113 (2008)], we obtain a reduced system of ordinary differential equations describing the mean-field dynamics, with significantly lower dimensionality compared with the complete set of dynamical equations for the system. We find that, for sufficiently large networks and degrees, the dynamical behavior of the reduced system agrees well with that of the full network. This dimensional reduction allows for an efficient characterization of system phase transitions and attractors. For networks with tightly peaked degree distributions, the macroscopic behavior closely resembles that of fully connected networks previously studied by others. In contrast, networks with highly skewed degree distributions exhibit different macroscopic dynamics due to the emergence of degree dependent behavior of different oscillators. For nonassortative networks (i.e., networks without degree correlations), we observe the presence of a synchronously firing phase that can be suppressed by the presence of either assortativity or disassortativity in the network. We show that the results derived here can be used to analyze the effects of network topology on macroscopic behavior in neuronal networks in a computationally efficient fashion.

Ott-Antonsen attractiveness for parameter-dependent oscillatory systems

The Ott-Antonsen (OA) ansatz [Ott and Antonsen, Chaos 18, 037113 (2008); Chaos 19, 023117 (2009)] has been widely used to describe large systems of coupled phase oscillators. If the coupling is sinusoidal and if the phase dynamics does not depend on the specific oscillator, then the macroscopic behavior of the systems can be fully described by a low-dimensional dynamics. Does the corresponding manifold remain attractive when introducing an intrinsic dependence between an oscillator's phase and its dynamics by additional, oscillator specific parameters? To answer this, we extended the OA ansatz and proved that parameter-dependent oscillatory systems converge to the OA manifold given certain conditions. Our proof confirms recent numerical findings that already hinted at this convergence. Furthermore, we offer a thorough mathematical underpinning for networks of so-called theta neurons, where the OA ansatz has just been applied. In a final step, we extend our proof by allowing for time-dependent and multi-dimensional parameters as well as for network topologies other than global coupling. This renders the OA ansatz an excellent starting point for the analysis of a broad class of realistic settings.

Average activity of excitatory and inhibitory neural populations

We develop an extension of the Ott-Antonsen method [E. Ott and T. M. Antonsen, Chaos 18(3), 037113 (2008)] that allows obtaining the mean activity (spiking rate) of a population of excitable units. By means of the Ott-Antonsen method, equations for the dynamics of the order parameters of coupled excitatory and inhibitory populations of excitable units are obtained, and their mean activities are computed. Two different excitable systems are studied: Adler units and theta neurons. The resulting bifurcation diagrams are compared with those obtained from studying the phenomenological Wilson-Cowan model in some regions of the parameter space. Compatible behaviors, as well as higher dimensional chaotic solutions, are observed. We study numerical simulations to further validate the equations.

Travelling waves in arrays of delay-coupled phase oscillators

We consider the effects of several forms of delays on the existence and stability of travelling waves in non-locally coupled networks of Kuramoto-type phase oscillators and theta neurons. By passing to the continuum limit and using the Ott/Antonsen ansatz, we derive evolution equations for a spatially dependent order parameter. For phase oscillator networks, the travelling waves take the form of uniformly twisted waves, and these can often be characterised analytically. For networks of theta neurons, the waves are studied numerically.

Macroscopic self-oscillations and aging transition in a network of synaptically coupled quadratic integrate-and-fire neurons

We analyze the dynamics of a large network of coupled quadratic integrate-and-fire neurons, which represent the canonical model for class I neurons near the spiking threshold. The network is heterogeneous in that it includes both inherently spiking and excitable neurons. The coupling is global via synapses that take into account the finite width of synaptic pulses. Using a recently developed reduction method based on the Lorentzian ansatz, we derive a closed system of equations for the neuron's firing rate and the mean membrane potential, which are exact in the infinite-size limit. The bifurcation analysis of the reduced equations reveals a rich scenario of asymptotic behavior, the most interesting of which is the macroscopic limit-cycle oscillations. It is shown that the finite width of synaptic pulses is a necessary condition for the existence of such oscillations. The robustness of the oscillations against aging damage, which transforms spiking neurons into nonspiking neurons, is analyzed. The validity of the reduced equations is confirmed by comparing their solutions with the solutions of microscopic equations for the finite-size networks.

Bumps in Small-World Networks

We consider a network of coupled excitatory and inhibitory theta neurons which is capable of supporting stable spatially-localized “bump” solutions. We randomly add long-range and simultaneously remove short-range connections within the network to form a small-world network and investigate the effects of this rewiring on the existence and stability of the bump solution. We consider two limits in which continuum equations can be derived; bump solutions are fixed points of these equations. We can thus use standard numerical bifurcation analysis to determine the stability of these bumps and to follow them as parameters (such as rewiring probabilities) are varied. We find that under some rewiring schemes bumps are quite robust, whereas in other schemes they can become unstable via Hopf bifurcation or even be destroyed in saddle-node bifurcations.

Heterogeneity of time delays determines synchronization of coupled oscillators

Network couplings of oscillatory large-scale systems, such as the brain, have a space-time structure composed of connection strengths and signal transmission delays. We provide a theoretical framework, which allows treating the spatial distribution of time delays with regard to synchronization, by decomposing it into patterns and therefore reducing the stability analysis into the tractable problem of a finite set of delay-coupled differential equations. We analyze delay-structured networks of phase oscillators and we find that, depending on the heterogeneity of the delays, the oscillators group in phase-shifted, anti-phase, steady, and non-stationary clusters, and analytically compute their stability boundaries. These results find direct application in the study of brain oscillations.

Dynamics of a population of oscillatory and excitable elements

We analyze a variant of a model proposed by Kuramoto, Shinomoto, and Sakaguchi for a large population of coupled oscillatory and excitable elements. Using the Ott-Antonsen ansatz, we reduce the behavior of the population to a two-dimensional dynamical system with three parameters. We present the stability diagram and calculate several of its bifurcation curves analytically, for both excitatory and inhibitory coupling. Our main result is that when the coupling function is broad, the system can display bistability between steady states of constant high and low activity, whereas when the coupling function is narrow and inhibitory, one of the states in the bistable regime can show persistent pulsations in activity.

Tweezers for Chimeras in Small Networks

We propose a control scheme which can stabilize and fix the position of chimera states in small networks. Chimeras consist of coexisting domains of spatially coherent and incoherent dynamics in systems of nonlocally coupled identical oscillators. Chimera states are generally difficult to observe in small networks due to their short lifetime and erratic drifting of the spatial position of the incoherent domain. The control scheme, like a tweezer, might be useful in experiments, where usually only small networks can be realized.

Regular and irregular patterns of self-localized excitation in arrays of coupled phase oscillators

We study a system of phase oscillators with nonlocal coupling in a ring that supports self-organized patterns of coherence and incoherence, called chimera states. Introducing a global feedback loop, connecting the phase lag to the order parameter, we can observe chimera states also for systems with a small number of oscillators. Numerical simulations show a huge variety of regular and irregular patterns composed of localized phase slipping events of single oscillators. Using methods of classical finite dimensional chaos and bifurcation theory, we can identify the emergence of chaotic chimera states as a result of transitions to chaos via period doubling cascades, torus breakup, and intermittency. We can explain the observed phenomena by a mechanism of self-modulated excitability in a discrete excitable medium.

Macroscopic complexity from an autonomous network of networks of theta neurons

We examine the emergence of collective dynamical structures and complexity in a network of interacting populations of neuronal oscillators. Each population consists of a heterogeneous collection of globally-coupled theta neurons, which are a canonical representation of Type-1 neurons. For simplicity, the populations are arranged in a fully autonomous driver-response configuration, and we obtain a full description of the asymptotic macroscopic dynamics of this network. We find that the collective macroscopic behavior of the response population can exhibit equilibrium and limit cycle states, multistability, quasiperiodicity, and chaos, and we obtain detailed bifurcation diagrams that clarify the transitions between these macrostates. Furthermore, we show that despite the complexity that emerges, it is possible to understand the complicated dynamical structure of this system by building on the understanding of the collective behavior of a single population of theta neurons. This work is a first step in the construction of a mathematically-tractable network-of-networks representation of neuronal network dynamics.