Periodic solutions in next generation neural field models.

We consider a next generation neural field model which describes the dynamics of a network of theta neurons on a ring. For some parameters the network supports stable time-periodic solutions. Using the fact that the dynamics at each spatial location are described by a complex-valued Riccati equation we derive a self-consistency equation that such periodic solutions must satisfy. We determine the stability of these solutions, and present numerical results to illustrate the usefulness of this technique. The generality of this approach is demonstrated through its application to several other systems involving delays, two-population architecture and networks of Winfree oscillators.

Chimeras on a ring of oscillator populations.

Chimeras occur in networks of coupled oscillators and are characterized by coexisting groups of synchronous oscillators and asynchronous oscillators. We consider a network formed from N equal-sized populations at equally spaced points around a ring. We use the Ott/Antonsen ansatz to derive coupled ordinary differential equations governing the level of synchrony within each population and describe chimeras using a self-consistency argument. For N=2 and 3, our results are compared with previously known ones. We obtain new results for the cases of 4,5,…,12 populations and a numerically based conjecture resulting from the behavior of larger numbers of populations. We find macroscopic chaos when more than five populations are considered, but conjecture that this behavior vanishes as the number of populations is increased.

Chimeras in phase oscillator networks locally coupled through an auxiliary field: Stability and bifurcations.

We study networks in the form of a lattice of nodes with a large number of phase oscillators and an auxiliary variable at each node. The only interactions between nodes are nearest-neighbor. The Ott/Antonsen ansatz is used to derive equations for the order parameters of the phase oscillators at each node, resulting in a set of coupled ordinary differential equations. Chimeras are steady states of these equations, and we follow them as parameters are varied, determining their stability and bifurcations. In two-dimensional domains, we find that spiral wave chimeras and rotating waves have significantly different properties than those in networks with nonlocal coupling.

Chimeras on annuli.

Chimeras occur in networks of coupled oscillators and are characterized by the coexistence of synchronous and asynchronous groups of oscillators in different parts of the network. We consider a network of nonlocally coupled phase oscillators on an annular domain. The Ott/Antonsen ansatz is used to derive a continuum level description of the oscillators' expected dynamics in terms of a complex-valued order parameter. The equations for this order parameter are numerically analyzed in order to investigate solutions with the same symmetry as the domain and chimeras which are analogous to the "multi-headed" chimeras observed on one-dimensional domains. Such solutions are stable only for domains with widths that are neither too large nor too small. We also study rotating waves with different winding numbers, which are similar to spiral wave chimeras seen in two-dimensional domains. We determine ranges of parameters, such as the size of the domain for which such solutions exist and are stable, and the bifurcations by which they lose stability. All of these bifurcations appear subcritical.

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.

Global and local reduced models for interacting, heterogeneous agents.

Large collections of coupled, heterogeneous agents can manifest complex dynamical behavior presenting difficulties for simulation and analysis. However, if the collective dynamics lie on a low-dimensional manifold, then the original agent-based model may be approximated with a simplified surrogate model on and near the low-dimensional space where the dynamics live. Analytically identifying such simplified models can be challenging or impossible, but here we present a data-driven coarse-graining methodology for discovering such reduced models. We consider two types of reduced models: globally based models that use global information and predict dynamics using information from the whole ensemble and locally based models that use local information, that is, information from just a subset of agents close (close in heterogeneity space, not physical space) to an agent, to predict the dynamics of an agent. For both approaches, we are able to learn laws governing the behavior of the reduced system on the low-dimensional manifold directly from time series of states from the agent-based system. These laws take the form of either a system of ordinary differential equations (ODEs), for the globally based approach, or a partial differential equation (PDE) in the locally based case. For each technique, we employ a specialized artificial neural network integrator that has been templated on an Euler time stepper (i.e., a ResNet) to learn the laws of the reduced model. As part of our methodology, we utilize the proper orthogonal decomposition (POD) to identify the low-dimensional space of the dynamics. Our globally based technique uses the resulting POD basis to define a set of coordinates for the agent states in this space and then seeks to learn the time evolution of these coordinates as a system of ODEs. For the locally based technique, we propose a methodology for learning a partial differential equation representation of the agents; the PDE law depends on the state variables and partial derivatives of the state variables with respect to model heterogeneities. We require that the state variables are smooth with respect to model heterogeneities, which permit us to cast the discrete agent-based problem as a continuous one in heterogeneity space. The agents in such a representation bear similarity to the discretization points used in typical finite element/volume methods. As an illustration of the efficacy of our techniques, we consider a simplified coupled neuron model for rhythmic oscillations in the pre-Bötzinger complex and demonstrate how our data-driven surrogate models are able to produce dynamics comparable to the dynamics of the full system. A nontrivial conclusion is that the dynamics can be equally well reproduced by an all-to-all coupled and by a locally coupled model of the same agents.

The effects of within-neuron degree correlations in networks of spiking neurons.

We consider the effects of correlations between the in- and out-degrees of individual neurons on the dynamics of a network of neurons. By using theta neurons, we can derive a set of coupled differential equations for the expected dynamics of neurons with the same in-degree. A Gaussian copula is used to introduce correlations between a neuron's in- and out-degree, and numerical bifurcation analysis is used determine the effects of these correlations on the network's dynamics. For excitatory coupling, we find that inducing positive correlations has a similar effect to increasing the coupling strength between neurons, while for inhibitory coupling it has the opposite effect. We also determine the propensity of various two- and three-neuron motifs to occur as correlations are varied and give a plausible explanation for the observed changes in dynamics.

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.

Chaos in small networks of theta neurons.

We consider small networks of instantaneously coupled theta neurons. For inhibitory coupling and fixed parameter values, some initial conditions give chaotic solutions while others give quasiperiodic ones. This behaviour seems to result from the reversibility of the equations governing the networks' dynamics. We investigate the robustness of the chaotic behaviour with respect to changes in initial conditions and parameters and find the behaviour to be quite robust as long as the reversibility of the system is preserved.

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.

Equation-free analysis of spike-timing-dependent plasticity.

Spike-timing-dependent plasticity is the process by which the strengths of connections between neurons are modified as a result of the precise timing of the action potentials fired by the neurons. We consider a model consisting of one integrate-and-fire neuron receiving excitatory inputs from a large number-here, 1000-of Poisson neurons whose synapses are plastic. When correlations are introduced between the firing times of these input neurons, the distribution of synaptic strengths shows interesting, and apparently low-dimensional, dynamical behaviour. This behaviour is analysed in two different parameter regimes using equation-free techniques, which bypass the explicit derivation of the relevant low-dimensional dynamical system. We demonstrate both coarse projective integration (which speeds up the time integration of a dynamical system) and the use of recently developed data mining techniques to identify the appropriate low-dimensional description of the complex dynamical systems in our model.

Partially coherent twisted states in arrays of coupled phase oscillators.

We consider a one-dimensional array of phase oscillators with non-local coupling and a Lorentzian distribution of natural frequencies. The primary objects of interest are partially coherent states that are uniformly "twisted" in space. To analyze these, we take the continuum limit, perform an Ott/Antonsen reduction, integrate over the natural frequencies, and study the resulting spatio-temporal system on an unbounded domain. We show that these twisted states and their stability can be calculated explicitly. We find that stable twisted states with different wave numbers appear for increasing coupling strength in the well-known Eckhaus scenario. Simulations of finite arrays of oscillators show good agreement with results of the analysis of the infinite system.

Disorder-induced dynamics in a pair of coupled heterogeneous phase oscillator networks.

We consider a pair of coupled heterogeneous phase oscillator networks and investigate their dynamics in the continuum limit as the intrinsic frequencies of the oscillators are made more and more disparate. The Ott/Antonsen Ansatz is used to reduce the system to three ordinary differential equations. We find that most of the interesting dynamics, such as chaotic behaviour, can be understood by analysing a gluing bifurcation of periodic orbits; these orbits can be thought of as "breathing chimeras" in the limit of identical oscillators. We also add Gaussian white noise to the oscillators' dynamics and derive a pair of coupled Fokker-Planck equations describing the dynamics in this case. Comparison with simulations of finite networks of oscillators is used to confirm many of the results.

Chimeras in random non-complete networks of phase oscillators.

We consider the simplest network of coupled non-identical phase oscillators capable of displaying a "chimera" state (namely, two subnetworks with strong coupling within the subnetworks and weaker coupling between them) and systematically investigate the effects of gradually removing connections within the network, in a random but systematically specified way. We average over ensembles of networks with the same random connectivity but different intrinsic oscillator frequencies and derive ordinary differential equations (ODEs), whose fixed points describe a typical chimera state in a representative network of phase oscillators. Following these fixed points as parameters are varied we find that chimera states are quite sensitive to such random removals of connections, and that oscillations of chimera states can be either created or suppressed in apparent bifurcation points, depending on exactly how the connections are gradually removed.

Chimeras with uniformly distributed heterogeneity: Two coupled populations.

Chimeras occur in networks of two coupled populations of oscillators when the oscillators in one population synchronize while those in the other are asynchronous. We consider chimeras of this form in networks of planar oscillators for which one parameter associated with the dynamics of an oscillator is randomly chosen from a uniform distribution. A generalization of the previous approach [Laing, Phys. Rev. E 100, 042211 (2019)2470-004510.1103/PhysRevE.100.042211], which dealt with identical oscillators, is used to investigate the existence and stability of chimeras for these heterogeneous networks in the limit of an infinite number of oscillators. In all cases, making the oscillators more heterogeneous destroys the stable chimera in a saddle-node bifurcation. The results help us understand the robustness of chimeras in networks of general oscillators to heterogeneity.

Collective states in a ring network of theta neurons.

We consider a ring network of theta neurons with non-local homogeneous coupling. We analyse the corresponding continuum evolution equation, analytically describing all possible steady states and their stability. By considering a number of different parameter sets, we determine the typical bifurcation scenarios of the network, and put on a rigorous footing some previously observed numerical results.

Learning emergent partial differential equations in a learned emergent space.

We propose an approach to learn effective evolution equations for large systems of interacting agents. This is demonstrated on two examples, a well-studied system of coupled normal form oscillators and a biologically motivated example of coupled Hodgkin-Huxley-like neurons. For such types of systems there is no obvious space coordinate in which to learn effective evolution laws in the form of partial differential equations. In our approach, we accomplish this by learning embedding coordinates from the time series data of the system using manifold learning as a first step. In these emergent coordinates, we then show how one can learn effective partial differential equations, using neural networks, that do not only reproduce the dynamics of the oscillator ensemble, but also capture the collective bifurcations when system parameters vary. The proposed approach thus integrates the automatic, data-driven extraction of emergent space coordinates parametrizing the agent dynamics, with machine-learning assisted identification of an emergent PDE description of the dynamics in this parametrization.

Effects of degree distributions in random networks of type-I neurons.

We consider large networks of theta neurons and use the Ott-Antonsen ansatz to derive degree-based mean-field equations governing the expected dynamics of the networks. Assuming random connectivity, we investigate the effects of varying the widths of the in- and out-degree distributions on the dynamics of excitatory or inhibitory synaptically coupled networks and gap junction coupled networks. For synaptically coupled networks, the dynamics are independent of the out-degree distribution. Broadening the in-degree distribution destroys oscillations in inhibitory networks and decreases the range of bistability in excitatory networks. For gap junction coupled neurons, broadening the degree distribution varies the values of parameters at which there is an onset of collective oscillations. Many of the results are shown to also occur in networks of more realistic neurons.

Dynamics of Structured Networks of Winfree Oscillators.

Winfree oscillators are phase oscillator models of neurons, characterized by their phase response curve and pulsatile interaction function. We use the Ott/Antonsen ansatz to study large heterogeneous networks of Winfree oscillators, deriving low-dimensional differential equations which describe the evolution of the expected state of networks of oscillators. We consider the effects of correlations between an oscillator's in-degree and out-degree, and between the in- and out-degrees of an "upstream" and a "downstream" oscillator (degree assortativity). We also consider correlated heterogeneity, where some property of an oscillator is correlated with a structural property such as degree. We finally consider networks with parameter assortativity, coupling oscillators according to their intrinsic frequencies. The results show how different types of network structure influence its overall dynamics.

Solvable model of spiral wave chimeras.

Spiral waves are ubiquitous in two-dimensional systems of chemical or biological oscillators coupled locally by diffusion. At the center of such spirals is a phase singularity, a topological defect where the oscillator amplitude drops to zero. But if the coupling is nonlocal, a new kind of spiral can occur, with a circular core consisting of desynchronized oscillators running at full amplitude. Here, we provide the first analytical description of such a spiral wave chimera and use perturbation theory to calculate its rotation speed and the size of its incoherent core.