Generic scaling at the onset of macroscopic mutual entrainment in limit-cycle oscillators with uniform all-to-all coupling

Synchronization, clustering, and weak chimeras in a densely coupled transcription-based oscillator model for split circadian rhythms

The synchronization dynamics for the circadian gene expression in the suprachiasmatic nucleus is investigated using a transcriptional circadian clock gene oscillator model. With global coupling in constant dark (DD) conditions, the model exhibits a one-cluster phase synchronized state, in dim light (dim LL), bistability between one- and two-cluster states and in bright LL, a two-cluster state. The two-cluster phase synchronized state, where some oscillator pairs synchronize in-phase, and some anti-phase, can explain the splitting of the circadian clock, i.e., generation of two bouts of daily activities with certain species, e.g., with hamsters. The one- and two-cluster states can be reached by transferring the animal from DD or bright LL to dim LL, i.e., the circadian synchrony has a memory effect. The stability of the one- and two-cluster states was interpreted analytically by extracting phase models from the ordinary differential equation models. In a modular network with two strongly coupled oscillator populations with weak intragroup coupling, with appropriate initial conditions, one group is synchronized to the one-cluster state and the other group to the two-cluster state, resulting in a weak-chimera state. Computational modeling suggests that the daily rhythms in sleep-wake depend on light intensity acting on bilateral networks of suprachiasmatic nucleus (SCN) oscillators. Addition of a network heterogeneity (coupling between the left and right SCN) allowed the system to exhibit chimera states. The simulations can guide experiments in the circadian rhythm research to explore the effect of light intensity on the complexities of circadian desynchronization.

Neuronal avalanches: Sandpiles of self-organized criticality or critical dynamics of brain waves?

Analytical expressions for scaling of brain wave spectra derived from the general nonlinear wave Hamiltonian form show excellent agreement with experimental "neuronal avalanche" data. The theory of the weakly evanescent nonlinear brain wave dynamics [Phys. Rev. Research 2, 023061 (2020); J. Cognitive Neurosci. 32, 2178 (2020)] reveals the underlying collective processes hidden behind the phenomenological statistical description of the neuronal avalanches and connects together the whole range of brain activity states, from oscillatory wave-like modes, to neuronal avalanches, to incoherent spiking, showing that the neuronal avalanches are just the manifestation of the different nonlinear side of wave processes abundant in cortical tissue. In a more broad way these results show that a system of wave modes interacting through all possible combinations of the third order nonlinear terms described by a general wave Hamiltonian necessarily produces anharmonic wave modes with temporal and spatial scaling properties that follow scale free power laws. To the best of our knowledge this has never been reported in the physical literature and may be applicable to many physical systems that involve wave processes and not just to neuronal avalanches.

Synchronization transitions in Kuramoto networks with higher-mode interaction

Synchronization is an omnipresent collective phenomenon in nature and technology, whose understanding is still elusive for real-world systems in particular. We study the synchronization transition in a phase oscillator system with two nonvanishing Fourier-modes in the interaction function, hence going beyond the Kuramoto paradigm. We show that the transition scenarios crucially depend on the interplay of the two coupling modes. We describe the multistability induced by the presence of a second coupling mode. By extending the collective coordinate approach, we describe the emergence of various states observed in the transition from incoherence to coherence. Remarkably, our analysis suggests that, in essence, the two-mode coupling gives rise to states characterized by two independent but interacting groups of oscillators. We believe that these findings will stimulate future research on dynamical systems, including complex interaction functions beyond the Kuramoto-type.

Inhibitory neurons control the consolidation of neural assemblies via adaptation to selective stimuli

Brain circuits display modular architecture at different scales of organization. Such neural assemblies are typically associated to functional specialization but the mechanisms leading to their emergence and consolidation still remain elusive. In this paper we investigate the role of inhibition in structuring new neural assemblies driven by the entrainment to various inputs. In particular, we focus on the role of partially synchronized dynamics for the creation and maintenance of structural modules in neural circuits by considering a network of excitatory and inhibitory [Formula: see text]-neurons with plastic Hebbian synapses. The learning process consists of an entrainment to temporally alternating stimuli that are applied to separate regions of the network. This entrainment leads to the emergence of modular structures. Contrary to common practice in artificial neural networks-where the acquired weights are typically frozen after the learning session-we allow for synaptic adaptation even after the learning phase. We find that the presence of inhibitory neurons in the network is crucial for the emergence and the post-learning consolidation of the modular structures. Indeed networks made of purely excitatory neurons or of neurons not respecting Dale's principle are unable to form or to maintain the modular architecture induced by the stimuli. We also demonstrate that the number of inhibitory neurons in the network is directly related to the maximal number of neural assemblies that can be consolidated, supporting the idea that inhibition has a direct impact on the memory capacity of the neural network.

Critical brain wave dynamics of neuronal avalanches

Analytical expressions for scaling of brain wave spectra derived from the general non-linear wave Hamiltonian form show excellent agreement with experimental "neuronal avalanche" data. The theory of the weakly evanescent non-linear brain wave dynamics reveals the underlying collective processes hidden behind the phenomenological statistical description of the neuronal avalanches and connects together the whole range of brain activity states, from oscillatory wave-like modes, to neuronal avalanches, to incoherent spiking, showing that the neuronal avalanches are just the manifestation of the different non-linear side of wave processes abundant in cortical tissue. In a more broad way these results show that a system of wave modes interacting through all possible combinations of the third order non-linear terms described by a general wave Hamiltonian necessarily produces anharmonic wave modes with temporal and spatial scaling properties that follow scale free power laws. To the best of our knowledge this has never been reported in the physical literature and may be applicable to many physical systems that involve wave processes and not just to neuronal avalanches.

Critical Parameters in Dynamic Network Modeling of Sepsis

In this work, we propose a dynamical systems perspective on the modeling of sepsis and its organ-damaging consequences. We develop a functional two-layer network model for sepsis based upon the interaction of parenchymal cells and immune cells via cytokines, and the coevolutionary dynamics of parenchymal, immune cells, and cytokines. By means of the simple paradigmatic model of phase oscillators in a two-layer system, we analyze the emergence of organ threatening interactions between the dysregulated immune system and the parenchyma. We demonstrate that the complex cellular cooperation between parenchyma and stroma (immune layer) either in the physiological or in the pathological case can be related to dynamical patterns of the network. In this way we explain sepsis by the dysregulation of the healthy homeostatic state (frequency synchronized) leading to a pathological state (desynchronized or multifrequency cluster) in the parenchyma. We provide insight into the complex stabilizing and destabilizing interplay of parenchyma and stroma by determining critical interaction parameters. The coupled dynamics of parenchymal cells (metabolism) and nonspecific immune cells (response of the innate immune system) is represented by nodes of a duplex layer. Cytokine interaction is modeled by adaptive coupling weights between nodes representing immune cells (with fast adaptation timescale) and parenchymal cells (slow adaptation timescale), and between pairs of parenchymal and immune cells in the duplex network (fixed bidirectional coupling). The proposed model allows for a functional description of organ dysfunction in sepsis and the recurrence risk in a plausible pathophysiological context.

Controlling permeation in electrically deforming liquid crystal network films: A dynamical Landau theory

Liquid crystal networks exploit the coupling between the responsivity of liquid crystalline mesogens, e.g., to electric fields, and the (visco)elastic properties of a polymer network. Because of this, these materials have been put forward for a wide array of applications, including responsive surfaces such as artificial skins and membranes. For such applications, the desired functional response must generally be realized under strict geometrical constraints, such as provided by supported thin films. To model such settings, we present a dynamical, spatially heterogeneous Landau-type theory for electrically actuated liquid crystal network films. We find that the response of the liquid crystal network permeates the film from top to bottom, and illustrate how this affects the timescale associated with macroscopic deformation. Finally, by linking our model parameters to experimental quantities, we suggest that the permeation rate can be controlled by varying the aspect ratio of the mesogens and their degree of orientational order when crosslinked into the polymer network, for which we predict a single optimum. Our results contribute specifically to the rational design of future applications involving transport or on-demand release of molecular cargo in liquid crystal network films.

Generalized splay states in phase oscillator networks

Networks of coupled phase oscillators play an important role in the analysis of emergent collective phenomena. In this article, we introduce generalized m-splay states constituting a special subclass of phase-locked states with vanishing mth order parameter. Such states typically manifest incoherent dynamics, and they often create high-dimensional families of solutions (splay manifolds). For a general class of phase oscillator networks, we provide explicit linear stability conditions for splay states and exemplify our results with the well-known Kuramoto-Sakaguchi model. Importantly, our stability conditions are expressed in terms of just a few observables such as the order parameter or the trace of the Jacobian. As a result, these conditions are simple and applicable to networks of arbitrary size. We generalize our findings to phase oscillators with inertia and adaptively coupled phase oscillator models.

Critical exponents in coupled phase-oscillator models on small-world networks

A coupled phase-oscillator model consists of phase oscillators, each of which has the natural frequency obeying a probability distribution and couples with other oscillators through a given periodic coupling function. This type of model is widely studied since it describes the synchronization transition, which emerges between the nonsynchronized state and partially synchronized states. The synchronization transition is characterized by several critical exponents, and we focus on the critical exponent defined by coupling strength dependence of the order parameter for revealing universality classes. In a typical interaction represented by the perfect graph, an infinite number of universality classes is yielded by dependency on the natural frequency distribution and the coupling function. Since the synchronization transition is also observed in a model on a small-world network, whose number of links is proportional to the number of oscillators, a natural question is whether the infinite number of universality classes remains in small-world networks irrespective of the order of links. Our numerical results suggest that the number of universality classes is reduced to one and the critical exponent is shared in the considered models having coupling functions up to second harmonics with unimodal and symmetric natural frequency distributions.

Connecting empirical phenomena and theoretical models of biological coordination across scales

Coordination in living systems-from cells to people-must be understood at multiple levels of description. Analyses and modelling of empirically observed patterns of biological coordination often focus either on ensemble-level statistics in large-scale systems with many components, or on detailed dynamics in small-scale systems with few components. The two approaches have proceeded largely independent of each other. To bridge this gap between levels and scales, we have recently conducted a human experiment of mid-scale social coordination specifically designed to reveal coordination at multiple levels (ensemble, subgroups and dyads) simultaneously. Based on this experiment, the present work shows that, surprisingly, a single system of equations captures key observations at all relevant levels. It also connects empirically validated models of large- and small-scale biological coordination-the Kuramoto and extended Haken-Kelso-Bunz (HKB) models-and the hallmark phenomena that each is known to capture. For example, it exhibits both multistability and metastability observed in small-scale empirical research (via the second-order coupling and symmetry breaking in extended HKB) and the growth of biological complexity as a function of scale (via the scalability of the Kuramoto model). Only by incorporating both of these features simultaneously can we reproduce the essential coordination behaviour observed in our experiment.

Bus bunching as a synchronisation phenomenon

Bus bunching is a perennial phenomenon that not only diminishes the efficiency of a bus system, but also prevents transit authorities from keeping buses on schedule. We present a physical theory of buses serving a loop of bus stops as a ring of coupled self-oscillators, analogous to the Kuramoto model. Sustained bunching is a repercussion of the process of phase synchronisation whereby the phases of the oscillators are locked to each other. This emerges when demand exceeds a critical threshold. Buses also bunch at low demand, albeit temporarily, due to frequency detuning arising from different human drivers’ distinct natural speeds. We calculate the critical transition when complete phase locking (full synchronisation) occurs for the bus system, and posit the critical transition to completely no phase locking (zero synchronisation). The intermediate regime is the phase where clusters of partially phase locked buses exist. Intriguingly, these theoretical results are in close correspondence to real buses in a university’s shuttle bus system.

Collective behavior of oscillating electric dipoles

We investigate the dynamics of a population of identical biomolecules mimicked as electric dipoles with random orientations and positions in space and oscillating with their intrinsic frequencies. The biomolecules, beyond being coupled among themselves via the dipolar interaction, are also driven by a common external energy supply. A collective mode emerges by decreasing the average distance among the molecules as testified by the emergence of a clear peak in the power spectrum of the total dipole moment. This is due to a coherent vibration of the most part of the molecules at a frequency definitely larger than their own frequencies corresponding to a partial cluster synchronization of the biomolecules. These results can be verified experimentally via spectroscopic investigations of the strength of the intermolecular electrodynamic interactions, thus being able to test the possible biological relevance of the observed macroscopic mode.

Landau-like theory for universality of critical exponents in quasistationary states of isolated mean-field systems

An external force dynamically drives an isolated mean-field Hamiltonian system to a long-lasting quasistationary state, whose lifetime increases with population of the system. For second order phase transitions in quasistationary states, two nonclassical critical exponents have been reported individually by using a linear and a nonlinear response theories in a toy model. We provide a simple way to compute the critical exponents all at once, which is an analog of the Landau theory. The present theory extends the universality class of the nonclassical exponents to spatially periodic one-dimensional systems and shows that the exponents satisfy a classical scaling relation inevitably by using a key scaling of momentum.

Interaction function of oscillating coupled neurons

Large scale simulations of electrically coupled neuronal oscillators often employ the phase coupled oscillator paradigm to understand and predict network behavior. We study the nature of the interaction between such coupled oscillators using weakly coupled oscillator theory. By employing piecewise linear approximations for phase response curves and voltage time courses and parametrizing their shapes, we compute the interaction function for all such possible shapes and express it in terms of discrete Fourier modes. We find that reasonably good approximation is achieved with four Fourier modes that comprise of both sine and cosine terms.

Nonstandard scaling law of fluctuations in finite-size systems of globally coupled oscillators

Universal scaling laws form one of the central issues in physics. A nonstandard scaling law or a breakdown of a standard scaling law, on the other hand, can often lead to the finding of a new universality class in physical systems. Recently, we found that a statistical quantity related to fluctuations follows a nonstandard scaling law with respect to the system size in a synchronized state of globally coupled nonidentical phase oscillators [I. Nishikawa et al., Chaos 22, 013133 (2012)]. However, it is still unclear how widely this nonstandard scaling law is observed. In the present paper, we discuss the conditions required for the unusual scaling law in globally coupled oscillator systems and validate the conditions by numerical simulations of several different models.

Autonomous and forced dynamics of oscillator ensembles with global nonlinear coupling: an experimental study

We perform experiments with 72 electronic limit-cycle oscillators, globally coupled via a linear or nonlinear feedback loop. While in the linear case we observe a standard Kuramoto-like synchronization transition, in the nonlinear case, with increase of the coupling strength, we first observe a transition to full synchrony and then a desynchronization transition to a quasiperiodic state. However, in this state the ensemble remains coherent so that the amplitude of the mean field is nonzero, but the frequency of the mean field is larger than frequencies of all oscillators. Next, we analyze effects of common periodic forcing of the linearly or nonlinearly coupled ensemble and demonstrate regimes when the mean field is entrained by the force whereas the oscillators are not.

Hierarchical synchrony of phase oscillators in modular networks

We study synchronization of sinusoidally coupled phase oscillators on networks with modular structure and a large number of oscillators in each community. Of particular interest is the hierarchy of local and global synchrony, i.e., synchrony within and between communities, respectively. Using the recent ansatz of Ott and Antonsen [Chaos 18, 037113 (2008)], we find that the degree of local synchrony can be determined from a set of coupled low-dimensional equations. If the number of communities in the network is large, a low-dimensional description of global synchrony can be also found. Using these results, we study bifurcations between different types of synchrony. We find that, depending on the relative strength of local and global coupling, the transition to synchrony in the network can be mediated by local or global effects.

Cluster synchrony in systems of coupled phase oscillators with higher-order coupling

We study the phenomenon of cluster synchrony that occurs in ensembles of coupled phase oscillators when higher-order modes dominate the coupling between oscillators. For the first time, we develop a complete analytic description of the dynamics in the limit of a large number of oscillators and use it to quantify the degree of cluster synchrony, cluster asymmetry, and switching. We use a variation of the recent dimensionality-reduction technique of Ott and Antonsen [Chaos 18, 037113 (2008)] and find an analytic description of the degree of cluster synchrony valid on a globally attracting manifold. Shaped by this manifold, there is an infinite family of steady-state distributions of oscillators, resulting in a high degree of multistability in the cluster asymmetry. We also show how through external forcing the degree of asymmetry can be controlled, and suggest that systems displaying cluster synchrony can be used to encode and store data.

Many-body theory of synchronization by long-range interactions

Synchronization of coupled oscillators on a d-dimensional lattice with the power-law coupling G(r) = g0/rα and randomly distributed intrinsic frequency is analyzed. A systematic perturbation theory is developed to calculate the order parameter profile and correlation functions in powers of ϵ = α/d-1. For α ≤ d, the system exhibits a sharp synchronization transition as described by the conventional mean-field theory. For α > d, the transition is smeared by the quenched disorder, and the macroscopic order parameter ψ decays slowly with g0 as |ψ| ∝ g(0)(2).

Simple model for bursting dynamics of neurons

Neuronal cells in isolation or as an assembly exhibit bursting behavior on two different time scales. We introduce a simple one-dimensional model which requires only one phase variable to describe the phenomenon of parabolic bursting. The analysis in the continuum limit reveals that for any unimodal distribution of frequencies, the qualitative properties of the full and the reduced model are identical. Further, we derive analytically an exact low-dimensional description of a globally coupled network of bursting oscillators for our model. Study of the stability for this low-dimensional model reveals different dynamical signatures in the parameter space. We demonstrate that the structure of the parameter space remains independent of the number of spikes per burst.

Renormalization of oscillator lattices with disorder

A real-space renormalization transformation is constructed for lattices of nonidentical oscillators with dynamics of the general form dvarphi_{k}/dt=omega_{k}+g summation operator_{l}f_{lk}(varphi_{l},varphi_{k}) . The transformation acts on ensembles of such lattices. Critical properties corresponding to a second-order phase transition toward macroscopic synchronization are deduced. The analysis is potentially exact but relies in part on unproven assumptions. Numerically, second-order phase transitions with the predicted properties are observed as g increases in two structurally different two-dimensional oscillator models. One model has smooth coupling f_{lk}(varphi_{l},varphi_{k})=phi(varphi_{l}-varphi_{k}) , where phi(x) is nonodd. The other model is pulse coupled, with f_{lk}(varphi_{l},varphi_{k})=delta(varphi_{l})phi(varphi_{k}) . Lower bounds for the critical dimensions for different types of coupling are obtained. For nonodd coupling, macroscopic synchronization cannot be ruled out for any dimension D> or =1 , whereas in the case of odd coupling, the well-known result that it can be ruled out for D<3 is regained.

Stability diagram for the forced Kuramoto model

We analyze the periodically forced Kuramoto model. This system consists of an infinite population of phase oscillators with random intrinsic frequencies, global sinusoidal coupling, and external sinusoidal forcing. It represents an idealization of many phenomena in physics, chemistry, and biology in which mutual synchronization competes with forced synchronization. In other words, the oscillators in the population try to synchronize with one another while also trying to lock onto an external drive. Previous work on the forced Kuramoto model uncovered two main types of attractors, called forced entrainment and mutual entrainment, but the details of the bifurcations between them were unclear. Here we present a complete bifurcation analysis of the model for a special case in which the infinite-dimensional dynamics collapse to a two-dimensional system. Exact results are obtained for the locations of Hopf, saddle-node, and Takens-Bogdanov bifurcations. The resulting stability diagram bears a striking resemblance to that for the weakly nonlinear forced van der Pol oscillator.

Resonances and entrainment breakup in Kuramoto models with multimodal frequency densities

We characterize some intriguing aspects of the entrainment behavior of coupled oscillators by means of a perturbation analysis of the partially synchronized solution of the classical Kuramoto-Sakaguchi model. The analysis reveals that partial entrainment may disappear with increasing coupling strength. It also predicts the occurrence of resonances: partial entrainment is induced in oscillators with natural frequencies in specific intervals not corresponding to high oscillator densities. The results are illustrated by simulations.

Continuous and discontinuous phase transitions and partial synchronization in stochastic three-state oscillators

We investigate both continuous (second-order) and discontinuous (first-order) transitions to macroscopic synchronization within a single class of discrete, stochastic (globally) phase-coupled oscillators. We provide analytical and numerical evidence that the continuity of the transition depends on the coupling coefficients and, in some nonuniform populations, on the degree of quenched disorder. Hence, in a relatively simple setting this class of models exhibits the qualitative behaviors characteristic of a variety of considerably more complicated models. In addition, we study the microscopic basis of synchronization above threshold and detail the counterintuitive subtleties relating measurements of time-averaged frequencies and mean-field oscillations. Most notably, we observe a state of suprathreshold partial synchronization in which time-averaged frequency measurements from individual oscillators do not correspond to the frequency of macroscopic oscillations observed in the population.

Chaotic spatial bifurcation by complex coupling

A spatial bifurcation (a transition from stationary to oscillatory regime) in a chain of unidirectionally coupled phase systems is studied. It is shown that complication of coupling terms can make this bifurcation spatially chaotic in contrast to the previously observed "regular" and "predictable" type. It is demonstrated that the found type of spatial bifurcation corresponds to a smooth (predictable) manifold in the parameter space, while its spatial location gets actually unpredictable being governed by regularities of chaotic behavior. We infer that complex collective dynamics may arise in networks with plain architecture and simple dynamics of individual elements if nontrivial coupling is realized.

Phase diagram of a generalized Winfree model

We study the phase diagram of a generalized Winfree model. The modification is such that the coupling depends on the fraction of synchronized oscillators, a situation which has been noted in some experiments on coupled Josephson junctions and mechanical systems. We let the global coupling k be a function of the Kuramoto order parameter r through an exponent z such that z=1 corresponds to the standard Winfree model, z<1 strengthens the coupling at low r (low amount of synchronization), and at z>1 , the coupling is weakened for low r . Using both analytical and numerical approaches, we find that z controls the size of the incoherent phase region and that one may make the incoherent behavior less typical by choosing z<1 . We also find that the original Winfree model is a rather special case; indeed, the partial locked behavior disappears for z>1 . At fixed k and varying gamma , the stability boundary of the locked phase corresponds to a transition that is continuous for z<1 and first order for z>1 . This change in the nature of the transition is in accordance with a previous study of a similarly modified Kuramoto model.

Chimera Ising walls in forced nonlocally coupled oscillators

Nonlocally coupled oscillator systems can exhibit an exotic spatiotemporal structure called a chimera, where the system splits into two groups of oscillators with sharp boundaries, one of which is phase locked and the other phase randomized. Two examples of chimera states are known: the first one appears in a ring of phase oscillators, and the second is associated with two-dimensional rotating spiral waves. In this paper, we report yet another example of the chimera state that is associated with the so-called Ising walls in one-dimensional spatially extended systems. This chimera state is exhibited by a nonlocally coupled complex Ginzburg-Landau equation with external forcing.

Singular unlocking transition in the Winfree model of coupled oscillators

The Winfree model consists of a population of globally coupled phase oscillators with randomly distributed natural frequencies. As the coupling strength and the spread of natural frequencies are varied, the various stable states of the model can undergo bifurcations, nearly all of which have been characterized previously. The one exception is the unlocking transition, in which the frequency-locked state disappears abruptly as the spread of natural frequencies exceeds a critical width. Viewed as a function of the coupling strength, this critical width defines a bifurcation curve in parameter space. For the special case where the frequency distribution is uniform, earlier work had uncovered a puzzling singularity in this bifurcation curve. Here we seek to understand what causes the singularity. Using the Poincaré-Lindstedt method of perturbation theory, we analyze the locked state and its associated unlocking transition, first for an arbitrary distribution of natural frequencies, and then for discrete systems of N oscillators. We confirm that the bifurcation curve becomes singular for a continuum uniform distribution, yet find that it remains well behaved for any finite N , suggesting that the continuum limit is responsible for the singularity.

Method for determining a coupling function in coupled oscillators with application to Belousov-Zhabotinsky oscillators

A coupling function that describes the interaction between self-sustained oscillators in a phase equation is derived and applied experimentally to Belousov-Zhabotinsky (BZ) oscillators. It is demonstrated that the synchronous behavior of coupled BZ reactors is explained extremely well in terms of the coupling function thus obtained. This method does not require comprehensive knowledge of either the oscillation mechanism or the interaction among the oscillators, both of these being often difficult to elucidate in an actual system. These facts enable us to accurately analyze the weakly coupled entrainment phenomenon through the direct measurement of the coupling function.

Determination of a coupling function in multicoupled oscillators

A new method to determine a coupling function in a phase model is theoretically derived for coupled self-sustained oscillators and applied to Belousov-Zhabotinsky (BZ) oscillators. The synchronous behavior of two coupled BZ reactors is explained extremely well in terms of the coupling function thus obtained. This method is expected to be applicable to weakly coupled multioscillator systems, in which mutual coupling among nearly identical oscillators occurs in a similar manner. The importance of higher-order harmonic terms involved in the coupling function is also discussed.

Coarse graining the dynamics of coupled oscillators

We present an equation-free computational approach to the study of the coarse-grained dynamics of finite assemblies of nonidentical coupled oscillators at and near full synchronization. We use coarse-grained observables which account for the (rapidly developing) correlations between phase angles and natural frequencies. Exploiting short bursts of appropriately initialized detailed simulations, we circumvent the derivation of closures for the long-term dynamics of the assembly statistics.

Repulsive synchronization in an array of phase oscillators

We study the dynamics of a repulsively coupled array of phase oscillators. For an array of globally coupled identical oscillators, repulsive coupling results in a family of synchronized regimes characterized by zero mean field. If the number of oscillators is sufficiently large, phase locking among oscillators is destroyed, independently of the coupling strength, when the oscillators' natural frequencies are not the same. In locally coupled networks, however, phase locking occurs even for nonidentical oscillators when the coupling strength is sufficiently strong.

Noise-controlled oscillations and their bifurcations in coupled phase oscillators

We derive in Gaussian approximation dynamical equations for the first two cumulants of the mean field fluctuations in a system of globally coupled stochastic phase oscillators. In these equations the intensity of noise serves as an explicit control parameter. Its variation generates transitions between three dynamical regimes: (i) stationary, (ii) rotatory and (iii) locally oscillatory (breathing). The latter regime has previously not been reported in studies of globally coupled noisy phase oscillators. Our detailed bifurcation analysis is supported by numerical simulations of an ensemble of coupled stochastic phase oscillators. Similar regimes are also found in simulations of globally coupled stochastic FitzHugh-Nagumo elements.

Phase diagram for the Winfree model of coupled nonlinear oscillators

In 1967 Winfree proposed a mean-field model for the spontaneous synchronization of chorusing crickets, flashing fireflies, circadian pacemaker cells, or other large populations of biological oscillators. Here we give the first bifurcation analysis of the model, for a tractable special case. The system displays rich collective dynamics as a function of the coupling strength and the spread of natural frequencies. Besides incoherence, frequency locking, and oscillator death, there exist hybrid solutions that combine two or more of these states. We present the phase diagram and derive several of the stability boundaries analytically.