Adaptation controls synchrony and cluster states of coupled threshold-model neurons

Rapid changes in synchronizability in conductance-based neuronal networks with conductance-based coupling

Real neurons connect to each other non-randomly. These connectivity graphs can potentially impact the ability of networks to synchronize, along with the dynamics of neurons and the dynamics of their connections. How the connectivity of networks of conductance-based neuron models like the classical Hodgkin-Huxley model or the Morris-Lecar model impacts synchronizability remains unknown. One powerful tool to resolve the synchronizability of these networks is the master stability function (MSF). Here, we apply and extend the MSF approach to networks of Morris-Lecar neurons with conductance-based coupling to determine under which parameters and for which graphs the synchronous solutions are stable. We consider connectivity graphs with a constant non-zero row sum, where the MSF approach can be readily extended to conductance-based synapses rather than the more well-studied diffusive connectivity case, which primarily applies to gap junction connectivity. In this formulation, the synchronous solution is a single, self-coupled, or "autaptic" neuron. We find that the primary determining parameter for the stability of the synchronous solution is, unsurprisingly, the reversal potential, as it largely dictates the excitatory/inhibitory potential of a synaptic connection. However, the change between "excitatory" and "inhibitory" synapses is rapid, with only a few millivolts separating stability and instability of the synchronous state for most graphs. We also find that for specific coupling strengths (as measured by the global synaptic conductance), islands of synchronizability in the MSF can emerge for inhibitory connectivity. We verified the stability of these islands by direct simulation of pairs of neurons coupled with eigenvalues in the matching spectrum.

Blinking coupling enhances network synchronization

This paper studies the synchronization of a network with linear diffusive coupling, which blinks between the variables periodically. The synchronization of the blinking network in the case of sufficiently fast blinking is analyzed by showing that the stability of the synchronous solution depends only on the averaged coupling and not on the instantaneous coupling. To illustrate the effect of the blinking period on the network synchronization, the Hindmarsh-Rose model is used as the dynamics of nodes. The synchronization is investigated by considering constant single-variable coupling, averaged coupling, and blinking coupling through a linear stability analysis. It is observed that by decreasing the blinking period, the required coupling strength for synchrony is reduced. It equals that of the averaged coupling model times the number of variables. However, in the averaged coupling, all variables participate in the coupling, while in the blinking model only one variable is coupled at any time. Therefore, the blinking coupling leads to an enhanced synchronization in comparison with the single-variable coupling. Numerical simulations of the average synchronization error of the network confirm the results obtained from the linear stability analysis.

Desynchronization Transitions in Adaptive Networks

Adaptive networks change their connectivity with time, depending on their dynamical state. While synchronization in structurally static networks has been studied extensively, this problem is much more challenging for adaptive networks. In this Letter, we develop the master stability approach for a large class of adaptive networks. This approach allows for reducing the synchronization problem for adaptive networks to a low-dimensional system, by decoupling topological and dynamical properties. We show how the interplay between adaptivity and network structure gives rise to the formation of stability islands. Moreover, we report a desynchronization transition and the emergence of complex partial synchronization patterns induced by an increasing overall coupling strength. We illustrate our findings using adaptive networks of coupled phase oscillators and FitzHugh-Nagumo neurons with synaptic plasticity.

Influence of Autapses on Synchronization in Neural Networks With Chemical Synapses

A great deal of research has been devoted on the investigation of neural dynamics in various network topologies. However, only a few studies have focused on the influence of autapses, synapses from a neuron onto itself via closed loops, on neural synchronization. Here, we build a random network with adaptive exponential integrate-and-fire neurons coupled with chemical synapses, equipped with autapses, to study the effect of the latter on synchronous behavior. We consider time delay in the conductance of the pre-synaptic neuron for excitatory and inhibitory connections. Interestingly, in neural networks consisting of both excitatory and inhibitory neurons, we uncover that synchronous behavior depends on their synapse type. Our results provide evidence on the synchronous and desynchronous activities that emerge in random neural networks with chemical, inhibitory and excitatory synapses where neurons are equipped with autapses.

Discontinuity-induced intermittent synchronization transitions in coupled non-smooth systems

The synchronization transition in coupled non-smooth systems is studied for increasing coupling strength. The average order parameter is calculated to diagnose synchronization of coupled non-smooth systems. It is found that the coupled non-smooth system exhibits an intermittent synchronization transition from the cluster synchronization state to the complete synchronization state, depending on the coupling strength and initial conditions. Detailed numerical analyses reveal that the discontinuity always plays an important role in the synchronization transition of the coupled non-smooth system. In addition, it is found that increasing the coupling strength leads to the coexistence of periodic cluster states. Detailed research illustrates that the periodic clusters consist of two or more coexisting periodic attractors. Their periodic trajectory passes from one region to another region that is divided by discontinuous boundaries in the phase space. The mutual interactions of the local nonlinearity and the spatial coupling ultimately result in a stable periodic trajectory.

Synchronization of Electrically Coupled Resonate-and-Fire Neurons

Electrical coupling between neurons is broadly present across brain areas and is typically assumed to synchronize network activity. However, intrinsic properties of the coupled cells can complicate this simple picture. Many cell types with electrical coupling show a diversity of post-spike subthreshold fluctuations, often linked to subthreshold resonance, which are transmitted through electrical synapses in addition to action potentials. Using the theory of weakly coupled oscillators, we explore the effect of both subthreshold and spike-mediated coupling on synchrony in small networks of electrically coupled resonate-and-fire neurons, a hybrid neuron model with damped subthreshold oscillations and a range of post-spike voltage dynamics. We calculate the phase response curve using an extension of the adjoint method that accounts for the discontinuous post-spike reset rule. We find that both spikes and subthreshold fluctuations can jointly promote synchronization. The subthreshold contribution is strongest when the voltage exhibits a significant post-spike elevation in voltage, or plateau potential. Additionally, we show that the geometry of trajectories approaching the spiking threshold causes a "reset-induced shear" effect that can oppose synchrony in the presence of network asymmetry, despite having no effect on the phase-locking of symmetrically coupled pairs.

Clusters in nonsmooth oscillator networks

For coupled oscillator networks with Laplacian coupling, the master stability function (MSF) has proven a particularly powerful tool for assessing the stability of the synchronous state. Using tools from group theory, this approach has recently been extended to treat more general cluster states. However, the MSF and its generalizations require the determination of a set of Floquet multipliers from variational equations obtained by linearization around a periodic orbit. Since closed form solutions for periodic orbits are invariably hard to come by, the framework is often explored using numerical techniques. Here, we show that further insight into network dynamics can be obtained by focusing on piecewise linear (PWL) oscillator models. Not only do these allow for the explicit construction of periodic orbits, their variational analysis can also be explicitly performed. The price for adopting such nonsmooth systems is that many of the notions from smooth dynamical systems, and in particular linear stability, need to be modified to take into account possible jumps in the components of Jacobians. This is naturally accommodated with the use of saltation matrices. By augmenting the variational approach for studying smooth dynamical systems with such matrices we show that, for a wide variety of networks that have been used as models of biological systems, cluster states can be explicitly investigated. By way of illustration, we analyze an integrate-and-fire network model with event-driven synaptic coupling as well as a diffusively coupled network built from planar PWL nodes, including a reduction of the popular Morris-Lecar neuron model. We use these examples to emphasize that the stability of network cluster states can depend as much on the choice of single node dynamics as it does on the form of network structural connectivity. Importantly, the procedure that we present here, for understanding cluster synchronization in networks, is valid for a wide variety of systems in biology, physics, and engineering that can be described by PWL oscillators.

Robustness of chimera states in nonlocally coupled networks of nonidentical logistic maps

We investigate the dynamics of nonlocally coupled time-discrete maps with emphasis on the occurrence and robustness of chimera states. These peculiar, hybrid states are characterized by a coexistence of coherent and incoherent regions. We consider logistic maps coupled on a one-dimensional ring with finite coupling radius. Domains of chimera existence form different tongues in the parameter space of coupling range and coupling strength. For a sufficiently large coupling strength, each tongue refers to a wave number describing the structure of the spatial profile. We also analyze the period-adding scheme within these tongues and multiplicity of period solutions. Furthermore, we study the robustness of chimeras with respect to parameter inhomogeneities and find that these states persist for different widths of the parameter distribution. Finally, we explore the spatial structure of the chimera using a spatial correlation function.

How adaptation currents change threshold, gain, and variability of neuronal spiking

Many types of neurons exhibit spike rate adaptation, mediated by intrinsic slow K(+) currents, which effectively inhibit neuronal responses. How these adaptation currents change the relationship between in vivo like fluctuating synaptic input, spike rate output, and the spike train statistics, however, is not well understood. In this computational study we show that an adaptation current that primarily depends on the subthreshold membrane voltage changes the neuronal input-output relationship (I-O curve) subtractively, thereby increasing the response threshold, and decreases its slope (response gain) for low spike rates. A spike-dependent adaptation current alters the I-O curve divisively, thus reducing the response gain. Both types of an adaptation current naturally increase the mean interspike interval (ISI), but they can affect ISI variability in opposite ways. A subthreshold current always causes an increase of variability while a spike-triggered current decreases high variability caused by fluctuation-dominated inputs and increases low variability when the average input is large. The effects on I-O curves match those caused by synaptic inhibition in networks with asynchronous irregular activity, for which we find subtractive and divisive changes caused by external and recurrent inhibition, respectively. Synaptic inhibition, however, always increases the ISI variability. We analytically derive expressions for the I-O curve and ISI variability, which demonstrate the robustness of our results. Furthermore, we show how the biophysical parameters of slow K(+) conductances contribute to the two different types of an adaptation current and find that Ca(2+)-activated K(+) currents are effectively captured by a simple spike-dependent description, while muscarine-sensitive or Na(+)-activated K(+) currents show a dominant subthreshold component.