Dynamics of the Kuramoto-Sakaguchi Oscillator Network with Asymmetric Order Parameter

We study the dynamics of a generalized version of the famous Kuramoto-Sakaguchi coupled oscillator model. In the classic version of this system, all oscillators are governed by the same ODE, which depends on the order parameter of the oscillator configuration. The order parameter is the arithmetic mean of the configuration of complex oscillator phases, multiplied by some constant complex coupling factor. In the generalized model we consider, the order parameter is allowed to be any complex linear combination of the complex oscillator phases, so the oscillators are no longer necessarily weighted identically in the order parameter. This asymmetric version of the K-S model exhibits a much richer variety of steady-state dynamical behavior than the classic symmetric version; in addition to stable synchronized states, the system may possess multiple stable (N-1,1) states, in which all but one of the oscillators are in sync, as well as multiple families of neutrally stable asynchronous states or closed orbits, in which no two oscillators are in sync. We present an exhaustive description of the possible steady state dynamical behaviors; our classification depends on the complex coefficients that determine the order parameter. We use techniques from group theory and hyperbolic geometry to reduce the dynamic analysis to a 2D flow on the unit disc, which has geometric significance relative to the hyperbolic metric. The geometric-analytic techniques we develop can in turn be applied to study even more general versions of Kuramoto oscillator networks

Hyperbolic Geometry of Kuramoto Oscillator Networks

Kuramoto oscillator networks have the special property that their trajectories are constrained to lie on the (at most) 3D orbits of the M\"obius group acting on the state space $T^N$ (the $N$-fold torus). This result has been used to explain the existence of the $N-3$ constants of motion discovered by Watanabe and Strogatz for Kuramoto oscillator networks. In this work we investigate geometric consequences of this M\"obius group action. The dynamics of Kuramoto phase models can be further reduced to 2D reduced group orbits, which have a natural geometry equivalent to the unit disk $\Delta$ with the hyperbolic metric. We show that in this metric the original Kuramoto phase model (with order parameter $Z_1$ equal to the centroid of the oscillator configuration of points on the unit circle) is a gradient flow and the model with order parameter $iZ_1$ (corresponding to cosine phase coupling) is a completely integrable Hamiltonian flow. We give necessary and sufficient conditions for general Kuramoto phase models to be gradient or Hamiltonian flows in this metric. This allows us to identify several new infinite families of hyperbolic gradient or Hamiltonian Kuramoto oscillator networks which therefore have simple dynamics with respect to this geometry. We prove that for the $Z_1$ model, a generic 2D reduced group orbit has a unique fixed point corresponding to the hyperbolic barycenter of the oscillator configuration, and therefore the dynamics are equivalent on different generic reduced group orbits. This is not always the case for more general hyperbolic gradient or Hamiltonian flows; the reduced group orbits may have multiple fixed points, which also may bifurcate as the reduced group orbits vary.Comment: Accepted for publication in Journal of Physics A: Mathematical and Theoretica

Dynamical Phase Transitions In Driven Integrate-And-Fire Neurons

We explore the dynamics of an integrate-and-fire neuron with an oscillatory stimulus. The frustration due to the competition between the neuron's natural firing period and that of the oscillatory rhythm, leads to a rich structure of asymptotic phase locking patterns and ordering dynamics. The phase transitions between these states can be classified as either tangent or discontinuous bifurcations, each with its own characteristic scaling laws. The discontinuous bifurcations exhibit a new kind of phase transition that may be viewed as intermediate between continuous and first order, while tangent bifurcations behave like continuous transitions with a diverging coherence scale.Comment: 4 pages, 5 figure

Physics of the rhythmic applause

We discuss in detail a human scale example of the synchronization phenomenon, namely the dynamics of the rhythmic applause. After a detailed experimental investigation, we describe the phenomenon with an approach based on the classical Kuramoto model. Computer simulations based on the theoretical assumptions, reproduce perfectly the observed dynamics. We argue that a frustration present in the system is responsible for the interesting interplay between synchronized and unsynchronized regimesComment: 5 pages, 5 figure

Stability Analysis of Asynchronous States in Neuronal Networks with Conductance-Based Inhibition

Oscillations in networks of inhibitory interneurons have been reported at various sites of the brain and are thought to play a fundamental role in neuronal processing. This Letter provides a self-contained analytical framework that allows numerically efficient calculations of the population activity of a network of conductance-based integrate-and-fire neurons that are coupled through inhibitory synapses. Based on a normalization equation this Letter introduces a novel stability criterion for a network state of asynchronous activity and discusses its perturbations. The analysis shows that, although often neglected, the reversal potential of synaptic inhibition has a strong influence on the stability as well as the frequency of network oscillations

Spatial patterns of desynchronization bursts in networks

We adapt a previous model and analysis method (the {\it master stability function}), extensively used for studying the stability of the synchronous state of networks of identical chaotic oscillators, to the case of oscillators that are similar but not exactly identical. We find that bubbling induced desynchronization bursts occur for some parameter values. These bursts have spatial patterns, which can be predicted from the network connectivity matrix and the unstable periodic orbits embedded in the attractor. We test the analysis of bursts by comparison with numerical experiments. In the case that no bursting occurs, we discuss the deviations from the exactly synchronous state caused by the mismatch between oscillators

Breaking Synchrony by Heterogeneity in Complex Networks

For networks of pulse-coupled oscillators with complex connectivity, we demonstrate that in the presence of coupling heterogeneity precisely timed periodic firing patterns replace the state of global synchrony that exists in homogenous networks only. With increasing disorder, these patterns persist until they reach a critical temporal extent that is of the order of the interaction delay. For stronger disorder these patterns cease to exist and only asynchronous, aperiodic states are observed. We derive self-consistency equations to predict the precise temporal structure of a pattern from the network heterogeneity. Moreover, we show how to design heterogenous coupling architectures to create an arbitrary prescribed pattern.Comment: 4 pages, 3 figure

Time-periodic phases in populations of nonlinearly coupled oscillators with bimodal frequency distributions

The mean field Kuramoto model describing the synchronization of a population of phase oscillators with a bimodal frequency distribution is analyzed (by the method of multiple scales) near regions in its phase diagram corresponding to synchronization to phases with a time periodic order parameter. The richest behavior is found near the tricritical point were the incoherent, stationarily synchronized, ``traveling wave'' and ``standing wave'' phases coexist. The behavior near the tricritical point can be extrapolated to the rest of the phase diagram. Direct Brownian simulation of the model confirms our findings.Comment: Revtex,16 pag.,10 fig., submitted to Physica

Coupled Oscillators with Chemotaxis

A simple coupled oscillator system with chemotaxis is introduced to study morphogenesis of cellular slime molds. The model successfuly explains the migration of pseudoplasmodium which has been experimentally predicted to be lead by cells with higher intrinsic frequencies. Results obtained predict that its velocity attains its maximum value in the interface region between total locking and partial locking and also suggest possible roles played by partial synchrony during multicellular development.Comment: 4 pages, 5 figures, latex using jpsj.sty and epsf.sty, to appear in J. Phys. Soc. Jpn. 67 (1998

Solvable model of a phase oscillator network on a circle with infinite-range Mexican-hat-type interaction

We describe a solvable model of a phase oscillator network on a circle with infinite-range Mexican-hat-type interaction. We derive self-consistent equations of the order parameters and obtain three non-trivial solutions characterized by the rotation number. We also derive relevant characteristics such as the location-dependent distributions of the resultant frequencies of desynchronized oscillators. Simulation results closely agree with the theoretical ones