Spontaneous Synchrony Breaking

Research on synchronization of coupled oscillators has helped explain how uniform behavior emerges in populations of non-uniform systems. But explaining how uniform populations engage in sustainable non-uniform synchronization may prove to be just as fascinating

Macroscopic coherent structures in a stochastic neural network: from interface dynamics to coarse-grained bifurcation analysis

We study coarse pattern formation in a cellular automaton modelling a spatially-extended stochastic neural network. The model, originally proposed by Gong and Robinson (Phys Rev E 85(5):055,101(R), 2012), is known to support stationary and travelling bumps of localised activity. We pose the model on a ring and study the existence and stability of these patterns in various limits using a combination of analytical and numerical techniques. In a purely deterministic version of the model, posed on a continuum, we construct bumps and travelling waves analytically using standard interface methods from neural field theory. In a stochastic version with Heaviside firing rate, we construct approximate analytical probability mass functions associated with bumps and travelling waves. In the full stochastic model posed on a discrete lattice, where a coarse analytic description is unavailable, we compute patterns and their linear stability using equation-free methods. The lifting procedure used in the coarse time-stepper is informed by the analysis in the deterministic and stochastic limits. In all settings, we identify the synaptic profile as a mesoscopic variable, and the width of the corresponding activity set as a macroscopic variable. Stationary and travelling bumps have similar meso- and macroscopic profiles, but different microscopic structure, hence we propose lifting operators which use microscopic motifs to disambiguate them. We provide numerical evidence that waves are supported by a combination of high synaptic gain and long refractory times, while meandering bumps are elicited by short refractory times

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 oscillatorsfals

Neural field models with threshold noise

The original neural field model of Wilson and Cowan is often interpreted as the averaged behaviour of a network of switch like neural elements with a distribution of switch thresholds, giving rise to the classic sigmoidal population firing-rate function so prevalent in large scale neuronal modelling. In this paper we explore the effects of such threshold noise without recourse to averaging and show that spatial correlations can have a strong effect on the behaviour of waves and patterns in continuum models. Moreover, for a prescribed spatial covariance function we explore the differences in behaviour that can emerge when the underlying stationary distribution is changed from Gaussian to non-Gaussian. For travelling front solutions, in a system with exponentially decaying spatial interactions, we make use of an interface approach to calculate the instantaneous wave speed analytically as a series expansion in the noise strength. From this we find that, for weak noise, the spatially averaged speed depends only on the choice of covariance function and not on the shape of the stationary distribution. For a system with a Mexican-hat spatial connectivity we further find that noise can induce localised bump solutions, and using an interface stability argument show that there can be multiple stable solution branches

PERIODIC SOLUTIONS FOR A PAIR OF DELAY-COUPLED ACTIVE THETA NEURONS

We consider a pair of identical theta neurons in the active regime, each coupled to the other via a delayed Dirac delta function. The network can support periodic solutions and we concentrate on solutions for which the neurons are half a period out of phase with one another, and also solutions for which the neurons are perfectly synchronous. The dynamics are analytically solvable, so we can derive explicit expressions for the existence and stability of both types of solutions. We find two branches of solutions, connected by symmetry-broken solutions which arise when the period of a solution as a function of delay is at a maximum or a minimum.fals

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.fals

Combining spatial and parametric working memory in a dynamic neural field model

We present a novel dynamic neural field model consisting of two coupled fields of Amari-type which supports the existence of localized activity patterns or “bumps” with a continuum of amplitudes. Bump solutions have been used in the past to model spatial working memory. We apply the model to explain input-specific persistent activity that increases monotonically with the time integral of the input (parametric working memory). In numerical simulations of a multi-item memory task, we show that the model robustly memorizes the strength and/or duration of inputs. Moreover, and important for adaptive behavior in dynamic environments, the memory strength can be changed at any time by new behaviorally relevant information. A direct comparison of model behaviors shows that the 2-field model does not suffer the problems of the classical Amari model when the inputs are presented sequentially as opposed to simultaneously

Periodic solutions for a pair of delay-coupled excitable theta neurons

We consider a pair of identical theta neurons in the excitable regime, each coupled to the other via a delayed Dirac delta function with the same delay. This simple network can support different periodic solutions, and we concentrate on two important types: those for which the neurons are perfectly synchronous, and those where the neurons are exactly half a period out of phase and fire alternatingly. Owing to the specific type of pulsatile feedback, we are able to determine these solutions and their stability analytically. More specifically, (infinitely many) branches of periodic solutions of either type are created at saddle-node bifurcations, and they gain stability at symmetry-breaking bifurcations when their period as a function of the delay is at its minimum. We also determine the respective branches of symmetry-broken periodic solutions and show that they are all unstable. We demonstrate by considering smoothed pulse-like coupling that the special case of the Dirac delta function can be seen as a sort of normal form: the basic structure of the different periodic solutions of the two theta neurons is preserved, but there may be additional changes of stability along the different branches.fals

Theta neuron subject to delayed feedback: a prototypical model for self-sustained pulsing

We consider a single theta neuron with delayed self-feedback in the form of a Dirac delta function in time. Because the dynamics of a theta neuron on its own can be solved explicitly—it is either excitable or shows self-pulsations—we are able to derive algebraic expressions for the existence and stability of the periodic solutions that arise in the presence of feedback. These periodic solutions are characterized by one or more equally spaced pulses per delay interval, and there is an increasing amount of multistability with increasing delay time. We present a complete description of where these self-sustained oscillations can be found in parameter space; in particular, we derive explicit expressions for the loci of their saddle-node bifurcations. We conclude that the theta neuron with delayed self-feedback emerges as a prototypical model: it provides an analytical basis for understanding pulsating dynamics observed in other excitable systems subject to delayed self-coupling.fals

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.fals