Bumps in Small-World Networks

Pulse Shape and Voltage-Dependent Synchronization in Spiking Neuron Networks

Pulse-coupled spiking neural networks are a powerful tool to gain mechanistic insights into how neurons self-organize to produce coherent collective behavior. These networks use simple spiking neuron models, such as the θ-neuron or the quadratic integrate-and-fire (QIF) neuron, that replicate the essential features of real neural dynamics. Interactions between neurons are modeled with infinitely narrow pulses, or spikes, rather than the more complex dynamics of real synapses. To make these networks biologically more plausible, it has been proposed that they must also account for the finite width of the pulses, which can have a significant impact on the network dynamics. However, the derivation and interpretation of these pulses are contradictory, and the impact of the pulse shape on the network dynamics is largely unexplored. Here, I take a comprehensive approach to pulse coupling in networks of QIF and θ-neurons. I argue that narrow pulses activate voltage-dependent synaptic conductances and show how to implement them in QIF neurons such that their effect can last through the phase after the spike. Using an exact low-dimensional description for networks of globally coupled spiking neurons, I prove for instantaneous interactions that collective oscillations emerge due to an effective coupling through the mean voltage. I analyze the impact of the pulse shape by means of a family of smooth pulse functions with arbitrary finite width and symmetric or asymmetric shapes. For symmetric pulses, the resulting voltage coupling is not very effective in synchronizing neurons, but pulses that are slightly skewed to the phase after the spike readily generate collective oscillations. The results unveil a voltage-dependent spike synchronization mechanism at the heart of emergent collective behavior, which is facilitated by pulses of finite width and complementary to traditional synaptic transmission in spiking neuron networks.

Patched patterns and emergence of chaotic interfaces in arrays of nonlocally coupled excitable systems

We disclose a new class of patterns, called patched patterns, in arrays of non-locally coupled excitable units with attractive and repulsive interactions. The self-organization process involves the formation of two types of patches, majority and minority ones, characterized by uniform average spiking frequencies. Patched patterns may be temporally periodic, quasiperiodic, or chaotic, whereby chaotic patterns may further develop interfaces comprised of units with average frequencies in between those of majority and minority patches. Using chaos and bifurcation theory, we demonstrate that chaos typically emerges via a torus breakup and identify the secondary bifurcation that gives rise to chaotic interfaces. It is shown that the maximal Lyapunov exponent of chaotic patched patterns does not decay, but rather converges to a finite value with system size. Patched patterns with a smaller wavenumber may exhibit diffusive motion of chaotic interfaces, similar to that of the incoherent part of chimeras.

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.

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.

Alternating activity patterns and a chimeralike state in a network of globally coupled excitable Morris-Lecar neurons

Spatiotemporal chaos collapses to either a rest state or a propagating pulse in a ring network of diffusively coupled, excitable Morris-Lecar neurons. Adding global varying synaptic coupling to the ring network reveals complex transient behavior. Spatiotemporal chaos collapses into a transient pulse that reinitiates spatiotemporal chaos to allow sequential pattern switching until a collapse to the rest state. A domain of irregular neuron activity coexists with a domain of inactive neurons forming a transient chimeralike state. Transient spatial localization of the chimeralike state is observed for stronger synapses.

Spike-burst chimera states in an adaptive exponential integrate-and-fire neuronal network

Chimera states are spatiotemporal patterns in which coherence and incoherence coexist. We observe the coexistence of synchronous (coherent) and desynchronous (incoherent) domains in a neuronal network. The network is composed of coupled adaptive exponential integrate-and-fire neurons that are connected by means of chemical synapses. In our neuronal network, the chimera states exhibit spatial structures both with spike and burst activities. Furthermore, those desynchronized domains not only have either spike or burst activity, but we show that the structures switch between spikes and bursts as the time evolves. Moreover, we verify the existence of multicluster chimera states.

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.

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.

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.