Synchronization in phase-coupled Kuramoto oscillator networks with axonal delay and synaptic plasticity

Hebbian and anti-Hebbian adaptation-induced dynamical states in adaptive networks

We investigate the interplay of an external forcing and an adaptive network, whose connection weights coevolve with the dynamical states of the phase oscillators. In particular, we consider the Hebbian and anti-Hebbian adaptation mechanisms for the evolution of the connection weights. The Hebbian adaptation manifests several interesting partially synchronized states, such as phase and frequency clusters, bump state, bump frequency phase clusters, and forced entrained clusters, in addition to the completely synchronized and forced entrained states. Anti-Hebbian adaptation facilitates the manifestation of the itinerant chimera characterized by randomly evolving coherent and incoherent domains along with some of the aforementioned dynamical states induced by the Hebbian adaptation. We introduce three distinct measures for the strength of incoherence based on the local standard deviations of the time-averaged frequency and the instantaneous phase of each oscillator, and the time-averaged mean frequency for each bin to corroborate the distinct dynamical states and to demarcate the two parameter phase diagrams. We also arrive at the existence and stability conditions for the forced entrained state using the linear stability analysis, which is found to be consistent with the simulation results.

Self-adaptation of networks of nonidentical pulse-coupled excitatory and inhibitory oscillators in the presence of distance-related delays to achieve frequency synchronization

We show that a network of nonidentical nodes, with excitable dynamics, pulse-coupled, with coupling delays depending on the Euclidean distance between nodes, is able to adapt the topology of its connections to obtain spike frequency synchronization. The adapted network exhibits remarkable properties: sparse, anticluster, necessary presence of a minimum of inhibitory nodes, predominance of connections from inhibitory nodes over those from excitatory nodes, and finally spontaneous spatial structuring of the inhibitory projections: the furthest are the most intense. In a second step, we discuss the possible implications of our findings to neural systems.

Chimera Patterns of Synchrony in a Frustrated Array of Hebb Synapses

The union of the Kuramoto–Sakaguchi model and the Hebb dynamics reproduces the Lisman switch through a bistability in synchronized states. Here, we show that, within certain ranges of the frustration parameter, the chimera pattern can emerge, causing a different, time-evolving, distribution in the Hebbian synaptic strengths. We study the stability range of the chimera as a function of the frustration (phase-lag) parameter. Depending on the range of the frustration, two different types of chimeras can appear spontaneously, i.e., from randomized initial conditions. In the first type, the oscillators in the coherent region rotate, on average, slower than those in the incoherent region; while in the second type, the average rotational frequencies of the two regions are reversed, i.e., the coherent region runs, on average, faster than the incoherent region. We also show that non-stationary behavior at finite N can be controlled by adjusting the natural frequency of a single pacemaker oscillator. By slowly cycling the frequency of the pacemaker, we observe hysteresis in the system. Finally, we discuss how we can have a model for learning and memory.

Exotic states induced by coevolving connection weights and phases in complex networks

We consider an adaptive network, whose connection weights coevolve in congruence with the dynamical states of the local nodes that are under the influence of an external stimulus. The adaptive dynamical system mimics the adaptive synaptic connections common in neuronal networks. The adaptive network under external forcing displays exotic dynamical states such as itinerant chimeras whose population density of coherent and incoherent domains coevolves with the synaptic connection, bump states, and bump frequency cluster states, which do not exist in adaptive networks without forcing. In addition, the adaptive network also exhibits partial synchronization patterns such as phase and frequency clusters, forced entrained, and incoherent states. We introduce two measures for the strength of incoherence based on the standard deviation of the temporally averaged (mean) frequency and on the mean frequency to classify the emergent dynamical states as well as their transitions. We provide a two-parameter phase diagram showing the wealth of dynamical states. We additionally deduce the stability condition for the frequency-entrained state. We use the paradigmatic Kuramoto model of phase oscillators, which is a simple generic model that has been widely employed in unraveling a plethora of cooperative phenomena in natural and man-made systems.

Path-dependent dynamics induced by rewiring networks of inertial oscillators

In networks of coupled oscillators, it is of interest to understand how interaction topology affects synchronization. Many studies have gained key insights into this question by studying the classic Kuramoto oscillator model on static networks. However, new questions arise when the network structure is time varying or when the oscillator system is multistable, the latter of which can occur when an inertial term is added to the Kuramoto model. While the consequences of evolving topology and multistability on collective behavior have been examined separately, real-world systems such as gene regulatory networks and the brain may exhibit these properties simultaneously. It is thus relevant to ask how time-varying network connectivity impacts synchronization in systems that can exhibit multistability. To address this question, we study how the dynamics of coupled Kuramoto oscillators with inertia are affected when the topology of the underlying network changes in time. We show that hysteretic synchronization behavior in networks of coupled inertial oscillators can be driven by changes in connection topology alone. Moreover, we find that certain fixed-density rewiring schemes induce significant changes to the level of global synchrony that remain even after the network returns to its initial configuration, and we show that these changes are robust to a wide range of network perturbations. Our findings highlight that the specific progression of network topology over time, in addition to its initial or final static structure, can play a considerable role in modulating the collective behavior of systems evolving on complex networks.

Modeling Frequency Reduction in Human Groups Performing a Joint Oscillatory Task

In human groups performing oscillatory tasks, it has been observed that the frequency of participants' oscillations reduces when compared to that acquired in solo. This experimental observation is not captured by the standard Kuramoto oscillators, often employed to model human synchronization. In this work, we aim at capturing this observed phenomenon by proposing three alternative modifications of the standard Kuramoto model that are based on three different biologically-relevant hypotheses underlying group synchronization. The three models are tuned, validated and compared against experiments on a group synchronization task, which is a multi-agent extension of the so-called mirror game.

Modeling Tumor Disease and Sepsis by Networks of Adaptively Coupled Phase Oscillators

In this study, we provide a dynamical systems perspective to the modelling of pathological states induced by tumors or infection. A unified disease model is established using the innate immune system as the reference point. We propose a two-layer network model for carcinogenesis and sepsis based upon the interaction of parenchymal cells and immune cells via cytokines, and the co-evolutionary dynamics of parenchymal, immune cells, and cytokines. Our aim is to show that the complex cellular cooperation between parenchyma and stroma (immune layer) in the physiological and pathological case can be qualitatively and functionally described by a simple paradigmatic model of phase oscillators. By this, we explain carcinogenesis, tumor progression, and sepsis by destabilization of the healthy homeostatic state (frequency synchronized), and emergence of a pathological state (desynchronized or multifrequency cluster). The coupled dynamics of parenchymal cells (metabolism) and nonspecific immune cells (reaction of innate immune system) are represented by nodes of a duplex layer. The cytokine interaction is modeled by adaptive coupling weights between the nodes representing the immune cells (with fast adaptation time scale) and the parenchymal cells (slow adaptation time scale) and between the pairs of parenchymal and immune cells in the duplex network (fixed bidirectional coupling). Thereby, carcinogenesis, organ dysfunction in sepsis, and recurrence risk can be described in a correct functional context.

Birth and Stabilization of Phase Clusters by Multiplexing of Adaptive Networks

We propose a concept to generate and stabilize diverse partial synchronization patterns (phase clusters) in adaptive networks which are widespread in neuroscience and social sciences, as well as biology, engineering, and other disciplines. We show by theoretical analysis and computer simulations that multiplexing in a multilayer network with symmetry can induce various stable phase cluster states in a situation where they are not stable or do not even exist in the single layer. Further, we develop a method for the analysis of Laplacian matrices of multiplex networks which allows for insight into the spectral structure of these networks enabling a reduction to the stability problem of single layers. We employ the multiplex decomposition to provide analytic results for the stability of the multilayer patterns. As local dynamics we use the paradigmatic Kuramoto phase oscillator, which is a simple generic model and has been successfully applied in the modeling of synchronization phenomena in a wide range of natural and technological systems.

Effect of repulsive links on frustration in attractively coupled networks

We investigate the impact of attractive-repulsive interaction in networks of limit cycle oscillators. Mainly we focus on the design principle for generating an antiphase state between adjacent nodes in a complex network. We establish that a partial negative control throughout the branches of a spanning tree inside the positively coupled limit cycle oscillators works efficiently well in comparison with randomly chosen negative links to establish zero frustration (antiphase synchronization) in bipartite graphs. Based on the emergence of zero frustration, we develop a universal 0-π rule to understand the antiphase synchronization in a bipartite graph. Further, this rule is used to construct a nonbipartite graph for a given nonzero frustrated value. We finally show the generality of 0-π rule by implementing it in arbitrary undirected nonbipartite graphs of attractive-repulsively coupled limit cycle oscillators and successfully calculate the nonzero frustration value, which matches with numerical data. The validation of the rule is checked through the bifurcation analysis of small networks. Our work may unveil the underlying mechanism of several synchronization phenomena that exist in a network of oscillators having a mixed type of coupling.

Hierarchical frequency clusters in adaptive networks of phase oscillators

Adaptive dynamical networks appear in various real-word systems. One of the simplest phenomenological models for investigating basic properties of adaptive networks is the system of coupled phase oscillators with adaptive couplings. In this paper, we investigate the dynamics of this system. We extend recent results on the appearance of hierarchical frequency multiclusters by investigating the effect of the time scale separation. We show that the slow adaptation in comparison with the fast phase dynamics is necessary for the emergence of the multiclusters and their stability. Additionally, we study the role of double antipodal clusters, which appear to be unstable for all considered parameter values. We show that such states can be observed for a relatively long time, i.e., they are metastable. A geometrical explanation for such an effect is based on the emergence of a heteroclinic orbit.

Effects of synaptic and myelin plasticity on learning in a network of Kuramoto phase oscillators

Models of learning typically focus on synaptic plasticity. However, learning is the result of both synaptic and myelin plasticity. Specifically, synaptic changes often co-occur and interact with myelin changes, leading to complex dynamic interactions between these processes. Here, we investigate the implications of these interactions for the coupling behavior of a system of Kuramoto oscillators. To that end, we construct a fully connected, one-dimensional ring network of phase oscillators whose coupling strength (reflecting synaptic strength) as well as conduction velocity (reflecting myelination) are each regulated by a Hebbian learning rule. We evaluate the behavior of the system in terms of structural (pairwise connection strength and conduction velocity) and functional connectivity (local and global synchronization behavior). We find that adaptive myelination is able to both functionally decouple structurally connected oscillators as well as to functionally couple structurally disconnected oscillators. With regard to the latter, we find that for conditions in which a system limited to synaptic plasticity develops two distinct clusters both structurally and functionally, additional adaptive myelination allows for functional communication across these structural clusters. These results confirm that network states following learning may be different when myelin plasticity is considered in addition to synaptic plasticity, pointing toward the relevance of integrating both factors in computational models of learning.

Beta-Rhythm Oscillations and Synchronization Transition in Network Models of Izhikevich Neurons: Effect of Topology and Synaptic Type

Despite their significant functional roles, beta-band oscillations are least understood. Synchronization in neuronal networks have attracted much attention in recent years with the main focus on transition type. Whether one obtains explosive transition or a continuous transition is an important feature of the neuronal network which can depend on network structure as well as synaptic types. In this study we consider the effect of synaptic interaction (electrical and chemical) as well as structural connectivity on synchronization transition in network models of Izhikevich neurons which spike regularly with beta rhythms. We find a wide range of behavior including continuous transition, explosive transition, as well as lack of global order. The stronger electrical synapses are more conducive to synchronization and can even lead to explosive synchronization. The key network element which determines the order of transition is found to be the clustering coefficient and not the small world effect, or the existence of hubs in a network. These results are in contrast to previous results which use phase oscillator models such as the Kuramoto model. Furthermore, we show that the patterns of synchronization changes when one goes to the gamma band. We attribute such a change to the change in the refractory period of Izhikevich neurons which changes significantly with frequency.

Development of structural correlations and synchronization from adaptive rewiring in networks of Kuramoto oscillators

Synchronization of non-identical oscillators coupled through complex networks is an important example of collective behavior, and it is interesting to ask how the structural organization of network interactions influences this process. Several studies have explored and uncovered optimal topologies for synchronization by making purposeful alterations to a network. On the other hand, the connectivity patterns of many natural systems are often not static, but are rather modulated over time according to their dynamics. However, this co-evolution and the extent to which the dynamics of the individual units can shape the organization of the network itself are less well understood. Here, we study initially randomly connected but locally adaptive networks of Kuramoto oscillators. In particular, the system employs a co-evolutionary rewiring strategy that depends only on the instantaneous, pairwise phase differences of neighboring oscillators, and that conserves the total number of edges, allowing the effects of local reorganization to be isolated. We find that a simple rule-which preserves connections between more out-of-phase oscillators while rewiring connections between more in-phase oscillators-can cause initially disordered networks to organize into more structured topologies that support enhanced synchronization dynamics. We examine how this process unfolds over time, finding a dependence on the intrinsic frequencies of the oscillators, the global coupling, and the network density, in terms of how the adaptive mechanism reorganizes the network and influences the dynamics. Importantly, for large enough coupling and after sufficient adaptation, the resulting networks exhibit interesting characteristics, including degree-frequency and frequency-neighbor frequency correlations. These properties have previously been associated with optimal synchronization or explosive transitions in which the networks were constructed using global information. On the contrary, by considering a time-dependent interplay between structure and dynamics, this work offers a mechanism through which emergent phenomena and organization can arise in complex systems utilizing local rules.

The stability of fixed points for a Kuramoto model with Hebbian interactions

We consider a variation of the Kuramoto model with dynamic coupling, where the coupling strengths are allowed to evolve in response to the phase difference between the oscillators, a model first considered by Ha, Noh, and Park. We demonstrate that the fixed points of this model, as well as their stability, can be completely expressed in terms of the fixed points and stability of the analogous classical Kuramoto problem where the coupling strengths are fixed to a constant (the same for all edges). In particular, for the "all-to-all" network, where the underlying graph is the complete graph, the problem reduces to the problem of understanding the fixed points and stability of the all-to-all Kuramoto model with equal edge weights, a problem that is well understood.

Describing synchronization and topological excitations in arrays of magnetic spin torque oscillators through the Kuramoto model

The collective dynamics in populations of magnetic spin torque oscillators (STO) is an intensely studied topic in modern magnetism. Here, we show that arrays of STO coupled via dipolar fields can be modeled using a variant of the Kuramoto model, a well-known mathematical model in non-linear dynamics. By investigating the collective dynamics in arrays of STO we find that the synchronization in such systems is a finite size effect and show that the critical coupling-for a complete synchronized state-scales with the number of oscillators. Using realistic values of the dipolar coupling strength between STO we show that this imposes an upper limit for the maximum number of oscillators that can be synchronized. Further, we show that the lack of long range order is associated with the formation of topological defects in the phase field similar to the two-dimensional XY model of ferromagnetism. Our results shed new light on the synchronization of STO, where controlling the mutual synchronization of several oscillators is considered crucial for applications.

Phase synchronization of coupled bursting neurons and the generalized Kuramoto model

Bursting neurons fire rapid sequences of action potential spikes followed by a quiescent period. The basic dynamical mechanism of bursting is the slow currents that modulate a fast spiking activity caused by rapid ionic currents. Minimal models of bursting neurons must include both effects. We considered one of these models and its relation with a generalized Kuramoto model, thanks to the definition of a geometrical phase for bursting and a corresponding frequency. We considered neuronal networks with different connection topologies and investigated the transition from a non-synchronized to a partially phase-synchronized state as the coupling strength is varied. The numerically determined critical coupling strength value for this transition to occur is compared with theoretical results valid for the generalized Kuramoto model.