Conformists and contrarians in a Kuramoto model with identical natural frequencies

Dynamical phase transitions in the nonreciprocal Ising model

Nonreciprocal interactions in many-body systems lead to time-dependent states, commonly observed in biological, chemical, and ecological systems. The stability of these states in the thermodynamic limit and the critical behavior of the phase transition from static to time-dependent states are not yet fully understood. To address these questions, we study a minimalistic system endowed with nonreciprocal interactions: an Ising model with two spin species having opposing goals. The mean-field equation predicts three stable phases: disorder, static order, and a time-dependent swap phase. Large-scale numerical simulations support the following: (i) in two dimensions, the swap phase is destabilized by defects; (ii) in three dimensions, the swap phase is stable and has the properties of a time crystal; (iii) the transition from disorder to swap in three dimensions is characterized by the critical exponents of the 3D XY model and corresponds to the breaking of a continuous symmetry, time translation invariance; (iv) when the two species have fully antisymmetric couplings, the static-order phase is unstable in any finite dimension due to droplet growth; and (v) in the general case of asymmetric couplings, static order can be restored by a droplet-capture mechanism preventing the droplets from growing indefinitely. We provide details on the full phase diagram, which includes first- and second-order-like phase transitions, and study how the system coarsens into swap and static-order states.

Introduction to Focus Issue: Data-driven models and analysis of complex systems

This Focus Issue highlights recent advances in the study of complex systems, with a particular emphasis on data-driven research. Our editorial explores a diverse array of topics, including financial markets, electricity pricing, power grids, lasers, the Earth's climate, hydrology, neuronal assemblies and the brain, biomedicine, complex networks, real-world hypergraphs, animal behavior, and social media. This diversity underscores the broad applicability of complex systems research. Here, we summarize the 47 published works under this Focus Issue, which employ state-of-the-art or novel methodologies in machine learning, higher-order correlations, control theory, embeddings, information theory, symmetry analysis, extreme event modeling, time series analysis, fractal techniques, Markov chains, and persistent homology, to name a few. These methods have substantially enhanced our understanding of the intricate dynamics of complex systems. Furthermore, the published works demonstrate the potential of data-driven approaches to revolutionize the study of complex systems, paving the way for future research directions and breakthroughs at the intersection of complexity science and the digital era of data.

Calibrating coupling for collaboration with Kuramoto

We calibrate a Kuramoto-model-inspired representation of peer-to-peer collaboration using data on the maximum team size where coordination breaks down. The Kuramoto model is modified, normalizing the coupling by the degree of input and output nodes, reflecting dispersion of cognitive resources in both absorbing incoming- and tracking outgoing-information. We find a critical point, with loss of synchronization as the number of nodes grows and analytically determine this point and calibrate the coupling with the known maximum team size. We test that against the known "span of control" for a leader/supervisor organization. Our results suggest larger maximum team sizes than early management science proposes, but are consistent with studies that focus only the relationship between supervisor and subordinates, excluding other internal interactions in the team.

Chimera-inspired dynamics: When higher-order interactions are expressed differently

The exploration of chimera-inspired dynamics in nonlocally coupled networks of Kuramoto oscillators with higher-order interactions is still in its nascent stages. Concurrently, the investigation of collective phenomena in higher-order interaction networks is gaining attraction. Here, we observe that hypergraph networks tend to synchronize through lower-order interactions, whereas simplicial complex networks exhibit a preference for higher-order interactions. This observation suggests that higher-order representations manifest substantial differences in chimera-inspired synchronization regions. Moreover, we introduce an explicit expression for identifying the chimera state. With a comprehensive basin stability analysis and the interplay of pairwise and higher-order interaction strengths, the emergence of the chimera state is inherent in high-order interaction networks. Our findings contribute to the understanding of chimera-inspired dynamics in higher-order interaction networks.

Phase holonomy underlies puzzling temporal patterns in Kuramoto models with two sub-populations

We present a geometric investigation of curious dynamical behaviors previously reported in Kuramoto models with two sub-populations. Our study demonstrates that chimeras and traveling waves in such models are associated with the birth of geometric phase. Although manifestations of geometric phase are frequent in various fields of physics, this is the first time (to our best knowledge) that such a phenomenon is exposed in ensembles of Kuramoto oscillators or, more broadly, in complex systems.

Nature of the Volcano Transition in the Fully Disordered Kuramoto Model

Randomly coupled phase oscillators may synchronize into disordered patterns of collective motion. We analyze this transition in a large, fully connected Kuramoto model with symmetric but otherwise independent random interactions. Using the dynamical cavity method, we reduce the dynamics to a stochastic single-oscillator problem with self-consistent correlation and response functions that we study analytically and numerically. We clarify the nature of the volcano transition and elucidate its relation to the existence of an oscillator glass phase.

Cluster formation due to repulsive spanning trees in attractively coupled networks

Ensembles of coupled nonlinear oscillators are a popular paradigm and an ideal benchmark for analyzing complex collective behaviors. The onset of cluster synchronization is found to be at the core of various technological and biological processes. The current literature has investigated cluster synchronization by focusing mostly on the case of attractive coupling among the oscillators. However, the case of two coexisting competing interactions is of practical interest due to their relevance in diverse natural settings, including neuronal networks consisting of excitatory and inhibitory neurons, the coevolving social model with voters of opposite opinions, and ecological plant communities with both facilitation and competition, to name a few. In the present article, we investigate the impact of repulsive spanning trees on cluster formation within a connected network of attractively coupled limit-cycle oscillators. We successfully predict which nodes belong to each cluster and the emergent frustration of the connected networks independent of the particular local dynamics at the network nodes. We also determine local asymptotic stability of the cluster states using an approach based on the formulation of a master stability function. We additionally validate the emergence of solitary states and antisynchronization for some specific choices of spanning trees and networks.

Chimeras and traveling waves in ensembles of Kuramoto oscillators off the Poisson manifold

We examine how perturbations off the Poisson manifold affect chimeras and traveling waves (TWs) in Kuramoto models with two sub-populations. Our numerical study is based on simulations on invariant manifolds, which contain von Mises probability distributions. Our study demonstrates that chimeras and TWs off the Poisson manifold always "breathe", and the effect of breathing is more pronounced further from the Poisson manifold. On the other side, TWs arising in similar models on the sphere always breathe moderately, no matter if the dynamics take place near the Poisson manifold or far away from it.

Directional synchrony among self-propelled particles under spatial influence

Synchronization is one of the emerging collective phenomena in interacting particle systems. Its ubiquitous presence in nature, science, and technology has fascinated the scientific community over the decades. Moreover, a great deal of research has been, and is still being, devoted to understand various physical aspects of the subject. In particular, the study of interacting active particles has led to exotic phase transitions in such systems which have opened up a new research front-line. Motivated by this line of work, in this paper, we study the directional synchrony among self-propelled particles. These particles move inside a bounded region, and crucially their directions are also coupled with spatial degrees of freedom. We assume that the directional coupling between two particles is influenced by the relative spatial distance which changes over time. Furthermore, the nature of the influence is considered to be both short and long-ranged. We explore the phase transition scenario in both the cases and propose an approximation technique which enables us to analytically find the critical transition point. The results are further supported with numerical simulations. Our results have potential importance in the study of active systems like bird flocks, fish schools, and swarming robots where spatial influence plays a pertinent role.

Deterministic correlations enhance synchronization in oscillator populations with heterogeneous coupling

Synchronization is a critical phenomenon that displays a pivotal role in a wealth of dynamical processes ranging from natural to artificial systems. Here, we untangle the synchronization optimization in a system of globally coupled phase oscillators incorporating heterogeneous interactions encoded by the deterministic-random coupling. We uncover that, within the given restriction, the added deterministic correlations can profoundly enhance the synchronizability in comparison with the uncorrelated scenario. The critical points manifesting the onset of synchronization and desynchronization transitions, as well as the level of phase coherence, are significantly shaped by the increment of deterministic correlations. In particular, we provide an analytical treatment to properly ground the mechanism underlying synchronization enhancement and substantiate that the analytical predictions are in fair agreement with the numerical simulations. This study is a step forward in highlighting the importance of heterogeneous coupling among dynamical agents, which provides insights for control strategies of synchronization in complex systems.

On the number of stable solutions in the Kuramoto model

We consider a system of n coupled oscillators described by the Kuramoto model with the dynamics given by θ˙=ω+Kf(θ). In this system, an equilibrium solution θ∗ is considered stable when ω+Kf(θ∗)=0, and the Jacobian matrix Df(θ∗) has a simple eigenvalue of zero, indicating the presence of a direction in which the oscillators can adjust their phases. Additionally, the remaining eigenvalues of Df(θ∗) are negative, indicating stability in orthogonal directions. A crucial constraint imposed on the equilibrium solution is that |Γ(θ∗)|≤π, where |Γ(θ∗)| represents the length of the shortest arc on the unit circle that contains the equilibrium solution θ∗. We provide a proof that there exists a unique solution satisfying the aforementioned stability criteria. This analysis enhances our understanding of the stability and uniqueness of these solutions, offering valuable insights into the dynamics of coupled oscillators in this system.

Abrupt symmetry-preserving transition from the chimera state

We consider two populations of the globally coupled Sakaguchi-Kuramoto model with the same intra- and interpopulations coupling strengths. The oscillators constituting the intrapopulation are identical whereas the interpopulations are nonidentical with a frequency mismatch. The asymmetry parameters ensure the permutation symmetry among the oscillators constituting the intrapopulation and a reflection symmetry among the oscillators constituting the interpopulation. We show that the chimera state manifests by spontaneously breaking the reflection symmetry and also exists in almost in the entire explored range of the asymmetry parameter without restricting to the near π/2 values of it. The saddle-node bifurcation mediates the abrupt transition from the symmetry breaking chimera state to the symmetry-preserving synchronized oscillatory state in the reverse trace, whereas the homoclinic bifurcation mediates the transition from the synchronized oscillatory state to synchronized steady state in the forward trace. We deduce the governing equations of motion for the macroscopic order parameters employing the finite-dimensional reduction by Watanabe and Strogatz. The analytical saddle-node and homoclinic bifurcation conditions agree well with the simulations results and the bifurcation curves.

`Friend or foe' and decision making initiative in complex conflict environments

We present a novel mathematical model of two adversarial forces in the vicinity of a non-combatant population in order to explore the impact of each force pursuing specific decision-making strategies. Each force has the opportunity to draw support by enabling the decision-making initiative of the population, in tension with maintaining tactical and organisational effectiveness over their adversary. Each dynamic model component of force, population and decision-making, is defined by the archetypal Lanchester, Lotka-Volterra and Kuramoto-Sakaguchi models, with feedback between each component adding heterogeneity. Developing a scheme where cultural factors determine decision-making strategies for each force, this work highlights the parametric and topological factors that influence favourable results in a non-linear system where physical outcomes are highly dependent on the non-physical and cognitive nature of each force's intended strategy.

Partial synchronization and community switching in phase-oscillator networks and its analysis based on a bidirectional, weighted chain of three oscillators

Complex networks often possess communities defined based on network connectivity. When dynamics undergo in a network, one can also consider dynamical communities, i.e., a group of nodes displaying a similar dynamical process. We have investigated both analytically and numerically the development of a dynamical community structure, where the community is referred to as a group of nodes synchronized in frequency, in networks of phase oscillators. We first demonstrate that using a few example networks, the community structure changes when network connectivity or interaction strength is varied. In particular, we found that community switching, i.e., a portion of oscillators change the group to which they synchronize, occurs for a range of parameters. We then propose a three-oscillator model: a bidirectional, weighted chain of three Kuramoto phase oscillators, as a theoretical framework for understanding the community formation and its variation. Our analysis demonstrates that the model shows a variety of partially synchronized patterns: oscillators with similar natural frequencies tend to synchronize for weak coupling, while tightly connected oscillators tend to synchronize for strong coupling. We obtain approximate expressions for the critical coupling strengths by employing a perturbative approach in a weak coupling regime and a geometric approach in strong coupling regimes. Moreover, we elucidate the bifurcation types of transitions between different patterns. Our theory might be useful for understanding the development of partially synchronized patterns in a wider class of complex networks than community structured networks.

First-order like phase transition induced by quenched coupling disorder

We investigate the collective dynamics of a population of X Y model-type oscillators, globally coupled via non-separable interactions that are randomly chosen from a positive or negative value and subject to thermal noise controlled by temperature T. We find that the system at T = 0 exhibits a discontinuous, first-order like phase transition from the incoherent to the fully coherent state; when thermal noise is present ( T > 0 ), the transition from incoherence to the partial coherence is continuous and the critical threshold is now larger compared to the deterministic case ( T = 0 ). We derive an exact formula for the critical transition from incoherent to coherent oscillations for the deterministic and stochastic case based on both stability analysis for finite oscillators as well as for the thermodynamic limit ( N → ∞) based on a rigorous mean-field theory using graphons, valid for heterogeneous graph structures. Our theoretical results are supported by extensive numerical simulations. Remarkably, the synchronization threshold induced by the type of random coupling considered here is identical to the one found in studies, which consider uniform input or output strengths for each oscillator node [H. Hong and S. H. Strogatz, Phys. Rev. E 84(4), 046202 (2011); Phys. Rev. Lett. 106(5), 054102 (2011)], which suggests that these systems display a "universal" character for the onset of synchronization.

Swarmalators on a ring with distributed couplings

We study a simple model of identical "swarmalators," generalizations of phase oscillators that swarm through space. We confine the movements to a one-dimensional (1D) ring and consider distributed (nonidentical) couplings; the combination of these two effects captures an aspect of the more realistic two-dimensional swarmalator model. We discover several collective states which we describe analytically. These states imitate the behavior of vinegar eels, catalytic microswimmers, and other swarmalators which move on quasi-1D rings.

Higher-order interactions promote chimera states

Since the discovery of chimera states, the presence of a nonzero phase lag parameter turns out to be an essential attribute for the emergence of chimeras in a nonlocally coupled identical Kuramoto phase oscillators' network with pairwise interactions. In this Letter, we report the emergence of chimeras without phase lag in a nonlocally coupled identical Kuramoto network owing to the introduction of nonpairwise interactions. The influence of added nonlinearity in the coupled system dynamics in the form of simplicial complexes mitigates the requisite of a nonzero phase lag parameter for the emergence of chimera states. Chimera states stimulated by the reciprocity of the pairwise and nonpairwise interaction strengths and their multistable nature are characterized with appropriate measures and are demonstrated in the parameter spaces.

Emerging chimera states under nonidentical counter-rotating oscillators

Frequency plays a crucial role in exhibiting various collective dynamics in the coexisting corotating and counter-rotating systems. To illustrate the impact of counter-rotating frequencies, we consider a network of nonidentical and globally coupled Stuart-Landau oscillators with additional perturbation. Primarily, we investigate the dynamical transitions in the absence of perturbation, demonstrating that the transition from desynchronized state to cluster oscillatory state occurs through an interesting partial synchronization state in the oscillatory regime. Following this, the system dynamics transits to amplitude death and oscillation death states. Importantly, we find that the observed dynamical states do not preserve the parity (P) symmetry in the absence of perturbation. When the perturbation is increased one can note that the system dynamics exhibits a kind of transition which corresponds to a change from incoherent mixed synchronization to coherent mixed synchronization through a chimera state. In particular, incoherent mixed synchronization and coherent mixed synchronization states completely preserve the P symmetry, whereas the chimera state preserves the P symmetry only partially. To demonstrate the occurrence of such partial symmetry-breaking (chimera) state, we use basin stability analysis and discover that partial symmetry breaking exists as a result of the coexistence of symmetry-preserving and symmetry-breaking behavior in the initial state space. Further, a measure of the strength of P symmetry is established to quantify the P symmetry in the observed dynamical states. Subsequently, the dynamical transitions are investigated in the parametric spaces. Finally, by increasing the network size, the robustness of the chimera state is also inspected, and we find that the chimera state is robust even in networks of larger sizes. We also show the generality of the above results in the related reduced phase. model as well as in other coupled models such as the globally coupled van der Pol and Rössler oscillators.

Contrariety and inhibition enhance synchronization in a small-world network of phase oscillators

We numerically study Kuramoto model synchronization consisting of the two groups of conformist-contrarian and excitatory-inhibitory phase oscillators with equal intrinsic frequency. We consider random and small-world (SW) topologies for the connectivity network of the oscillators. In random networks, regardless of the contrarian to conformist connection strength ratio, we found a crossover from the π-state to the blurred π-state and then a continuous transition to the incoherent state by increasing the fraction of contrarians. However, for the excitatory-inhibitory model in a random network, we found that for all the values of the fraction of inhibitors, the two groups remain in phase and the transition point of fully synchronized to an incoherent state reduced by strengthening the ratio of inhibitory to excitatory links. In the SW networks we found that the order parameters for both models do not show monotonic behavior in terms of the fraction of contrarians and inhibitors. Up to the optimal fraction of contrarians and inhibitors, the synchronization rises by introducing the number of contrarians and inhibitors and then falls. We discuss that the nonmonotonic behavior in synchronization is due to the weakening of the defects already formed in the pure conformist and excitatory agent model in SW networks. We found that in SW networks, the optimal fraction of contrarians and inhibitors remain unchanged for the rewiring probabilities up to ∼0.15, above which synchronization falls monotonically, like the random network. We also showed that in the conformist-contrarian model, the optimal fraction of contrarians is independent of the strength of contrarian links. However, in the excitatory-inhibitory model, the optimal fraction of inhibitors is approximately proportional to the inverse of inhibition strength.

Contrarians Synchronize beyond the Limit of Pairwise Interactions

We give evidence that a population of pure contrarian globally coupled D-dimensional Kuramoto oscillators reaches a collective synchronous state when the interplay between the units goes beyond the limit of pairwise interactions. Namely, we will show that the presence of higher-order interactions may induce the appearance of a coherent state even when the oscillators are coupled negatively to the mean field. An exact solution for the description of the microscopic dynamics for forward and backward transitions is provided, which entails imperfect symmetry breaking of the population into a frequency-locked state featuring two clusters of different instantaneous phases. Our results contribute to a better understanding of the powerful potential of group interactions entailing multidimensional choices and novel dynamical states in many circumstances, such as in social systems.

Catalytic feed-forward explosive synchronization in multilayer networks

Inhibitory couplings are crucial for the normal functioning of many real-world complex systems. Inhibition in one layer has been shown to induce explosive synchronization in another excitatory (or positive) layer of duplex networks. By extending this framework to multiplex networks, this article shows that inhibition in a single layer can act as a catalyst, leading to explosive synchronization transitions in the rest of the layers feed-forwarded through intermediate layer(s). Considering a multiplex network of coupled Kuramoto oscillators, we demonstrate that the characteristics of the transition emergent in a layer can be entirely controlled by the intra-layer coupling of other layers and the multiplexing strengths. The results presented here are essential to fathom the synchronization behavior of coupled dynamical units in multi-layer systems possessing inhibitory coupling in one of its layers, representing the importance of multiplexing.

Binary mixtures of locally coupled mobile oscillators

Synchronized behavior in a system of coupled dynamic objects is a fascinating example of an emerged cooperative phenomena which has been observed in systems as diverse as a group of insects, neural networks, or networks of computers. In many instances, however, the synchronization is undesired because it may lead to system malfunctioning, as in the case of Alzheimer's and Parkinson's diseases, for example. Recent studies of static networks of oscillators have shown that the presence of a small fraction of so-called contrarian oscillators can suppress the undesired network synchronization. On the other hand, it is also known that the mobility of the oscillators can significantly impact their synchronization dynamics. Here, we combine these two ideas-the oscillator mobility and the presence of heterogeneous interactions-and study numerically binary mixtures of phase oscillators performing two-dimensional random walks. Within the framework of a generalized Kuramoto model, we introduce two phase-coupling schemes. The first one is invariant when the types of any two oscillators are swapped, while the second model is not. We demonstrate that the symmetric model does not allow for a complete suppression of the synchronized state. However, it provides means for a robust control of the synchronization timescale by varying the overall number density and the composition of the mixture and the strength of the off-diagonal Kuramoto coupling constant. Instead, the asymmetric model predicts that the coherent state can be eliminated within a subpopulation of normal oscillators and evoked within a subpopulation of the contrarians.

Synchronization of coupled Kuramoto oscillators under resource constraints

A fundamental understanding of synchronized behavior in multiagent systems can be acquired by studying analytically tractable Kuramoto models. However, such models typically diverge from many real systems whose dynamics evolve under nonnegligible resource constraints. Here we construct a system of coupled Kuramoto oscillators that consume or produce resources as a function of their oscillation frequency. At high coupling, we observe strongly synchronized dynamics, whereas at low coupling, we observe independent oscillator dynamics as expected from standard Kuramoto models. For intermediate coupling, which typically induces a partially synchronized state, we empirically observe that (and theoretically explain why) the system can exist in either: (i) a state in which the order parameter oscillates in time, or (ii) a state in which multiple synchronization states are simultaneously stable. Whether (i) or (ii) occurs depends upon whether the oscillators consume or produce resources, respectively. Relevant for systems as varied as coupled neurons and social groups, our paper lays important groundwork for future efforts to develop quantitative predictions of synchronized dynamics for systems embedded in environments marked by sparse resources.

Transient circulant clusters in two-population network of Kuramoto oscillators with different rules of coupling adaptation

We considered a network consisting of two populations of phase oscillators, the interaction of which is determined by different rules for the coupling adaptation. The introduction of various adaptation rules leads to the suppression of splay states and the emergence of each population complex non-stationary behavior called transient circulant clusters. In such states, each population contains a pair of anti-phase clusters whose size and composition slowly change over time as a result of successive transitions of oscillators between clusters. We show that an increase in the mismatch of the adaptation rules makes it possible to stop the process of rearrangement of clusters in one or both populations of the network. Transitions to such modes are always preceded by the appearance of solitary states in one of the populations.

Discordant synchronization patterns on directed networks of identical phase oscillators with attractive and repulsive couplings

We study the collective dynamics of identical phase oscillators on globally coupled networks whose interactions are asymmetric and mediated by positive and negative couplings. We split the set of oscillators into two interconnected subpopulations. In this setup, oscillators belonging to the same group interact via symmetric couplings while the interaction between subpopulations occurs in an asymmetric fashion. By employing the dimensional reduction scheme of the Ott-Antonsen (OA) theory, we verify the existence of traveling wave and π-states, in addition to the classical fully synchronized and incoherent states. Bistability between all collective states is reported. Analytical results are generally in excellent agreement with simulations; for some parameters and initial conditions, however, we numerically detect chimera-like states which are not captured by the OA theory.

Non-reciprocal phase transitions

Out of equilibrium, a lack of reciprocity is the rule rather than the exception. Non-reciprocity occurs, for instance, in active matter^1-6, non-equilibrium systems^7-9, networks of neurons^10,11, social groups with conformist and contrarian members¹², directional interface growth phenomena^13-15 and metamaterials^16-20. Although wave propagation in non-reciprocal media has recently been closely studied^1,16-20, less is known about the consequences of non-reciprocity on the collective behaviour of many-body systems. Here we show that non-reciprocity leads to time-dependent phases in which spontaneously broken continuous symmetries are dynamically restored. We illustrate this mechanism with simple robotic demonstrations. The resulting phase transitions are controlled by spectral singularities called exceptional points²¹. We describe the emergence of these phases using insights from bifurcation theory^22,23 and non-Hermitian quantum mechanics^24,25. Our approach captures non-reciprocal generalizations of three archetypal classes of self-organization out of equilibrium: synchronization, flocking and pattern formation. Collective phenomena in these systems range from active time-(quasi)crystals to exceptional-point-enforced pattern formation and hysteresis. Our work lays the foundation for a general theory of critical phenomena in systems whose dynamics is not governed by an optimization principle.

Two-community noisy Kuramoto model with general interaction strengths. II

We generalize the study of the noisy Kuramoto model, considered on a network of two interacting communities, to the case where the interaction strengths within and across communities are taken to be different in general. Using a geometric interpretation of the self-consistency equations developed in Paper I of this series as well as perturbation arguments, we are able to identify all solution boundaries in the phase diagram. This allows us to completely classify the phase diagram in the four-dimensional parameter space and identify all possible bifurcation points. Furthermore, we analyze the asymptotic behavior of the solution boundaries. To illustrate these results and the rich behavior of the model, we present phase diagrams for selected regions of the parameter space.

Two-community noisy Kuramoto model with general interaction strengths. I

We generalize the study of the noisy Kuramoto model, considered on a network of two interacting communities, to the case where the interaction strengths within and across communities are taken to be different in general. By developing a geometric interpretation of the self-consistency equations, we are able to separate the parameter space into ten regions in which we identify the maximum number of solutions in the steady state. Furthermore, we prove that in the steady state, only the angles 0 and π are possible between the average phases of the two communities and derive the solution boundary for the unsynchronized solution. Last, we identify the equivalence class relation in the parameter space corresponding to the symmetrically synchronized solution.

Dynamics in the Sakaguchi-Kuramoto model with bimodal frequency distribution

In this work, we study the Sakaguchi-Kuramoto model with natural frequency following a bimodal distribution. By using Ott-Antonsen ansatz, we reduce the globally coupled phase oscillators to low dimensional coupled ordinary differential equations. For symmetrical bimodal frequency distribution, we analyze the stabilities of the incoherent state and different partial synchronous states. Different types of bifurcations are identified and the effect of the phase lag on the dynamics is investigated. For asymmetrical bimodal frequency distribution, we observe the revival of the incoherent state, and then the conditions for the revival are specified.

Spatial correlations in a finite-range Kuramoto model

We study spatial correlations between oscillator phases in the steady state of a Kuramoto model, in which phase oscillators that are randomly distributed in space interact with constant strength but within a limited range. Such a model could be relevant, for example, in the synchronization of gene expression oscillations in cells, where only oscillations of neighboring cells are coupled through cell-cell contacts. We analytically infer spatial phase-phase correlation functions from the known steady-state distribution of oscillators for the case of homogenous frequencies and show that these can contain information about the range and strength of interactions, provided the noise in the system can be estimated. We suggest a method for the latter, and also explore when correlations appear to be ergodic in this model, which would enable an experimental measurement of correlation functions to utilize temporal averages. Simulations show that our techniques also give qualitative results for the model with heterogenous frequencies. We illustrate our results by comparison with experimental data on genetic oscillations in the segmentation clock of zebrafish embryos.

Understanding the dynamics of biological and neural oscillator networks through exact mean-field reductions: a review

Many biological and neural systems can be seen as networks of interacting periodic processes. Importantly, their functionality, i.e., whether these networks can perform their function or not, depends on the emerging collective dynamics of the network. Synchrony of oscillations is one of the most prominent examples of such collective behavior and has been associated both with function and dysfunction. Understanding how network structure and interactions, as well as the microscopic properties of individual units, shape the emerging collective dynamics is critical to find factors that lead to malfunction. However, many biological systems such as the brain consist of a large number of dynamical units. Hence, their analysis has either relied on simplified heuristic models on a coarse scale, or the analysis comes at a huge computational cost. Here we review recently introduced approaches, known as the Ott-Antonsen and Watanabe-Strogatz reductions, allowing one to simplify the analysis by bridging small and large scales. Thus, reduced model equations are obtained that exactly describe the collective dynamics for each subpopulation in the oscillator network via few collective variables only. The resulting equations are next-generation models: Rather than being heuristic, they exactly link microscopic and macroscopic descriptions and therefore accurately capture microscopic properties of the underlying system. At the same time, they are sufficiently simple to analyze without great computational effort. In the last decade, these reduction methods have become instrumental in understanding how network structure and interactions shape the collective dynamics and the emergence of synchrony. We review this progress based on concrete examples and outline possible limitations. Finally, we discuss how linking the reduced models with experimental data can guide the way towards the development of new treatment approaches, for example, for neurological disease.

The role of timescale separation in oscillatory ensembles with competitive coupling

We study a heterogeneous population consisting of two groups of oscillatory elements, one with attractive and one with repulsive coupling. Moreover, we set different internal timescales for the oscillators of the two groups and concentrate on the role of this timescale separation in the collective behavior. Our results demonstrate that it may significantly modify synchronization properties of the system, and the implications are fundamentally different depending on the ratio between the group timescales. For the slower attractive group, synchronization properties are similar to the case of equal timescales. However, when the attractive group is faster, these properties significantly change and bistability appears. The other collective regimes such as frozen states and solitary states are also shown to be crucially influenced by timescale separation.

Characterizing nonstationary coherent states in globally coupled conformist and contrarian oscillators

For decades, the description and characterization of nonstationary coherent states in coupled oscillators have not been available. We here consider the Kuramoto model consisting of conformist and contrarian oscillators. In the model, contrarians are chosen from a bimodal Lorentzian frequency distribution and flipped into conformists at random. A complete and systematic analytical treatment of the model is provided based on the Ott-Antonsen ansatz. In particular, we predict and analyze not only the stability of all stationary states (such as the incoherent, the π, and the traveling-wave states), but also that of the two nonstationary states: the Bellerophon and the oscillating-π state. The theoretical predictions are fully supported by extensive numerical simulations.

Nonlinear phase coupling functions: a numerical study

Phase reduction is a general tool widely used to describe forced and interacting self-sustained oscillators. Here, we explore the phase coupling functions beyond the usual first-order approximation in the strength of the force. Taking the periodically forced Stuart-Landau oscillator as the paradigmatic model, we determine and numerically analyse the coupling functions up to the fourth order in the force strength. We show that the found nonlinear phase coupling functions can be used for predicting synchronization regions of the forced oscillator. This article is part of the theme issue 'Coupling functions: dynamical interaction mechanisms in the physical, biological and social sciences'.

Solitary states and partial synchrony in oscillatory ensembles with attractive and repulsive interactions

We numerically and analytically analyze transitions between different synchronous states in a network of globally coupled phase oscillators with attractive and repulsive interactions. The elements within the attractive or repulsive group are identical, but natural frequencies of the groups differ. In addition to a synchronous two-cluster state, the system exhibits a solitary state, when a single oscillator leaves the cluster of repulsive elements, as well as partially synchronous quasiperiodic dynamics. We demonstrate how the transitions between these states occur when the repulsion starts to prevail over attraction.

Chimera state on a spherical surface of nonlocally coupled oscillators with heterogeneous phase lags

We consider a network of coupled oscillators embedded in the surface of a sphere with nonlocal coupling strength and heterogeneous phase lags. A nonlocal coupling scheme with heterogeneous phase lags that allows the system to be solved analytically is suggested and the main effects of heterogeneity in the phase lags on the existence and stability of steady states are analyzed. We explore the stability of solutions along the Ott-Antonsen invariant manifold and present a complete bifurcation diagram for stationary patterns including the coherent, incoherent, and modulated drift states as well as chimera state. The stability analysis shows that a continuum of uniform drift states and the modulated drift state could become stable only due to the heterogeneity of the phase lags and that the chimera state is bifurcated from the modulated drift state. Our theoretical results are verified by using the direct numerical simulations of the model system.

Dynamical effects of breaking rotational symmetry in counter-rotating Stuart-Landau oscillators

Stuart-Landau oscillators can be coupled so as to either preserve or destroy the rotational symmetry that the uncoupled system possesses. We examine some of the simplest cases of such couplings for a system of two nonidentical oscillators. When the coupling breaks the rotational invariance, there is a qualitative difference between oscillators wherein the phase velocity has the same sign (termed co-rotation) or opposite signs (termed counter-rotation). In the regime of oscillation death the relative sense of the phase rotations plays a major role. In particular, when rotational invariance is broken, counter-rotation or phase velocities of opposite signs appear to destabilize existing fixed points, thereby preserving and possibly extending the range of oscillatory behavior. The dynamical "frustration" induced by counter-rotations can thus suppress oscillation quenching when coupling breaks the symmetry.

Rhythmic synchronization and hybrid collective states of globally coupled oscillators

Macroscopic rhythms are often signatures of healthy functioning in living organisms, but they are still poorly understood on their microscopic bases. Globally interacting oscillators with heterogeneous couplings are here considered. Thorough theoretical and numerical analyses indicate the presence of multiple phase transitions between different collective states, with regions of bi-stability. Novel coherent phases are unveiled, and evidence is given of the spontaneous emergence of macroscopic rhythms where oscillators’ phases are always found to be self-organized as in Bellerophon states, i.e. in multiple clusters with quantized values of their average frequencies. Due to their rather unconditional appearance, the circumstance is paved that the Bellerophon states grasp the microscopic essentials behind collective rhythms in more general systems of interacting oscillators.

Volcano Transition in a Solvable Model of Frustrated Oscillators

In 1992, a puzzling transition was discovered in simulations of randomly coupled limit-cycle oscillators. This so-called volcano transition has resisted analysis ever since. It was originally conjectured to mark the emergence of an oscillator glass, but here we show it need not. We introduce and solve a simpler model with a qualitatively identical volcano transition and find that its supercritical state is not glassy. We discuss the implications for the original model and suggest experimental systems in which a volcano transition and oscillator glass may appear.

Scale-freeness or partial synchronization in neural mass phase oscillator networks: Pick one of two?

Modeling and interpreting (partial) synchronous neural activity can be a challenge. We illustrate this by deriving the phase dynamics of two seminal neural mass models: the Wilson-Cowan firing rate model and the voltage-based Freeman model. We established that the phase dynamics of these models differed qualitatively due to an attractive coupling in the first and a repulsive coupling in the latter. Using empirical structural connectivity matrices, we determined that the two dynamics cover the functional connectivity observed in resting state activity. We further searched for two pivotal dynamical features that have been reported in many experimental studies: (1) a partial phase synchrony with a possibility of a transition towards either a desynchronized or a (fully) synchronized state; (2) long-term autocorrelations indicative of a scale-free temporal dynamics of phase synchronization. Only the Freeman phase model exhibited scale-free behavior. Its repulsive coupling, however, let the individual phases disperse and did not allow for a transition into a synchronized state. The Wilson-Cowan phase model, by contrast, could switch into a (partially) synchronized state, but it did not generate long-term correlations although being located close to the onset of synchronization, i.e. in its critical regime. That is, the phase-reduced models can display one of the two dynamical features, but not both.

The Dynamics of Networks of Identical Theta Neurons

We consider finite and infinite all-to-all coupled networks of identical theta neurons. Two types of synaptic interactions are investigated: instantaneous and delayed (via first-order synaptic processing). Extensive use is made of the Watanabe/Strogatz (WS) ansatz for reducing the dimension of networks of identical sinusoidally-coupled oscillators. As well as the degeneracy associated with the constants of motion of the WS ansatz, we also find continuous families of solutions for instantaneously coupled neurons, resulting from the reversibility of the reduced model and the form of the synaptic input. We also investigate a number of similar related models. We conclude that the dynamics of networks of all-to-all coupled identical neurons can be surprisingly complicated.

Chimera and modulated drift states in a ring of nonlocally coupled oscillators with heterogeneous phase lags

We consider a ring of phase oscillators with nonlocal coupling strength and heterogeneous phase lags. We analyze the effects of heterogeneity in the phase lags on the existence and stability of a variety of steady states. A nonlocal coupling with heterogeneous phase lags that allows the system to be solved analytically is suggested and the stability of solutions along the Ott-Antonsen invariant manifold is explored. We present a complete bifurcation diagram for stationary patterns including the uniform drift and modulated drift states as well as chimera state, which reveals that the stable modulated drift state and a continuum of metastable drift states could occur due to the heterogeneity of the phase lags. We verify our theoretical results using the direct numerical simulations of the model system.

Long-period clocks from short-period oscillators

We analyze repulsively coupled Kuramoto oscillators, which are exposed to a distribution of natural frequencies. This source of disorder leads to closed orbits of repetitive temporary patterns of phase-locked motion, providing clocks on macroscopic time scales. The periods can be orders of magnitude longer than the periods of individual oscillators. By construction, the attractor space is quite rich. This may cause long transients until the deterministic trajectories find their stationary orbits. The smaller the width of the distribution about the common natural frequency, the longer are the emerging time scales on average. Among the long-period orbits, we find self-similar sequences of temporary phase-locked motion on different time scales. The ratio of time scales is determined by the ratio of widths of the distributions. The results illustrate a mechanism for how simple systems can provide rather flexible dynamics, with a variety of periods even without external entrainment.

Spontaneous symmetry breaking due to the trade-off between attractive and repulsive couplings

Spontaneous symmetry breaking is an important phenomenon observed in various fields including physics and biology. In this connection, we here show that the trade-off between attractive and repulsive couplings can induce spontaneous symmetry breaking in a homogeneous system of coupled oscillators. With a simple model of a system of two coupled Stuart-Landau oscillators, we demonstrate how the tendency of attractive coupling in inducing in-phase synchronized (IPS) oscillations and the tendency of repulsive coupling in inducing out-of-phase synchronized oscillations compete with each other and give rise to symmetry breaking oscillatory states and interesting multistabilities. Further, we provide explicit expressions for synchronized and antisynchronized oscillatory states as well as the so called oscillation death (OD) state and study their stability. If the Hopf bifurcation parameter (λ) is greater than the natural frequency (ω) of the system, the attractive coupling favors the emergence of an antisymmetric OD state via a Hopf bifurcation whereas the repulsive coupling favors the emergence of a similar state through a saddle-node bifurcation. We show that an increase in the repulsive coupling not only destabilizes the IPS state but also facilitates the reentrance of the IPS state.

Multistable states in a system of coupled phase oscillators with inertia

We investigate the generalized Kuramoto model of globally coupled oscillators with inertia, in which oscillators with positive coupling strength are conformists and oscillators with negative coupling strength are contrarians. We consider the correlation between the coupling strengths of oscillators and the distributions of natural frequencies. Two different types of correlations are studied. It is shown that the model supports multistable synchronized states such as different types of travelling wave states, π state and another type of nonstationary state: an oscillating π state. The phase distribution oscillates in a confined region and the phase difference between conformists and contrarians oscillates around π periodically in the oscillating π state. The different types of travelling wave state may be characterized by the speed of travelling wave and the effective frequencies of oscillators. Finally, the bifurcation diagrams of the model in the parameter space are presented.

Intermittent Bellerophon state in frequency-weighted Kuramoto model

Recently, the Bellerophon state, which is a quantized, time dependent, clustering state, was revealed in globally coupled oscillators [Bi et al., Phys. Rev. Lett. 117, 204101 (2016)]. The most important characteristic is that in such a state, the oscillators split into multiple clusters. Within each cluster, the instantaneous frequencies of the oscillators are not the same, but their average frequencies lock to a constant. In this work, we further characterize an intermittent Bellerophon state in the frequency-weighted Kuramoto model with a biased Lorentzian frequency distribution. It is shown that the evolution of oscillators exhibits periodical intermittency, following a synchronous pattern of bursting in a short period and resting in a long period. This result suggests that the Bellerophon state might be generic in Kuramoto-like models regardless of different arrangements of natural frequencies.

Synchronization and Bellerophon states in conformist and contrarian oscillators

The study of synchronization in generalized Kuramoto models has witnessed an intense boost in the last decade. Several collective states were discovered, such as partially synchronized, chimera, π or traveling wave states. We here consider two populations of globally coupled conformist and contrarian oscillators (with different, randomly distributed frequencies), and explore the effects of a frequency-dependent distribution of the couplings on the collective behaviour of the system. By means of linear stability analysis and mean-field theory, a series of exact solutions is extracted describing the critical points for synchronization, as well as all the emerging stationary coherent states. In particular, a novel non-stationary state, here named as Bellerophon state, is identified which is essentially different from all other coherent states previously reported in the Literature. A robust verification of the rigorous predictions is supported by extensive numerical simulations.

Excitation and suppression of chimera states by multiplexing

We study excitation and suppression of chimera states in an ensemble of nonlocally coupled oscillators arranged in a framework of multiplex network. We consider the homogeneous network (all identical oscillators) with different parametric cases and interlayer heterogeneity by introducing parameter mismatch between the layers. We show the feasibility to suppress chimera states in the multiplex network via moderate interlayer interaction between a layer exhibiting chimera state and other layers which are in a coherent or incoherent state. On the contrary, for larger interlayer coupling, we observe the emergence of identical chimera states in both layers which we call an interlayer chimera state. We map the spatiotemporal behavior in a wide range of parameters, varying interlayer coupling strength and phase lag in two and three multiplexing layers. We also prove the emergence of interlayer chimera states in a multiplex network via evaluation of a continuous model. Furthermore, we consider the two-layered network of Hindmarsh-Rose neurons and reveal that in such a system multiplex interaction between layers is capable of exciting not only the synchronous interlayer chimera state but also nonidentical chimera patterns.

Incoherent chimera and glassy states in coupled oscillators with frustrated interactions

We suggest a site disorder model that describes the population of identical oscillators with quenched random interactions for both the coupling strength and coupling phase. We obtain the reduced equations for the suborder parameters, on the basis of Ott-Antonsen ansatz theory, and present a complete bifurcation analysis of the reduced system. New effects include the appearance of the incoherent chimera and glassy state, both of which are caused by heterogeneity of the coupling phases. In the incoherent chimera state, the system displays an exotic symmetry-breaking behavior in spite of the apparent structural symmetry where the oscillators for both of the two subpopulations are in a frustrated state, while the phase distribution for each subpopulation approaches a steady state that differs from each other. When the incoherent chimera undergoes Hopf bifurcation, the system displays a breathing incoherent chimera. The glassy state that occurs on a surface of three-dimensional parameter space exhibits a continuum of metastable states with zero value of the global order parameter. Explicit formulas are derived for the system's Hopf, saddle-node, and transcritical bifurcation curves, as well as the codimension-2 crossing points, including the Takens-Bogdanov point.