Conformists and contrarians in a Kuramoto model with identical natural frequencies

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.

Onset of time dependence in ensembles of excitable elements with global repulsive coupling

We consider the effect of global repulsive coupling on an ensemble of identical excitable elements. An increase of the coupling strength destabilizes the synchronous equilibrium and replaces it with many attracting oscillatory states, created in the transcritical heteroclinic bifurcation. The period of oscillations is inversely proportional to the distance from the critical parameter value. If the elements interact with the global field via the first Fourier harmonics of their phases, the stable equilibrium is in one step replaced by the attracting continuum of periodic motions.

Landau damping effects in the synchronization of conformist and contrarian oscillators

Two decades ago, a phenomenon resembling Landau damping was described in the synchronization of globally coupled oscillators: the evidence of a regime where the order parameter decays when linear theory predicts neutral stability for the incoherent state. We here show that such an effect is far more generic, as soon as phase oscillators couple to their mean field according to their natural frequencies, being then grouped into two distinct populations of conformists and contrarians. We report the analytical solution of this latter situation, which allows determining the critical coupling strength and the stability of the incoherent state, together with extensive numerical simulations that fully support all theoretical predictions. The relevance of our results is discussed in relationship to collective phenomena occurring in social and economical systems.

From Neural and Social Cooperation to the Global Emergence of Cognition

The recent article (Turalska et al., 2012) discusses the emergence of intelligence via criticality as a consequence of locality breakdown. Herein, we use criticality for the foundation of a novel generation of game theory making the local interaction between players yield long-range effects. We first establish that criticality is not confined to the Ising-like structure of the sociological model of (Turalska et al., 2012), called the decision making model (DMM), through the study of the emergence of altruism using the altruism-selfishness model (ASM). Both models generate criticality, one by imitation of opinion (DMM) and the other by imitation of behavior (ASM). The dynamics of a sociological network 𝒮 influences the behavioral network ℱ through two game theoretic paradigms: (i) the value of altruism; (ii) the benefit of rapid consensus. In (i), the network 𝒮 debates the moral issue of altruism by means of the DMM, while at the level ℱ the individuals operate according to the ASM. The individuals of the level 𝒮, through a weak influence on the individuals of the level ℱ, exert a societal control on ℱ, fitting the principle of complexity management and complexity matching. In (ii), the benefit to society is the rapid attainment of consensus in the 𝒮 level. The agents of the level ℱ operate according to the prisoner's dilemma prescription, with the defectors acting as DMM contrarians at the level 𝒮. The contrarians, acting as the inhibitory links of neural networks, exert on society the same beneficial effect of maintaining the criticality-induced resilience that they generate in neural networks. The conflict between personal and social benefit makes the networks evolve toward criticality. Finally, we show that the theory of this article is compatible with recent discoveries in the burgeoning field of social neuroscience.

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.

Mouse hair cycle expression dynamics modeled as coupled mesenchymal and epithelial oscillators

The hair cycle is a dynamic process where follicles repeatedly move through phases of growth, retraction, and relative quiescence. This process is an example of temporal and spatial biological complexity. Understanding of the hair cycle and its regulation would shed light on many other complex systems relevant to biological and medical research. Currently, a systematic characterization of gene expression and summarization within the context of a mathematical model is not yet available. Given the cyclic nature of the hair cycle, we felt it was important to consider a subset of genes with periodic expression. To this end, we combined several mathematical approaches with high-throughput, whole mouse skin, mRNA expression data to characterize aspects of the dynamics and the possible cell populations corresponding to potentially periodic patterns. In particular two gene clusters, demonstrating properties of out-of-phase synchronized expression, were identified. A mean field, phase coupled oscillator model was shown to quantitatively recapitulate the synchronization observed in the data. Furthermore, we found only one configuration of positive-negative coupling to be dynamically stable, which provided insight on general features of the regulation. Subsequent bifurcation analysis was able to identify and describe alternate states based on perturbation of system parameters. A 2-population mixture model and cell type enrichment was used to associate the two gene clusters to features of background mesenchymal populations and rapidly expanding follicular epithelial cells. Distinct timing and localization of expression was also shown by RNA and protein imaging for representative genes. Taken together, the evidence suggests that synchronization between expanding epithelial and background mesenchymal cells may be maintained, in part, by inhibitory regulation, and potential mediators of this regulation were identified. Furthermore, the model suggests that impairing this negative regulation will drive a bifurcation which may represent transition into a pathological state such as hair miniaturization.

The hair cycle represents a complex process of particular interest in the study of regulated proliferation, apoptosis and differentiation. While various modeling strategies are presented in the literature, none attempt to link extensive molecular details, provided by high-throughput experiments, with high-level, system properties. Thus, we re-analyzed a previously published mRNA expression time course study and found that we could readily identify a sizeable subset of genes that was expressed in synchrony with the hair cycle itself. The data is summarized in a dynamic, mathematical model of coupled oscillators. We demonstrate that a particular coupling scheme is sufficient to explain the observed synchronization. Further analysis associated specific expression patterns to general yet distinct cell populations, background mesenchymal and rapidly expanding follicular epithelial cells. Experimental imaging results are presented to show the localization of candidate genes from each population. Taken together, the results describe a possible mechanism for regulation between epithelial and mesenchymal populations. We also described an alternate state similar to hair miniaturization, which is predicted by the oscillator model. This study exemplifies the strengths of combining systems-level analysis with high-throughput experimental data to obtain a novel view of a complex system such as the hair cycle.

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.

Dynamics in the Sakaguchi-Kuramoto model with two subpopulations [corrected]

The dynamics in a variant of globally coupled Sakaguchi-Kuramoto [corrected]. phase oscillators is studied. The model consists of two subpopulations, each with a different phase lag and interaction strength. Using Ott-Antonson ansatz, we analyze the dynamics in the model and present the numerical results. There exist stationary synchronous states which are generalized π states and two types of traveling wave states. We find that the traveling wave states are the dominant dynamics in comparison with the stationary states. Particularly, we find that the stationary and traveling wave states can be smoothly connected through the properly chosen parameter paths.

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.

Solitary state at the edge of synchrony in ensembles with attractive and repulsive interactions

We discuss the desynchronization transition in networks of globally coupled identical oscillators with attractive and repulsive interactions. We show that, if attractive and repulsive groups act in antiphase or close to that, a solitary state emerges with a single repulsive oscillator split up from the others fully synchronized. With further increase of the repulsing strength, the synchronized cluster becomes fuzzy and the dynamics is given by a variety of stationary states with zero common forcing. Intriguingly, solitary states represent the natural link between coherence and incoherence. The phenomenon is described analytically for phase oscillators with sine coupling and demonstrated numerically for more general amplitude models.

Physical aging of classical oscillators

Aging is a familiar phenomenon from glassy systems like spin glasses and materials with slow relaxation processes, breaking of time-translation invariance, and dynamical scaling. We study aging in active rotators and Kuramoto oscillators that are coupled with frustrated bonds. The induced multiplicity of attractors of fixed-point or limit-cycle solutions leads to a rough potential landscape. When the system is exposed to noise, the oscillator phases migrate through this landscape and generate a multitude of different escape times from one metastable state to the next. When the system is quenched from the regime of a unique fixed point toward the regime of multistable limit-cycle solutions, the autocorrelation functions depend on the waiting time after the quench and show dynamical scaling. In this way we uncover a common mechanism behind aging in quite different realizations.

Transition to synchronization in a Kuramoto model with the first- and second-order interaction terms

We investigate a Kuramoto model incorporated with the first-order and the second-order interaction terms. We show that the model displays the coexistence of multiattractors and different attractors may be characterized by the phase distributions of oscillators. By investigating the transition diagrams in both forward continuation and backward continuation, we find that the synchronous state with unimodal phase distribution is the most stable one while the state in cluster synchrony with evenly distributed bimodal phase distribution is the least stable one. We also present the phase diagram of the model in the parameter space.

Dynamics of the Kuramoto model in the presence of correlation between distributions of frequencies and coupling strengths

As a paradigmatic model, the Kuramoto model has provided a platform for investigating synchronization among nonidentical oscillators. In this work, we consider the Kuramoto model consisting of conformists with positive coupling strength and contrarians with negative coupling strength. We introduce the correlation between the distributions of natural frequencies and the coupling strengths of oscillators. Three different types of correlations are considered. We find rich dynamics result from the correlation such as different types of traveling wave states and, most interestingly, another type of nonstationary state: an oscillating π state.

Complex dynamics of a nonlinear voter model with contrarian agents

We investigate mean-field dynamics of a nonlinear opinion formation model with congregator and contrarian agents. Each agent assumes one of the two possible states. Congregators imitate the state of other agents with a rate that increases with the number of other agents in the opposite state, as in the linear voter model and nonlinear majority voting models. Contrarians flip the state with a rate that increases with the number of other agents in the same state. The nonlinearity controls the strength of the majority voting and is used as a main bifurcation parameter. We show that the model undergoes a rich bifurcation scenario comprising the egalitarian equilibrium, two symmetric lopsided equilibria, limit cycle, and coexistence of different types of stable equilibria with intertwining attractive basins.

An efficient approach to suppress the negative role of contrarian oscillators in synchronization

It has been found that contrarian oscillators usually take a negative role in the collective behaviors formed by conformist oscillators. However, experiments revealed that it is also possible to achieve a strong coherence even when there are contrarians in the system such as neuron networks with both excitable and inhibitory neurons. To understand the underlying mechanism of this abnormal phenomenon, we here consider a complex network of coupled Kuramoto oscillators with mixed positive and negative couplings and present an efficient approach, i.e., tit-for-tat strategy, to suppress the negative role of contrarian oscillators in synchronization and thus increase the order parameter of synchronization. Two classes of contrarian oscillators are numerically studied and a brief theoretical analysis is provided to explain the numerical results.

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.

How to suppress undesired synchronization

Examples of synchronization can be found in a wide range of phenomena such as neurons firing, lasers cascades, chemical reactions, and opinion formation. However, in many situations the formation of a coherent state is not pleasant and should be mitigated. For example, the onset of synchronization can be the root of epileptic seizures, traffic congestion in networks, and the collapse of constructions. Here we propose the use of contrarians to suppress undesired synchronization. We perform a comparative study of different strategies, either requiring local or total knowledge, and show that the most efficient one solely requires local information. Our results also reveal that, even when the distribution of neighboring interactions is narrow, significant improvement is observed when contrarians sit at the highly connected elements. The same qualitative results are obtained for artificially generated networks and two real ones, namely, the Routers of the Internet and a neuronal network.

Onset of synchronization in complex networks of noisy oscillators

We study networks of noisy phase oscillators whose nodes are characterized by random degrees counting the number of their connections. Both these degrees and the natural frequencies of the oscillators are distributed according to a given probability density. Replacing the randomly connected network by an all-to-all coupled network with weighted edges allows us to formulate the dynamics of a single oscillator coupled to the mean field and to derive the corresponding Fokker-Planck equation. From the latter we calculate the critical coupling strength for the onset of synchronization as a function of the noise intensity, the frequency distribution, and the first two moments of the degree distribution. Our approach is applied to a dense small-world network model, for which we calculate the degree distribution. Numerical simulations prove the validity of the replacement. We also test the applicability to more sparsely connected networks and formulate homogeneity and absence of correlations in the degree distribution as limiting factors of our approach.

Collective synchronization in the presence of reactive coupling and shear diversity

We analyze the synchronization dynamics of a model obtained from the phase reduction of the mean-field complex Ginzburg-Landau equation with heterogeneity. We present exact results that uncover the role of dissipative and reactive couplings on the synchronization transition when shears and natural frequencies are independently distributed. As it occurs in the purely dissipative case, an excess of shear diversity prevents the onset of synchronization, but this does not hold true if coupling is purely reactive. In this case, the synchronization threshold turns out to depend on the mean of the shear distribution, but not on all the other distribution's moments.