Mean-field behavior in coupled oscillators with attractive and repulsive interactions

Realization of logic gates in bi-directionally coupled nonlinear oscillators

Implementation of logic gates has been investigated in nonlinear dynamical systems from various perspectives over the years. Specifically, logic gates have been implemented in both single nonlinear systems and coupled nonlinear oscillators. The majority of the works in the literature have been done on the evolution of single oscillators into OR/AND or NOR/NAND logic gates. In the present study, we demonstrate the design of logic gates in bi-directionally coupled double-well Duffing oscillators by applying two logic inputs to the drive system alone along with a fixed bias. The nonlinear system, comprising both bi-directional components, exhibits varied logic behaviors within an optimal range of coupling strength. Both attractive and repulsive couplings yield similar and complementary logic behaviors in the first and second oscillators. These couplings play a major role in exhibiting fundamental and universal logic gates in simple nonlinear systems. Under a positive bias, both the first and second oscillators demonstrate OR logic gate for the attractive coupling, while exhibiting OR and NOR logic gates, respectively, for the repulsive coupling. Conversely, under a negative bias, both the first and second oscillators display AND logic gate for the attractive coupling, and AND and NAND logical outputs for the repulsive coupling. Furthermore, we confirm the robustness of the bi-directional oscillators against moderate noise in maintaining the desired logical outputs.

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.

Asymmetry in the Kuramoto model with nonidentical coupling

Synchronization and desynchronization of coupled oscillators appear to be the key property of many physical systems. It is believed that to predict a synchronization (or desynchronization) event, the knowledge on the exact structure of the oscillatory network is required. However, natural sciences often deal with observations where the coupling coefficients are not available. In the present paper we suggest a way to characterize synchronization of two oscillators without the reconstruction of coupling. Our method is based on the Kuramoto chain with three oscillators with constant but nonidentical coupling. We characterize coupling in this chain by two parameters: the coupling strength s and disparity σ. We give an analytical expression of the boundary s_{max} of synchronization occurred when s>s_{max}. We propose asymmetry A of the generalized order parameter induced by the coupling disparity as a new characteristic of the synchronization between two oscillators. For the chain model with three oscillators we present the self-consistent inverse problem. We explore scaling properties of the asymmetry A constructed for the inverse problem. We demonstrate that the asymmetry A in the chain model is maximal when the coupling strength in the model reaches the boundary of synchronization s_{max}. We suggest that the asymmetry A may be derived from the phase difference of any two oscillators if one pretends that they are edges of an abstract chain with three oscillators. Performing such a derivation with the general three-oscillator Kuramoto model, we show that the crossover from the chain to general network of oscillators keeps the interrelation between the asymmetry A and synchronization. Finally, we apply the asymmetry A to describe synchronization of the solar magnetic field proxies and discuss its potential use for the forecast of solar cycle anomalies.

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.

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 of multilayer networks with amplification

We study the dynamics of a multilayer network of chaotic oscillators subject to amplification. Previous studies have proven that multilayer networks present phenomena such as synchronization, cluster, and chimera states. Here, we consider a network with two layers of Rössler chaotic oscillators as well as applications to multilayer networks of the chaotic jerk and Liénard oscillators. Intra-layer coupling is considered to be all to all in the case of Rössler oscillators, a ring for jerk oscillators and global mean field coupling in the case of Liénard, inter-layer coupling is unidirectional in all these three cases. The second layer has an amplification coefficient. An in-depth study on the case of a network of Rössler oscillators using a master stability function and order parameter leads to several phenomena such as complete synchronization, generalized, cluster, and phase synchronization with amplification. For the case of Rössler oscillators, we note that there are also certain values of coupling parameters and amplification where the synchronization does not exist or the synchronization can exist but without amplification. Using other systems with different topologies, we obtain some interesting results such as chimera state with amplification, cluster state with amplification, and complete synchronization with amplification.

Nonreciprocity as a generic route to traveling states

We examine a nonreciprocally coupled dynamical model of a mixture of two diffusing species. We demonstrate that nonreciprocity, which is encoded in the model via antagonistic cross-diffusivities, provides a generic mechanism for the emergence of traveling patterns in purely diffusive systems with conservative dynamics. In the absence of nonreciprocity, the binary fluid mixture undergoes a phase transition from a homogeneous mixed state to a demixed state with spatially separated regions rich in one of the two components. Above a critical value of the parameter tuning nonreciprocity, the static demixed pattern acquires a finite velocity, resulting in a state that breaks both spatial and time-reversal symmetry, as well as the reflection parity of the static pattern. We elucidate the generic nature of the transition to traveling patterns using a minimal model that can be studied analytically. Our work has direct relevance to nonequilibrium assembly in mixtures of chemically interacting colloids that are known to exhibit nonreciprocal effective interactions, as well as to mixtures of active and passive agents where traveling states of the type predicted here have been observed in simulations. It also provides insight on transitions to traveling and oscillatory states seen in a broad range of nonreciprocal systems with nonconservative dynamics, from reaction-diffusion and prey-predators models to multispecies mixtures of microorganisms with antagonistic interactions.

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.

Static and dynamic attractive-repulsive interactions in two coupled nonlinear oscillators

Many systems exhibit both attractive and repulsive types of interactions, which may be dynamic or static. A detailed understanding of the dynamical properties of a system under the influence of dynamically switching attractive or repulsive interactions is of practical significance. However, it can also be effectively modeled with two coexisting competing interactions. In this work, we investigate the effect of time-varying attractive-repulsive interactions as well as the hybrid model of coexisting attractive-repulsive interactions in two coupled nonlinear oscillators. The dynamics of two coupled nonlinear oscillators, specifically limit cycles as well as chaotic oscillators, are studied in detail for various dynamical transitions for both cases. Here, we show that dynamic or static attractive-repulsive interactions can induce an important transition from the oscillatory to steady state in identical nonlinear oscillators due to competitive effects. The analytical condition for the stable steady state in dynamic interactions at the low switching time period and static coexisting interactions are calculated using linear stability analysis, which is found to be in good agreement with the numerical results. In the case of a high switching time period, oscillations are revived for higher interaction strength.

Diversity of dynamical behaviors due to initial conditions: Extension of the Ott-Antonsen ansatz for identical Kuramoto-Sakaguchi phase oscillators

The Ott-Antonsen ansatz is a powerful tool to extract the behaviors of coupled phase oscillators, but it imposes a strong restriction on the initial condition. Herein, an extension of the Ott-Antonsen ansatz is proposed to relax the restriction, enabling the systematic approximation of the behavior of a globally coupled phase oscillator system with an arbitrary initial condition. The proposed method is validated on the Kuramoto-Sakaguchi model of identical phase oscillators. The method yields cluster and chimera-like solutions that are not obtained by the conventional ansatz.

Additional repulsion reduces the dynamical resilience in the damaged networks

In this paper, we investigate the dynamical robustness of diffusively coupled oscillatory networks under the influence of an additional repulsive link. Such a dynamical resilience property is realized through the aging process of the damaged network of active and inactive oscillators. The aging process is one type of phase transition, mainly appearing at a critical threshold of a fraction of the inactive oscillator node where the mean oscillation amplitude of the entire network suddenly vanishes. These critical fractions of the failure nodes in the network are broadly used as a measure of network resilience. Here, we analytically derived the critical fraction of the aging process in the dynamical network. We find that the addition of the repulsive link enhances the critical threshold of the aging transition of diffusively coupled oscillators, which indicated that the dynamical robustness of the coupled network decreases with the presence of the repulsive interaction. Furthermore, we investigate the dynamical robustness of the network against the number of deteriorating repulsive links. We observed that a certain percentage of the repulsive link is enabled to produce the aging process in the entire network. Finally, the effect of symmetry-breaking coupling and the targeted inactivation process on the dynamical robustness property of damaged networks were investigated. The analytically obtained results are verified numerically in the network of coupled Stuart-Landau oscillators. These findings may help us to better understand the role of the coupling mechanism on the phase transition in the damaged network.

A coupled oscillator model for the origin of bimodality and multimodality

Perhaps because of the elegance of the central limit theorem, it is often assumed that distributions in nature will approach singly-peaked, unimodal shapes reminiscent of the Gaussian normal distribution. However, many systems behave differently, with variables following apparently bimodal or multimodal distributions. Here, we argue that multimodality may emerge naturally as a result of repulsive or inhibitory coupling dynamics, and we show rigorously how it emerges for a broad class of coupling functions in variants of the paradigmatic Kuramoto model.

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 ordinary differential equation (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 that all oscillators are still governed by the same ODE, but 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 synchronized, as well as multiple families of neutrally stable states or closed orbits, in which no two oscillators are synchronized. 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.

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.

Kuramoto Model for Excitation-Inhibition-Based Oscillations

The Kuramoto model (KM) is a theoretical paradigm for investigating the emergence of rhythmic activity in large populations of oscillators. A remarkable example of rhythmogenesis is the feedback loop between excitatory (E) and inhibitory (I) cells in large neuronal networks. Yet, although the EI-feedback mechanism plays a central role in the generation of brain oscillations, it remains unexplored whether the KM has enough biological realism to describe it. Here we derive a two-population KM that fully accounts for the onset of EI-based neuronal rhythms and that, as the original KM, is analytically solvable to a large extent. Our results provide a powerful theoretical tool for the analysis of large-scale neuronal oscillations.

Superslow relaxation in identical phase oscillators with random and frustrated interactions

This paper is concerned with the relaxation dynamics of a large population of identical phase oscillators, each of which interacts with all the others through random couplings whose parameters obey the same Gaussian distribution with the average equal to zero and are mutually independent. The results obtained by numerical simulation suggest that for the infinite-size system, the absolute value of Kuramoto's order parameter exhibits superslow relaxation, i.e., 1/ln t as time t increases. Moreover, the statistics on both the transient time T for the system to reach a fixed point and the absolute value of Kuramoto's order parameter at t = T are also presented together with their distribution densities over many realizations of the coupling parameters.

Finite-size scaling in the system of coupled oscillators with heterogeneity in coupling strength

We consider a mean-field model of coupled phase oscillators with random heterogeneity in the coupling strength. The system that we investigate here is a minimal model that contains randomness in diverse values of the coupling strength, and it is found to return to the original Kuramoto model [Y. Kuramoto, Prog. Theor. Phys. Suppl. 79, 223 (1984)10.1143/PTPS.79.223] when the coupling heterogeneity disappears. According to one recent paper [H. Hong, H. Chaté, L.-H. Tang, and H. Park, Phys. Rev. E 92, 022122 (2015)10.1103/PhysRevE.92.022122], when the natural frequency of the oscillator in the system is "deterministically" chosen, with no randomness in it, the system is found to exhibit the finite-size scaling exponent ν[over ¯]=5/4. Also, the critical exponent for the dynamic fluctuation of the order parameter is found to be given by γ=1/4, which is different from the critical exponents for the Kuramoto model with the natural frequencies randomly chosen. Originally, the unusual finite-size scaling behavior of the Kuramoto model was reported by Hong et al. [H. Hong, H. Chaté, H. Park, and L.-H. Tang, Phys. Rev. Lett. 99, 184101 (2007)10.1103/PhysRevLett.99.184101], where the scaling behavior is found to be characterized by the unusual exponent ν[over ¯]=5/2. On the other hand, if the randomness in the natural frequency is removed, it is found that the finite-size scaling behavior is characterized by a different exponent, ν[over ¯]=5/4 [H. Hong, H. Chaté, L.-H. Tang, and H. Park, Phys. Rev. E 92, 022122 (2015)10.1103/PhysRevE.92.022122]. Those findings brought about our curiosity and led us to explore the effects of the randomness on the finite-size scaling behavior. In this paper, we pay particular attention to investigating the finite-size scaling and dynamic fluctuation when the randomness in the coupling strength is considered.

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.

Correlated disorder in the Kuramoto model: Effects on phase coherence, finite-size scaling, and dynamic fluctuations

We consider a mean-field model of coupled phase oscillators with quenched disorder in the natural frequencies and coupling strengths. A fraction p of oscillators are positively coupled, attracting all others, while the remaining fraction 1-p are negatively coupled, repelling all others. The frequencies and couplings are deterministically chosen in a manner which correlates them, thereby correlating the two types of disorder in the model. We first explore the effect of this correlation on the system's phase coherence. We find that there is a critical width γ_c in the frequency distribution below which the system spontaneously synchronizes. Moreover, this γ_c is independent of p. Hence, our model and the traditional Kuramoto model (recovered when p = 1) have the same critical width γ_c. We next explore the critical behavior of the system by examining the finite-size scaling and the dynamic fluctuation of the traditional order parameter. We find that the model belongs to the same universality class as the Kuramoto model with deterministically (not randomly) chosen natural frequencies for the case of p < 1.

Synchronization unveils the organization of ecological networks with positive and negative interactions

Network science has helped to understand the organization principles of the interactions among the constituents of large complex systems. However, recently, the high resolution of the data sets collected has allowed to capture the different types of interactions coexisting within the same system. A particularly important example is that of systems with positive and negative interactions, a usual feature appearing in social, neural, and ecological systems. The interplay of links of opposite sign presents natural difficulties for generalizing typical concepts and tools applied to unsigned networks and, moreover, poses some questions intrinsic to the signed nature of the network, such as how are negative interactions balanced by positive ones so to allow the coexistence and survival of competitors/foes within the same system? Here, we show that synchronization phenomenon is an ideal benchmark for uncovering such balance and, as a byproduct, to assess which nodes play a critical role in the overall organization of the system. We illustrate our findings with the analysis of synthetic and real ecological networks in which facilitation and competitive interactions coexist.

Statistics of Lyapunov exponent spectrum in randomly coupled Kuramoto oscillators

Characterization of spatiotemporal dynamics of coupled oscillatory systems can be done by computing the Lyapunov exponents. We study the spatiotemporal dynamics of randomly coupled network of Kuramoto oscillators and find that the spectral statistics obtained from the Lyapunov exponent spectrum show interesting sensitivity to the coupling matrix. Our results indicate that in the weak coupling limit the gap distribution of the Lyapunov spectrum is Poissonian, while in the limit of strong coupling the gap distribution shows level repulsion. Moreover, the oscillators settle to an inhomogeneous oscillatory state, and it is also possible to infer the random network properties from the Lyapunov exponent spectrum.

Phase coherence induced by correlated disorder

We consider a mean-field model of coupled phase oscillators with quenched disorder in the coupling strengths and natural frequencies. When these two kinds of disorder are uncorrelated (and when the positive and negative couplings are equal in number and strength), it is known that phase coherence cannot occur and synchronization is absent. Here we explore the effects of correlating the disorder. Specifically, we assume that any given oscillator either attracts or repels all the others, and that the sign of the interaction is deterministically correlated with the given oscillator's natural frequency. For symmetrically correlated disorder with zero mean, we find that the system spontaneously synchronizes, once the width of the frequency distribution falls below a critical value. For asymmetrically correlated disorder, the model displays coherent traveling waves: the complex order parameter becomes nonzero and rotates with constant frequency different from the system's mean natural frequency. Thus, in both cases, correlated disorder can trigger phase coherence.

Optimal synchronization of Kuramoto oscillators: A dimensional reduction approach

A recently proposed dimensional reduction approach for studying synchronization in the Kuramoto model is employed to build optimal network topologies to favor or to suppress synchronization. The approach is based in the introduction of a collective coordinate for the time evolution of the phase locked oscillators, in the spirit of the Ott-Antonsen ansatz. We show that the optimal synchronization of a Kuramoto network demands the maximization of the quadratic function ω(T)Lω, where ω stands for the vector of the natural frequencies of the oscillators and L for the network Laplacian matrix. Many recently obtained numerical results can be reobtained analytically and in a simpler way from our maximization condition. A computationally efficient hill climb rewiring algorithm is proposed to generate networks with optimal synchronization properties. Our approach can be easily adapted to the case of the Kuramoto models with both attractive and repulsive interactions, and again many recent numerical results can be rederived in a simpler and clearer analytical manner.

Collective dynamics in two populations of noisy oscillators with asymmetric interactions

We study two intertwined globally coupled networks of noisy Kuramoto phase oscillators that have the same natural frequency but differ in their perception of the mean field and their contribution to it. Such a give-and-take mechanism is given by asymmetric in- and out-coupling strengths which can be both positive and negative. We uncover in this minimal network of networks intriguing patterns of discordance, where the ensemble splits into two clusters separated by a constant phase lag. If it differs from π, then traveling wave solutions emerge. We observe a second route to traveling waves via traditional one-cluster states. Bistability is found between the various collective states. Analytical results and bifurcation diagrams are derived with a reduced system.

Traveling wave in a three-dimensional array of conformist and contrarian oscillators

We consider a system of conformist and contrarian oscillators coupled locally in a three-dimensional cubic lattice and explore collective behavior of the system. The conformist oscillators attractively interact with the neighbor oscillators and therefore tend to be aligned with the neighbors' phase. The contrarian oscillators interact repulsively with the neighbors and therefore tend to be out of phase with them. In this paper, we investigate whether many peculiar dynamics that have been observed in the mean-field system with global coupling can emerge even with local coupling. In particular, we pay attention to the possibility that a traveling wave may arise. We find that the traveling wave occurs due to coupling asymmetry and not by global coupling; this observation confirms that the global coupling is not essential to the occurrence of a traveling wave in the system. The traveling wave can be a mechanism for the coherent rhythm generation of the circadian clock or of hormone secretion in biological systems under local coupling.

Susceptibility of large populations of coupled oscillators

It is an important and interesting problem to elucidate how the degree of phase order in a large population of coupled oscillators responds to a synchronizing periodic force from the outside. Here this problem is studied analytically as well as numerically by introducing the concept of susceptibility for globally coupled phase oscillators with either nonrandom or random interactions. It is shown that the susceptibility diverges at the critical point in the nonrandom case with Widom's equality satisfied, while it exhibits a cusp in the most random case.

Bifurcation study of phase oscillator systems with attractive and repulsive interaction

We study a model of globally coupled phase oscillators that contains two groups of oscillators with positive (synchronizing) and negative (desynchronizing) incoming connections for the first and second groups, respectively. This model was previously studied by Hong and Strogatz (the Hong-Strogatz model) in the case of a large number of oscillators. We consider a generalized Hong-Strogatz model with a constant phase shift in coupling. Our approach is based on the study of invariant manifolds and bifurcation analysis of the system. In the case of zero phase shift, various invariant manifolds are analytically described and a new dynamical mode is found. In the case of a nonzero phase shift we obtained a set of bifurcation diagrams for various systems with three or four oscillators. It is shown that in these cases system dynamics can be complex enough and include multistability and chaotic oscillations.

Glassy states and super-relaxation in populations of coupled phase oscillators

Large networks of coupled oscillators appear in many branches of science, so that the kinds of phenomena they exhibit are not only of intrinsic interest but also of very wide importance. In 1975, Kuramoto proposed an analytically tractable model to describe these systems, which has since been successfully applied in many contexts and remains a subject of intensive research. Some related problems, however, remain unclarified for decades, such as the existence and properties of the oscillator glass state. Here we present a detailed analysis of a very general form of the Kuramoto model. In particular, we find the conditions when it can exhibit glassy behaviour, which represents a kind of synchronous disorder in the present case. Furthermore, we discover a new and intriguing phenomenon that we refer to as super-relaxation where the oscillators feel no interaction at all while relaxing to incoherence. Our findings offer the possibility of creating glassy states and observing super-relaxation in real systems.

Periodic synchronization and chimera in conformist and contrarian oscillators

We consider a system of phase oscillators that couple with both attractive and repulsive interaction under a pinning force and explore collective behavior of the system. The oscillators can be divided into two subpopulations of "conformist" oscillators with attractive interaction and "contrarian" ones with repulsive interaction. We find that the interplay between the pinning force and the opposite relationship of the conformist and contrarian oscillators induce peculiar dynamic states: periodic synchronization, breathing chimera, and fully pinned state depending on the fraction of the conformists. Using the Watanabe-Strogatz transformation, we reduce the dynamics into a low-dimensional one and find that the above dynamic states are generated from the reduced dynamics.

Symmetry and symmetry breaking in a Kuramoto model induced on a Möbius strip

The concept of bounded confidence in social dynamics is introduced into the Kuramoto model. Then a model with projective symmetry is naturally induced by a principal Z(2) bundle over a circle. A Möbius strip is constructed as the manifold of interactions such that two "sides" of a local section correspond to antipodal couplings with opposite signs. We show analytically that symmetric polarization of synchronized states (conformist vs contrarian) can emerge from the collective behavior of homogeneously coupled oscillators. In the continuum limit, a generalized asymmetric model reverts to the symmetric case under uniform initial condition, whereas for networks of finite size, it exhibits richer dynamics or the ordinary behavior of the classic model via symmetry breaking, depending on the relevant parameters values. The symmetry and symmetry breaking of the synchronization state are analogous to the concepts of equality and majority of opinions in sociophysical models of opinion formation. The synchronized clusters can be one (consensus), two (polarization), or more (fragmentation), as opinion clusters.

Fashion, cooperation, and social interactions

Fashion plays such a crucial rule in the evolution of culture and society that it is regarded as a second nature to the human being. Also, its impact on economy is quite nontrivial. On what is fashionable, interestingly, there are two viewpoints that are both extremely widespread but almost opposite: conformists think that what is popular is fashionable, while rebels believe that being different is the essence. Fashion color is fashionable in the first sense, and Lady Gaga in the second. We investigate a model where the population consists of the afore-mentioned two groups of people that are located on social networks (a spatial cellular automata network and small-world networks). This model captures two fundamental kinds of social interactions (coordination and anti-coordination) simultaneously, and also has its own interest to game theory: it is a hybrid model of pure competition and pure cooperation. This is true because when a conformist meets a rebel, they play the zero sum matching pennies game, which is pure competition. When two conformists (rebels) meet, they play the (anti-) coordination game, which is pure cooperation. Simulation shows that simple social interactions greatly promote cooperation: in most cases people can reach an extraordinarily high level of cooperation, through a selfish, myopic, naive, and local interacting dynamic (the best response dynamic). We find that degree of synchronization also plays a critical role, but mostly on the negative side. Four indices, namely cooperation degree, average satisfaction degree, equilibrium ratio and complete ratio, are defined and applied to measure people’s cooperation levels from various angles. Phase transition, as well as emergence of many interesting geographic patterns in the cellular automata network, is also observed.