Basin of attraction determines hysteresis in explosive synchronization

Canard explosions in turbulent thermo-fluid systems

A sudden transition to a state of high-amplitude periodic oscillations is catastrophic in a thermo-fluid system. Conventionally, upon varying the control parameter, a sudden transition is observed as an abrupt jump in the amplitude of the fluctuations in these systems. In contrast, we present an experimental discovery of a canard explosion in a turbulent reactive flow system where we observe a continuous bifurcation with a rapid rise in the amplitude of the fluctuations within a narrow range of control parameters. The observed transition is facilitated via a state of bursting, consisting of the epochs of large amplitude periodic oscillations amidst the epochs of low-amplitude periodic oscillations. The amplitude of the bursts is higher than the amplitude of the bursts of an intermittency state in a conventional gradual transition, as reported in turbulent reactive flow systems. During the bursting state, we observe that temperature fluctuations of the exhaust gas vary at a slower time scale in correlation with the amplitude envelope of the bursts. We also present a phenomenological model for thermoacoustic systems to describe the observed canard explosion. Using the model, we explain that the large amplitude bursts occur due to the slow-fast dynamics at the bifurcation regime of the canard explosion.

Explosive synchronization in a turbulent reactive flow system

The occurrence of abrupt dynamical transitions in the macroscopic state of a system has received growing attention. We present experimental evidence for abrupt transition via explosive synchronization in a real-world complex system, namely, a turbulent reactive flow system. In contrast to the paradigmatic continuous transition to a synchronized state from an initially desynchronized state, the system exhibits a discontinuous synchronization transition with a hysteresis. We consider the fluctuating heat release rate from the turbulent flames at each spatial location as locally coupled oscillators that are coupled to the global acoustic field in the confined system. We analyze the synchronization between these two subsystems during the transition to a state of oscillatory instability and discover that explosive synchronization occurs at the onset of oscillatory instability. Further, we explore the underlying mechanism of interaction between the subsystems and construct a mathematical model of the same.

Exploring the intensity, distribution and evolution of teleconnections using climate network analysis

Teleconnections refer to long-range climate system linkages occurring over typically thousands of kilometers. Generally speaking, most teleconnections are attributed to the transmission of energy and propagation of waves although the physical complexity and characteristics behind these waves are not fully understood. To address this knowledge gap, we develop a climate network-based approach to reveal their directions and distribution patterns, evaluate the intensity of teleconnections, and identify sensitive regions using global daily surface air temperature data. Our results reveal a stable average intensity distribution pattern for teleconnections across a substantial spatiotemporal scale from 1948 to 2021, with the extent and intensity of teleconnection impacts increasing more prominently in the Southern Hemisphere over the past 37 years. Furthermore, we pinpoint climate-sensitive regions, such as southeastern Australia, which are likely to face increasing impacts due to global warming. Our proposed method offers new insights into the dynamics of global climate patterns and can inform strategies to address climate change and extreme events.

Topologically induced suppression of explosive synchronization

Nowadays, explosive synchronization is a well-documented phenomenon consisting in a first-order transition that may coexist with classical synchronization. Typically, explosive synchronization occurs when the network structure is represented by the classical graph Laplacian, and the node frequency and its degree are correlated. Here, we answer the question on whether this phenomenon can be observed in networks when the oscillators are coupled via degree-biased Laplacian operators. We not only observe that this is the case but also that this new representation naturally controls the transition from explosive to standard synchronization in a network. We prove analytically that explosive synchronization emerges when using this theoretical setting in star-like networks. As soon as this star-like network is topologically converted into a network containing cycles, the explosive synchronization gives rise to classical synchronization. Finally, we hypothesize that this mechanism may play a role in switching from normal to explosive states in the brain, where explosive synchronization has been proposed to be related to some pathologies like epilepsy and fibromyalgia.

Local Dirac Synchronization on networks

We propose Local Dirac Synchronization that uses the Dirac operator to capture the dynamics of coupled nodes and link signals on an arbitrary network. In Local Dirac Synchronization, the harmonic modes of the dynamics oscillate freely while the other modes are interacting non-linearly, leading to a collectively synchronized state when the coupling constant of the model is increased. Local Dirac Synchronization is characterized by discontinuous transitions and the emergence of a rhythmic coherent phase. In this rhythmic phase, one of the two complex order parameters oscillates in the complex plane at a slow frequency (called emergent frequency) in the frame in which the intrinsic frequencies have zero average. Our theoretical results obtained within the annealed approximation are validated by extensive numerical results on fully connected networks and sparse Poisson and scale-free networks. Local Dirac Synchronization on both random and real networks, such as the connectome of Caenorhabditis Elegans, reveals the interplay between topology (Betti numbers and harmonic modes) and non-linear dynamics. This unveils how topology might play a role in the onset of brain rhythms.

Phase frustration induced remote synchronization

Remote synchronization (RS) may take an important role in brain functioning and its study has attracted much attention in recent years. So far, most studies of RS are focused on the Stuart-Landau oscillators with mean-field coupling. However, realistic cases may have more complicated couplings and behaviors, such as the brain networks. To make the study of RS a substantial progress toward realistic situations, we here present a model of RS with phase frustration and show that RS can be induced for those systems where no RS exists when there is no phase frustration. By numerical simulations on both the Stuart-Landau and Kuramoto oscillators, we find that the optimal range of RS depends on the match of phase frustrations between the hub and leaf nodes and a fixed relationship of this match is figured out. While for the non-optimal range of RS, we find that RS exists only in a linear band between the phase frustrations of the hub and leaf nodes. A brief theoretical analysis is provided to explain these results.

Metamorphoses and explosively remote synchronization in dynamical networks

We uncover a phenomenon in coupled nonlinear networks with a symmetry: as a bifurcation parameter changes through a critical value, synchronization among a subset of nodes can deteriorate abruptly, and, simultaneously, perfect synchronization emerges suddenly among a different subset of nodes that are not directly connected. This is a synchronization metamorphosis leading to an explosive transition to remote synchronization. The finding demonstrates that an explosive onset of synchrony and remote synchronization, two phenomena that have been studied separately, can arise in the same system due to symmetry, providing another proof that the interplay between nonlinear dynamics and symmetry can lead to a surprising phenomenon in physical systems.

The transition to synchronization on branching hierarchical lattices

We study the transition to synchronization on hierarchical lattices using the evolution of Chaté-Manneville maps placed on a triangular lattice. Connections are generated between the levels of the triangular lattice, assuming that each site is connected to its neighbors on the level below with probability half. The maps are diffusively coupled, and the map parameters increase hierarchically, depending on the map parameters at the sites they are coupled to in the previous level. The system shows a transition to synchronization, which is second order in nature, with associated critical exponents. However, the V-lattice, which is a special realization of this lattice, shows a transition to synchronization that is discontinuous with accompanying hysteretic behavior. This transition can thus be said to belong to the class of explosive synchronization with the explosive nature depending on the nature of the substrate. We carry out finite-size-finite-time scaling for the continuous transition and analyze the scaling of the jump size for the discontinuous case. We discuss the implications of our results and draw parallels with avalanche statistics on branching hierarchical lattices.

Basins with Tentacles

To explore basin geometry in high-dimensional dynamical systems, we consider a ring of identical Kuramoto oscillators. Many attractors coexist in this system; each is a twisted periodic orbit characterized by a winding number q, with basin size proportional to e^{-kq^{2}}. We uncover the geometry behind this size distribution and find the basins are octopuslike, with nearly all their volume in the tentacles, not the head of the octopus (the ball-like region close to the attractor). We present a simple geometrical reason why basins with tentacles should be common in high-dimensional systems.

Edges of inter-layer synchronization in multilayer networks with time-switching links

We investigate the transition to synchronization in a two-layer network of oscillators with time-switching inter-layer links. We focus on the role of the number of inter-layer links and the timescale of topological changes. Initially, we observe a smooth transition to complete synchronization for the static inter-layer topology by increasing the number of inter-layer links. Next, for a dynamic topology with the existent inter-layer links randomly changing among identical oscillators in the layers, we observe a significant improvement in the system synchronizability; i.e., the layers synchronize with lower inter-layer connectivity. More interestingly, we find that, for a critical switching time, the transition from the network state of low inter-layer synchronization to high inter-layer synchronization occurs abruptly as the number of inter-layer links increases. We interpret this phenomenon as shrinking and ultimately the disappearance of the basin of attraction of a desynchronized network state.

Coherence resonance in influencer networks

Complex networks are abundant in nature and many share an important structural property: they contain a few nodes that are abnormally highly connected (hubs). Some of these hubs are called influencers because they couple strongly to the network and play fundamental dynamical and structural roles. Strikingly, despite the abundance of networks with influencers, little is known about their response to stochastic forcing. Here, for oscillatory dynamics on influencer networks, we show that subjecting influencers to an optimal intensity of noise can result in enhanced network synchronization. This new network dynamical effect, which we call coherence resonance in influencer networks, emerges from a synergy between network structure and stochasticity and is highly nonlinear, vanishing when the noise is too weak or too strong. Our results reveal that the influencer backbone can sharply increase the dynamical response in complex systems of coupled oscillators.

Co-adaptation enhances the resilience of mutualistic networks

Mutualistic networks, which describe the ecological interactions between multiple types of species such as plants and pollinators, play a paramount role in the generation of Earth's biodiversity. The resilience of a mutualistic network denotes its ability to retain basic functionality when errors and failures threaten the persistence of the community. Under the disturbances of mass extinctions and human-induced disasters, it is crucial to understand how mutualistic networks respond to changes, which enables the system to increase resilience and tolerate further damages. Despite recent advances in the modelling of the structure-based adaptation, we lack mathematical and computational models to describe and capture the co-adaptation between the structure and dynamics of mutualistic networks. In this paper, we incorporate dynamic features into the adaptation of structure and propose a co-adaptation model that drastically enhances the resilience of non-adaptive and structure-based adaptation models. Surprisingly, the reason for the enhancement is that the co-adaptation mechanism simultaneously increases the heterogeneity of the mutualistic network significantly without changing its connectance. Owing to the broad applications of mutualistic networks, our findings offer new ways to design mechanisms that enhance the resilience of many other systems, such as smart infrastructures and social-economical systems.

Hysteresis and criticality in hybrid percolation transitions

Phase transitions (PTs) are generally classified into second-order and first-order transitions, each exhibiting different intrinsic properties. For instance, a first-order transition exhibits latent heat and hysteresis when a control parameter is increased and then decreased across a transition point, whereas a second-order transition does not. Recently, hybrid percolation transitions (HPTs) are issued in diverse complex systems, in which the features of first-order and second-order PTs occur at the same transition point. Thus, the question whether hysteresis appears in an HPT arises. Herein, we investigate this fundamental question with a so-called restricted Erdős-Rényi random network model, in which a cluster fragmentation process is additionally proposed. A hysteresis curve of the order parameter was obtained. Depending on when the reverse process is initiated, the shapes of hysteresis curves change, and the critical behavior of the HPT is conserved throughout the forward and reverse processes.

Amplification of explosive width in complex networks

We present an adaptive coupling strategy to induce hysteresis/explosive synchronization in complex networks of phase oscillators (Sakaguchi-Kuramoto model). The coupling strategy ensures explosive synchronization with significant explosive width enhancement. Results show the robustness of the strategy, and the strategy can diminish (by inducing enhanced hysteresis loop) the contrarian impact of phase frustration in the network, irrespective of the network structure or frequency distributions. Additionally, we design a set of frequency for the oscillators, which eventually ensure complete in-phase synchronization behavior among these oscillators (with enhanced explosive width) in the case of adaptive-coupling scheme. Based on a mean-field analysis, we develop a semi-analytical formalism, which can accurately predict the backward transition of the synchronization order parameter.

Reduction of oscillator dynamics on complex networks to dynamics on complete graphs through virtual frequencies

We consider the synchronization of oscillators in complex networks where there is an interplay between the oscillator dynamics and the network topology. Through a remarkable transformation in parameter space and the introduction of virtual frequencies we show that Kuramoto oscillators on annealed networks, with or without frequency-degree correlation, and Kuramoto oscillators on complete graphs with frequency-weighted coupling can be transformed to Kuramoto oscillators on complete graphs with a rearranged, virtual frequency distribution and uniform coupling. The virtual frequency distribution encodes both the natural frequency distribution (dynamics) and the degree distribution (topology). We apply this transformation to give direct explanations to a variety of phenomena that have been observed in complex networks, such as explosive synchronization and vanishing synchronization onset.

Effects of network robustness on explosive synchronization

Current studies have shown that there is a positive correlation between the network assortativity and robustness and that the assortativity also plays an important role in explosive synchronization. In this paper, taking the network robustness as a global property, we investigate its significance as well as the influence of its interaction with the assortativity on explosive synchronization. Our numerical results demonstrate that explosive synchronization is suppressed in extreme situations of both the robustness and assortativity. In addition, through appropriate adjustments of them, a maximum hysteresis area between the forward and backward transitions can be reached. Furthermore, our results might also provide reference for those who are interested in effects of network structure on synchronization, though this problem is still challenging as we show in the discussion.

Synchronization in starlike networks of phase oscillators

We fully describe the mechanisms underlying synchronization in starlike networks of phase oscillators. In particular, the routes to synchronization and the critical points for the associated phase transitions are determined analytically. In contrast to the classical Kuramoto theory, we unveil that relaxation rates to each equilibrium state indeed exist and remain invariant under three levels of descriptions corresponding to different geometric implications. The special symmetry in the coupling determines a quasi-Hamiltonian property, which is further unveiled on the basis of singular perturbation theory. Since starlike coupling configurations constitute the building blocks of technological and biological real world networks, our paper paves the way towards the understanding of the functioning of such real world systems in many practical situations.

Phase transition to synchronization in generalized Kuramoto model with low-pass filter

A second-order continuous synchronization has been well documented for the classic Kuramoto model. Here we generalize the classic Kuramoto model by incorporating a low-pass filter (LPF) in the coupling, which serves as a simple form of indirect coupling through a common external dynamic environment. We uncover that a first-order explosive synchronization turns out to be a very generic phenomenon in this generalized Kuramoto model with LPF. We establish theoretical results by providing a rigorous analytical treatment, which is validated by conducting extensive numerical simulations. Our study provides a new root for the emergence of first-order explosive synchronization, which could substantially deepen the understanding of the underlying mechanism of a first-order phase transition towards synchronization in coupled dynamical networks.

Synchronization of interconnected heterogeneous networks: The role of network sizes

Increasing evidence shows that real networks interact with each other, forming a network of networks (NONs). Synchronization, a ubiquitous process in natural and engineering systems, has fascinatingly gained rising attentions in the context of NONs. Despite efforts to study the synchronization of NONs, it is still a challenge to understand how do the network sizes affect the synchronization and its phase diagram of NONs coupled with nonlinear dynamics. Here, we model such NONs as star-like motifs to analytically derive the critical values of both the internal and the external coupling strengths, at which a phase transition from synchronization to incoherence occurs. Our results show that the critical values strongly depend on the network sizes. Reducing the difference between network sizes will enhance the synchronization of the whole system, which indicates the irrationality of previous studies that assume the network sizes to be the same. The optimal connection strategy also changes as the network sizes change, a discovery contradicting to the previous conclusion that connecting the high-degree nodes of each network is always the most effective strategy to achieve synchronization unchangeably. This finding emphasizes the crucial role of network sizes which has been neglected in the previous studies and could contribute to the design of a global synchronized system.

Time dependent stability margin in multistable systems

We propose a novel technique to analyze multistable, non-linear dynamical systems. It enables one to characterize the evolution of a time-dependent stability margin along stable periodic orbits. By that, we are able to indicate the moments along the trajectory when the stability surplus is minimal, and even relatively small perturbation can lead to a tipping point. We explain the proposed approach using two paradigmatic dynamical systems, i.e., Rössler and Duffing oscillators. Then, the method is validated experimentally using the rig with a double pendulum excited parametrically. Both numerical and experimental results reveal significant fluctuations of sensitivity to perturbations along the considered periodic orbits. The proposed concept can be used in multiple applications including engineering, fluid dynamics, climate research, and photonics.

Explosive death of conjugate coupled Van der Pol oscillators on networks

Explosive death phenomenon has been gradually gaining attention of researchers due to the research boom of explosive synchronization, and it has been observed recently for the identical or nonidentical coupled systems in all-to-all network. In this work, we investigate the emergence of explosive death in networked Van der Pol (VdP) oscillators with conjugate variables coupling. It is demonstrated that the network structures play a crucial role in identifying the types of explosive death behaviors. We also observe that the damping coefficient of the VdP system not only can determine whether the explosive death state is generated but also can adjust the forward transition point. We further show that the backward transition point is independent of the network topologies and the damping coefficient, which is well confirmed by theoretical analysis. Our results reveal the generality of explosive death phenomenon in different network topologies and are propitious to promote a better comprehension for the oscillation quenching behaviors.

Antagonistic Phenomena in Network Dynamics

Recent research on the network modeling of complex systems has led to a convenient representation of numerous natural, social, and engineered systems that are now recognized as networks of interacting parts. Such systems can exhibit a wealth of phenomena that not only cannot be anticipated from merely examining their parts, as per the textbook definition of complexity, but also challenge intuition even when considered in the context of what is now known in network science. Here we review the recent literature on two major classes of such phenomena that have far-reaching implications: (i) antagonistic responses to changes of states or parameters and (ii) coexistence of seemingly incongruous behaviors or properties-both deriving from the collective and inherently decentralized nature of the dynamics. They include effects as diverse as negative compressibility in engineered materials, rescue interactions in biological networks, negative resistance in fluid networks, and the Braess paradox occurring across transport and supply networks. They also include remote synchronization, chimera states and the converse of symmetry breaking in brain, power-grid and oscillator networks as well as remote control in biological and bio-inspired systems. By offering a unified view of these various scenarios, we suggest that they are representative of a yet broader class of unprecedented network phenomena that ought to be revealed and explained by future research.

The role of community structure on the nature of explosive synchronization

In this paper, we analyze explosive synchronization in networks with a community structure. The results of our study indicate that the mesoscopic structure of the networks could affect the synchronization of coupled oscillators. With the variation of three parameters, the degree probability distribution exponent, the community size probability distribution exponent, and the mixing parameter, we could have a fast or slow phase transition. Besides, in some cases, we could have communities which are synchronized inside but not with other communities and vice versa. We also show that there is a limit in these mesoscopic structures which suppresses the transition from the second-order phase transition and results in explosive synchronization. This could be considered as a tuning parameter changing the transition of the system from the second order to the first order.

Influence of stochastic perturbations on the cluster explosive synchronization of second-order Kuramoto oscillators on networks

Explosive synchronization in networked second-order Kuramoto oscillators has been well studied recently and it is revealed that the synchronization process is featured by cluster explosive synchronization. However, little attention has been paid to the influence of noise or perturbation. We here study this problem and discuss the influences of noise and perturbation. For the former, we interestingly find that noise has significant influence on the cluster explosive synchronization of those nodes with smaller degrees, i.e., their synchronization will change from the first-order to second-order transition and the critical points for both the forward and backward synchronization depend on the strength of noise. Especially, when the strength of noise is in an optimal range, a synchronization of the nodes with smaller degrees will be induced in the region of coupling strength where they do not display synchronization in the absence of noise. For the latter, we find that the effect of perturbation is similar to that of noise when its duration W is small. However, the perturbation will induce a change from cluster explosive synchronization to explosive synchronization when W is large. Furthermore, a brief theory is provided to explain the influence of perturbations on the critical points.

Functional Brain Network Mechanism of Hypersensitivity in Chronic Pain

Fibromyalgia (FM) is a chronic widespread pain condition characterized by augmented multi-modal sensory sensitivity. Although the mechanisms underlying this sensitivity are thought to involve an imbalance in excitatory and inhibitory activity throughout the brain, the underlying neural network properties associated with hypersensitivity to pain stimuli are largely unknown. In network science, explosive synchronization (ES) was introduced as a mechanism of hypersensitivity in diverse biological and physical systems that display explosive and global propagations with small perturbations. We hypothesized that ES may also be a mechanism of the hypersensitivity in FM brains. To test this hypothesis, we analyzed resting state electroencephalogram (EEG) of 10 FM patients. First, we examined theoretically well-known ES conditions within functional brain networks reconstructed from EEG, then tested whether a brain network model with ES conditions identified in the EEG data is sensitive to an external perturbation. We demonstrate for the first time that the FM brain displays characteristics of ES conditions, and that these factors significantly correlate with chronic pain intensity. The simulation data support the conclusion that networks with ES conditions are more sensitive to perturbation compared to non-ES network. The model and empirical data analysis provide convergent evidence that ES may be a network mechanism of FM hypersensitivity.

Self-similarity in explosive synchronization of complex networks

We report that explosive synchronization of networked oscillators (a process through which the transition to coherence occurs without intermediate stages but is rather characterized by a sudden and abrupt jump from the network's asynchronous to synchronous motion) is related to self-similarity of synchronous clusters of different size. Self-similarity is revealed by destructing the network synchronous state during the backward transition and observed with the decrease of the coupling strength between the nodes of the network. As illustrative examples, networks of Kuramoto oscillators with different topologies of links have been considered. For each one of such topologies, the destruction of the synchronous state goes step by step with self-similar configurations of interacting oscillators. At the critical point, the invariance of the phase distribution in the synchronized cluster with respect to the cluster size is reported.

Connection adaption for control of networked mobile chaotic agents

In this paper, we propose a strategy for the control of mobile chaotic oscillators by adaptively rewiring connections between nearby agents with local information. In contrast to the dominant adaptive control schemes where coupling strength is adjusted continuously according to the states of the oscillators, our method does not request adaption of coupling strength. As the resulting interaction structure generated by this proposed strategy is strongly related to unidirectional chains, by investigating synchronization property of unidirectional chains, we reveal that there exists a certain coupling range in which the agents could be controlled regardless of the length of the chain. This feature enables the adaptive strategy to control the mobile oscillators regardless of their moving speed. Compared with existing adaptive control strategies for networked mobile agents, our proposed strategy is simpler for implementation where the resulting interaction networks are kept unweighted at all time.

Basin stability measure of different steady states in coupled oscillators

In this report, we investigate the stabilization of saddle fixed points in coupled oscillators where individual oscillators exhibit the saddle fixed points. The coupled oscillators may have two structurally different types of suppressed states, namely amplitude death and oscillation death. The stabilization of saddle equilibrium point refers to the amplitude death state where oscillations are ceased and all the oscillators converge to the single stable steady state via inverse pitchfork bifurcation. Due to multistability features of oscillation death states, linear stability theory fails to analyze the stability of such states analytically, so we quantify all the states by basin stability measurement which is an universal nonlocal nonlinear concept and it interplays with the volume of basins of attractions. We also observe multi-clustered oscillation death states in a random network and measure them using basin stability framework. To explore such phenomena we choose a network of coupled Duffing-Holmes and Lorenz oscillators which are interacting through mean-field coupling. We investigate how basin stability for different steady states depends on mean-field density and coupling strength. We also analytically derive stability conditions for different steady states and confirm by rigorous bifurcation analysis.

A small change in neuronal network topology can induce explosive synchronization transition and activity propagation in the entire network

We here study explosive synchronization transitions and network activity propagation in networks of coupled neurons to provide a new understanding of the relationship between network topology and explosive dynamical transitions as in epileptic seizures and their propagations in the brain. We model local network motifs and configurations of coupled neurons and analyze the activity propagations between a group of active neurons to their inactive neuron neighbors in a variety of network configurations. We find that neuronal activity propagation is limited to local regions when network is highly clustered with modular structures as in the normal brain networks. When the network cluster structure is slightly changed, the activity propagates to the entire network, which is reminiscent of epileptic seizure propagation in the brain. Finally, we analyze intracranial electroencephalography (IEEG) recordings of a seizure episode from a epilepsy patient and uncover that explosive synchronization-like transition occurs around the clinically defined onset of seizure. These findings may provide a possible mechanism for the recurrence of epileptic seizures, which are known to be the results of aberrant neuronal network structure and/or function in the brain.

Synchronization in slowly switching networks of coupled oscillators

Networks whose structure of connections evolves in time constitute a big challenge in the study of synchronization, in particular when the time scales for the evolution of the graph topology are comparable with (or even longer than) those pertinent to the units’ dynamics. We here focus on networks with a slow-switching structure, and show that the necessary conditions for synchronization, i.e. the conditions for which synchronization is locally stable, are determined by the time average of the largest Lyapunov exponents of transverse modes of the switching topologies. Comparison between fast- and slow-switching networks allows elucidating that slow-switching processes prompt synchronization in the cases where the Master Stability Function is concave, whereas fast-switching schemes facilitate synchronization for convex curves. Moreover, the condition of slow-switching enables the introduction of a control strategy for inducing synchronization in networks with arbitrary structure and coupling strength, which is of evident relevance for broad applications in real world systems.

Collective dynamics of identical phase oscillators with high-order coupling

In this paper, we propose a framework to investigate the collective dynamics in ensembles of globally coupled phase oscillators when higher-order modes dominate the coupling. The spatiotemporal properties of the attractors in various regions of parameter space are analyzed. Furthermore, a detailed linear stability analysis proves that the stationary symmetric distribution is only neutrally stable in the marginal regime which stems from the generalized time-reversal symmetry. Moreover, the critical parameters of the transition among various regimes are determined analytically by both the Ott-Antonsen method and linear stability analysis, the transient dynamics are further revealed in terms of the characteristic curves method. Finally, for the more general initial condition the symmetric dynamics could be reduced to a rigorous three-dimensional manifold which shows that the neutrally stable chaos could also occur in this model for particular parameters. Our theoretical analysis and numerical results are consistent with each other, which can help us understand the dynamical properties in general systems with higher-order harmonics couplings.

Model bridging chimera state and explosive synchronization

Global synchronization and partial synchronization are the two distinctive forms of synchronization in coupled oscillators and have been well studied in recent decades. Recent attention on synchronization is focused on the chimera state (CS) and explosive synchronization (ES), but little attention has been paid to their relationship. Here we study this topic by presenting a model to bridge these two phenomena, which consists of two groups of coupled oscillators, and its coupling strength is adaptively controlled by a local order parameter. We find that this model displays either CS or ES in two limits. In between the two limits, this model exhibits both CS and ES, where CS can be observed for a fixed coupling strength and ES appears when the coupling is increased adiabatically. Moreover, we show both theoretically and numerically that there are a variety of CS basin patterns for the case of identical oscillators, depending on the distributions of both the initial order parameters and the initial average phases. This model suggests a way to easily observe CS, in contrast to other models having some (weak or strong) dependence on initial conditions.

Explosive synchronization coexists with classical synchronization in the Kuramoto model

Explosive synchronization has recently been reported in a system of adaptively coupled Kuramoto oscillators, without any conditions on the frequency or degree of the nodes. Here, we find that, in fact, the explosive phase coexists with the standard phase of the Kuramoto oscillators. We determine this by extending the mean-field theory of adaptively coupled oscillators with full coupling to the case with partial coupling of a fraction f. This analysis shows that a metastable region exists for all finite values of f > 0, and therefore explosive synchronization is expected for any perturbation of adaptively coupling added to the standard Kuramoto model. We verify this theory with GPU-accelerated simulations on very large networks (N ∼ 10(6)) and find that, in fact, an explosive transition with hysteresis is observed for all finite couplings. By demonstrating that explosive transitions coexist with standard transitions in the limit of f → 0, we show that this behavior is far more likely to occur naturally than was previously believed.

Synchronous slowing down in coupled logistic maps via random network topology

The speed and paths of synchronization play a key role in the function of a system, which has not received enough attention up to now. In this work, we study the synchronization process of coupled logistic maps that reveals the common features of low-dimensional dissipative systems. A slowing down of synchronization process is observed, which is a novel phenomenon. The result shows that there are two typical kinds of transient process before the system reaches complete synchronization, which is demonstrated by both the coupled multiple-period maps and the coupled multiple-band chaotic maps. When the coupling is weak, the evolution of the system is governed mainly by the local dynamic, i.e., the node states are attracted by the stable orbits or chaotic attractors of the single map and evolve toward the synchronized orbit in a less coherent way. When the coupling is strong, the node states evolve in a high coherent way toward the stable orbit on the synchronized manifold, where the collective dynamics dominates the evolution. In a mediate coupling strength, the interplay between the two paths is responsible for the slowing down. The existence of different synchronization paths is also proven by the finite-time Lyapunov exponent and its distribution.

Synchronization of phase oscillators with frequency-weighted coupling

Recently, the first-order synchronization transition has been studied in systems of coupled phase oscillators. In this paper, we propose a framework to investigate the synchronization in the frequency-weighted Kuramoto model with all-to-all couplings. A rigorous mean-field analysis is implemented to predict the possible steady states. Furthermore, a detailed linear stability analysis proves that the incoherent state is only neutrally stable below the synchronization threshold. Nevertheless, interestingly, the amplitude of the order parameter decays exponentially (at least for short time) in this regime, resembling the Landau damping effect in plasma physics. Moreover, the explicit expression for the critical coupling strength is determined by both the mean-field method and linear operator theory. The mechanism of bifurcation for the incoherent state near the critical point is further revealed by the amplitude expansion theory, which shows that the oscillating standing wave state could also occur in this model for certain frequency distributions. Our theoretical analysis and numerical results are consistent with each other, which can help us understand the synchronization transition in general networks with heterogenous couplings.

Basin stability for burst synchronization in small-world networks of chaotic slow-fast oscillators

The impact of connectivity and individual dynamics on the basin stability of the burst synchronization regime in small-world networks consisting of chaotic slow-fast oscillators is studied. It is shown that there are rewiring probabilities corresponding to the largest basin stabilities, which uncovers a reason for finding small-world topologies in real neuronal networks. The impact of coupling density and strength as well as the nodal parameters of relaxation or excitability are studied. Dynamic mechanisms are uncovered that most strongly influence basin stability of the burst synchronization regime.

Explosive or Continuous: Incoherent state determines the route to synchronization

Abrupt and continuous spontaneous emergence of collective synchronization of coupled oscillators have attracted much attention. In this paper, we propose a dynamical ensemble order parameter equation that enables us to grasp the essential low-dimensional dynamical mechanism of synchronization in networks of coupled oscillators. Different solutions of the dynamical ensemble order parameter equation build correspondences with diverse collective states, and different bifurcations reveal various transitions among these collective states. The structural relationship between the incoherent state and the synchronous state leads to different routes of transitions to synchronization, either continuous or discontinuous. The explosive synchronization is determined by the bistable state where the measure of each state and the critical points are obtained analytically by using the dynamical ensemble order parameter equation. Our method and results hold for heterogeneous networks with star graph motifs such as scale-free networks, and hence, provide an effective approach in understanding the routes to synchronization in more general complex networks.

Explosive synchronization with asymmetric frequency distribution

In this work, we study the synchronization in a generalized Kuramoto model with frequency-weighted coupling. In particular, we focus on the situations in which the frequency distributions of oscillators are asymmetric. For typical unimodal frequency distributions, such as Lorentzian, Gaussian, triangle, and even special Rayleigh, we find that the synchronization transition in the model generally converts from the first order to the second order as the central frequency shifts toward positive direction. We characterize two interesting coherent states in the system: In the former, two phase-locking clusters are formed, rotating with the same frequency. They correspond to those oscillators with relatively high frequencies while the oscillators with relatively small frequencies are not entrained. In the latter, two phase-locking clusters rotate with different frequencies, leading to the oscillation of the order parameter. We further conduct theoretical analysis to reveal the relation between the asymmetric frequency distribution and the conversion of synchronization type, and justify the coherent states observed in the system.

Practical synchronization on complex dynamical networks via optimal pinning control

We consider practical synchronization on complex dynamical networks under linear feedback control designed by optimal control theory. The control goal is to minimize global synchronization error and control strength over a given finite time interval, and synchronization error at terminal time. By utilizing the Pontryagin's minimum principle, and based on a general complex dynamical network, we obtain an optimal system to achieve the control goal. The result is verified by performing some numerical simulations on Star networks, Watts-Strogatz networks, and Barabási-Albert networks. Moreover, by combining optimal control and traditional pinning control, we propose an optimal pinning control strategy which depends on the network's topological structure. Obtained results show that optimal pinning control is very effective for synchronization control in real applications.

Explosive synchronization is discontinuous

Spontaneous explosive is an abrupt transition to collective behavior taking place in heterogeneous networks when the frequencies of the nodes are positively correlated with the node degree. This explosive transition was conjectured to be discontinuous. Indeed, numerical investigations reveal a hysteresis behavior associated with the transition. Here, we analyze explosive synchronization in star graphs. We show that in the thermodynamic limit the transition to (and out of) collective behavior is indeed discontinuous. The discontinuous nature of the transition is related to the nonlinear behavior of the order parameter, which in the thermodynamic limit exhibits multiple fixed points. Moreover, we unravel the hysteresis behavior in terms of the graph parameters. Our numerical results show that finite-size graphs are well described by our predictions.

Effects of assortative mixing in the second-order Kuramoto model

In this paper we analyze the second-order Kuramoto model in the presence of a positive correlation between the heterogeneity of the connections and the natural frequencies in scale-free networks. We numerically show that discontinuous transitions emerge not just in disassortative but also in strongly assortative networks, in contrast with the first-order model. We also find that the effect of assortativity on network synchronization can be compensated by adjusting the phase damping. Our results show that it is possible to control collective behavior of damped Kuramoto oscillators by tuning the network structure or by adjusting the dissipation related to the phases' movement.

Growing optimal scale-free networks via likelihood

Preferential attachment, by which new nodes attach to existing nodes with probability proportional to the existing nodes' degree, has become the standard growth model for scale-free networks, where the asymptotic probability of a node having degree k is proportional to k^{-γ}. However, the motivation for this model is entirely ad hoc. We use exact likelihood arguments and show that the optimal way to build a scale-free network is to attach most new links to nodes of low degree. Curiously, this leads to a scale-free network with a single dominant hub: a starlike structure we call a superstar network. Asymptotically, the optimal strategy is to attach each new node to one of the nodes of degree k with probability proportional to 1/N+ζ(γ)(k+1)(γ) (in a N node network): a stronger bias toward high degree nodes than exhibited by standard preferential attachment. Our algorithm generates optimally scale-free networks (the superstar networks) as well as randomly sampling the space of all scale-free networks with a given degree exponent γ. We generate viable realization with finite N for 1≪γ<2 as well as γ>2. We observe an apparently discontinuous transition at γ≈2 between so-called superstar networks and more treelike realizations. Gradually increasing γ further leads to reemergence of a superstar hub. To quantify these structural features, we derive a new analytic expression for the expected degree exponent of a pure preferential attachment process and introduce alternative measures of network entropy. Our approach is generic and can also be applied to an arbitrary degree distribution.

Analytical description for the critical fixations of evolutionary coordination games on finite complex structured populations

Evolutionary game theory is crucial to capturing the characteristic interaction patterns among selfish individuals. In a population of coordination games of two strategies, one of the central problems is to determine the fixation probability that the system reaches a state of networkwide of only one strategy, and the corresponding expectation times. The deterministic replicator equations predict the critical value of initial density of one strategy, which separates the two absorbing states of the system. However, numerical estimations of this separatrix show large deviations from the theory in finite populations. Here we provide a stochastic treatment of this dynamic process on complex networks of finite sizes as Markov processes, showing the evolutionary time explicitly. We describe analytically the effects of network structures on the intermediate fixations as observed in numerical simulations. Our theoretical predictions are validated by various simulations on both random and scale free networks. Therefore, our stochastic framework can be helpful in dealing with other networked game dynamics.

Effects of degree correlations on the explosive synchronization of scale-free networks

We study the organization of finite-size, large ensembles of phase oscillators networking via scale-free topologies in the presence of a positive correlation between the oscillators' natural frequencies and the network's degrees. Under those circumstances, abrupt transitions to synchronization are known to occur in growing scale-free networks, while the transition has a completely different nature for static random configurations preserving the same structure-dynamics correlation. We show that the further presence of degree-degree correlations in the network structure has important consequences on the nature of the phase transition characterizing the passage from the phase-incoherent to the phase-coherent network state. While high levels of positive and negative mixing consistently induce a second-order phase transition, moderate values of assortative mixing, such as those ubiquitously characterizing social networks in the real world, greatly enhance the irreversible nature of explosive synchronization in scale-free networks. The latter effect corresponds to a maximization of the area and of the width of the hysteretic loop that differentiates the forward and backward transitions to synchronization.

Self-organized correlations lead to explosive synchronization

Very recently, a first-order phase transition, named explosive synchronization (ES), has attracted great attention due to its remarkable novelty in theory and significant impact in applications. However, so far, all observations of ES have been associated with various correlation constraints on system parameters, which restrict its generality and applications. Here we consider heterogeneous networks around Hopf bifurcation point described by chemical reaction-diffusion systems and also by their reduced order parameter versions, the complex Ginzburg-Landau equations, and demonstrate that explosive synchronization can appear as an emergent feature of oscillatory networks, and the restrictions on specific parameter correlations used so far for ES can be lifted entirely. Theoretical analyses and numerical simulations show with a perfect agreement that explosive synchronization can appear in networks with nodes having identical natural frequencies, and necessary correlation conditions for ES can be realized in a self-organized manner by network evolution.

Quantification of causal couplings via dynamical effects: a unifying perspective

Quantitative characterization of causal couplings from time series is crucial in studies of complex systems of different origin. Various statistical tools for that exist and new ones are still being developed with a tendency to creating a single, universal, model-free quantifier of coupling strength. However, a clear and generally applicable way of interpreting such universal characteristics is lacking. This work suggests a general conceptual framework for causal coupling quantification, which is based on state space models and extends the concepts of virtual interventions and dynamical causal effects. Namely, two basic kinds of interventions (state space and parametric) and effects (orbital or transient and stationary or limit) are introduced, giving four families of coupling characteristics. The framework provides a unifying view of apparently different well-established measures and allows us to introduce new characteristics, always with a definite "intervention-effect" interpretation. It is shown that diverse characteristics cannot be reduced to any single coupling strength quantifier and their interpretation is inevitably model based. The proposed set of dynamical causal effect measures quantifies different aspects of "how the coupling manifests itself in the dynamics," reformulating the very question about the "causal coupling strength."

Exact solution for first-order synchronization transition in a generalized Kuramoto model

First-order, or discontinuous, synchronization transition, i.e. an abrupt and irreversible phase transition with hysteresis to the synchronized state of coupled oscillators, has attracted much attention along the past years. We here report the analytical solution of a generalized Kuramoto model, and derive a series of exact results for the first-order synchronization transition, including i) the exact, generic, solutions for the critical coupling strengths for both the forward and backward transitions, ii) the closed form of the forward transition point and the linear stability analysis for the incoherent state (for a Lorentzian frequency distribution), and iii) the closed forms for both the stable and unstable coherent states (and their stabilities) for the backward transition. Our results, together with elucidating the first-order nature of the transition, provide insights on the mechanisms at the basis of such a synchronization phenomenon.

Analysis of cluster explosive synchronization in complex networks

Correlations between intrinsic dynamics and local topology have become a new trend in the study of synchronization in complex networks. In this paper, we investigate the influence of topology on the dynamics of networks made up of second-order Kuramoto oscillators. In particular, based on mean-field calculations, we provide a detailed investigation of cluster explosive synchronization (CES) [Phys. Rev. Lett. 110, 218701 (2013)] in scale-free networks as a function of several topological properties. Moreover, we investigate the robustness of discontinuous transitions by including an additional quenched disorder, and we show that the phase coherence decreases with increasing strength of the quenched disorder. These results complement the previous findings regarding CES and also fundamentally deepen the understanding of the interplay between topology and dynamics under the constraint of correlating natural frequencies and local structure.

Synchronization in time-varying networks

We study the stability of the synchronized state in time-varying complex networks using the concept of basin stability, which is a nonlocal and nonlinear measure of stability that can be easily applied to high-dimensional systems [P. J. Menck, J. Heitzig, N. Marwan, and J. Kurths, Nature Phys. 9, 89 (2013)]. The time-varying character is included by stochastically rewiring each link with the average frequency f. We find that the time taken to reach synchronization is lowered and the stability range of the synchronized state increases considerably in dynamic networks. Further we uncover that small-world networks are much more sensitive to link changes than random ones, with the time-varying character of the network having a significant effect at much lower rewiring frequencies. At very high rewiring frequencies, random networks perform better than small-world networks and the synchronized state is stable over a much wider window of coupling strengths. Lastly we show that the stability range of the synchronized state may be quite different for small and large perturbations, and so the linear stability analysis and the basin stability criterion provide complementary indicators of stability.