Reexamination of explosive synchronization in scale-free networks: the effect of disassortativity

Hysteresis and synchronization processes of Kuramoto oscillators on high-dimensional simplicial complexes with competing simplex-encoded couplings

Recent studies of dynamic properties in complex systems point out the profound impact of hidden geometry features known as simplicial complexes, which enable geometrically conditioned many-body interactions. Studies of collective behaviors on the controlled-structure complexes can reveal the subtle interplay of geometry and dynamics. Here we investigate the phase synchronization (Kuramoto) dynamics under the competing interactions embedded on 1-simplex (edges) and 2-simplex (triangles) faces of a homogeneous four-dimensional simplicial complex. Its underlying network is a 1-hyperbolic graph with the assortative correlations among the node's degrees and the spectral dimension that exceeds d_{s}=4. By numerically solving the set of coupled equations for the phase oscillators associated with the network nodes, we determine the time-averaged system's order parameter to characterize the synchronization level. Our results reveal a variety of synchronization and desynchronization scenarios, including partially synchronized states and nonsymmetrical hysteresis loops, depending on the sign and strength of the pairwise interactions and the geometric frustrations promoted by couplings on triangle faces. For substantial triangle-based interactions, the frustration effects prevail, preventing the complete synchronization and the abrupt desynchronization transition disappears. These findings shed new light on the mechanisms by which the high-dimensional simplicial complexes in natural systems, such as human connectomes, can modulate their native synchronization processes.

Assortativity-induced explosive synchronization in a complex neuronal network

In this study, we consider a scale-free network of nonidentical Chialvo neurons, coupled through electrical synapses. For sufficiently strong coupling, the system undergoes a transition from completely out of phase synchronized to phase synchronized state. The principal focus of this study is to investigate the effect of the degree of assortativity over the synchronization transition process. It is observed that, depending on assortativity, bistability between two asymptotically stable states allows one to develop a hysteresis loop which transforms the phase transition from second order to first order. An expansion in the area of hysteresis loop is noticeable with increasing degree-degree correlation in the network. Our study also reveals that effective frequencies of nodes simultaneously go through a continuous or sudden transition to the synchronized state with the corresponding phases. Further, we examine the robustness of the results under the effect of network size and average degree, as well as diverse frequency setup. Finally, we investigate the dynamical mechanism in the process of generating explosive synchronization. We observe a significant impact of lower degree nodes behind such phenomena: in a positive assortative network the low degree nodes delay the synchronization transition.

Mechanism for explosive synchronization of neural networks

Here we investigate the mechanism for explosive synchronization (ES) of a complex neural network composed of nonidentical neurons and coupled by Newman-Watts small-world matrices. We find a range of nonlocal connection probabilities for which the network displays an abrupt transition to phase synchronization, characterizing ES. The mechanism behind the ES is the following: As the coupling parameter is varied in a network of distinct neurons, ES is likely to occur due to a bistable regime, namely a chaotic nonsynchronized and a regular phase-synchronized state in the phase space. In this case, even small coupling changes make possible a transition between them. The onset of ES occurs via a saddle-node bifurcation of a periodic orbit that leads the network dynamics to display a locally stable phase-synchronized state. The presence of this regime is accompanied by a hysteresis loop on the network dynamics as the coupling parameter is adiabatically increased and decreased. The end of the hysteresis loop is marked by a frontier crisis of the chaotic attractor which also determines the end of the coupling strength interval where ES is possible.

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.

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.

Relationship of Topology, Multiscale Phase Synchronization, and State Transitions in Human Brain Networks

How the brain reconstitutes consciousness and cognition after a major perturbation like general anesthesia is an important question with significant neuroscientific and clinical implications. Recent empirical studies in animals and humans suggest that the recovery of consciousness after anesthesia is not random but ordered. Emergence patterns have been classified as progressive and abrupt transitions from anesthesia to consciousness, with associated differences in duration and electroencephalogram (EEG) properties. We hypothesized that the progressive and abrupt emergence patterns from the unconscious state are associated with, respectively, continuous and discontinuous synchronization transitions in functional brain networks. The discontinuous transition is explainable with the concept of explosive synchronization, which has been studied almost exclusively in network science. We used the Kuramato model, a simple oscillatory network model, to simulate progressive and abrupt transitions in anatomical human brain networks acquired from diffusion tensor imaging (DTI) of 82 brain regions. To facilitate explosive synchronization, distinct frequencies for hub nodes with a large frequency disassortativity (i.e., higher frequency nodes linking with lower frequency nodes, or vice versa) were applied to the brain network. In this simulation study, we demonstrated that both progressive and abrupt transitions follow distinct synchronization processes at the individual node, cluster, and global network levels. The characteristic synchronization patterns of brain regions that are “progressive and earlier” or “abrupt but delayed” account for previously reported behavioral responses of gradual and abrupt emergence from the unconscious state. The characteristic network synchronization processes observed at different scales provide new insights into how regional brain functions are reconstituted during progressive and abrupt emergence from the unconscious state. This theoretical approach also offers a principled explanation of how the brain reconstitutes consciousness and cognitive functions after physiologic (sleep), pharmacologic (anesthesia), and pathologic (coma) perturbations.

Coexistence of Quantized, Time Dependent, Clusters in Globally Coupled Oscillators

We report on a novel collective state, occurring in globally coupled nonidentical oscillators in the proximity of the point where the transition from the system's incoherent to coherent phase converts from explosive to continuous. In such a state, the oscillators form quantized clusters, where neither their phases nor their instantaneous frequencies are locked. The oscillators' instantaneous speeds are different within the clusters, but they form a characteristic cusped pattern and, more importantly, they behave periodically in time so that their average values are the same. Given its intrinsic specular nature with respect to the recently introduced Chimera states, the phase is termed the Bellerophon state. We provide an analytical and numerical description of Bellerophon states, and furnish practical hints on how to seek them in a variety of experimental and natural systems.

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.

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.

Functional and Topological Conditions for Explosive Synchronization Develop in Human Brain Networks with the Onset of Anesthetic-Induced Unconsciousness

Sleep, anesthesia, and coma share a number of neural features but the recovery profiles are radically different. To understand the mechanisms of reversibility of unconsciousness at the network level, we studied the conditions for gradual and abrupt transitions in conscious and anesthetized states. We hypothesized that the conditions for explosive synchronization (ES) in human brain networks would be present in the anesthetized brain just over the threshold of unconsciousness. To test this hypothesis, functional brain networks were constructed from multi-channel electroencephalogram (EEG) recordings in seven healthy subjects across conscious, unconscious, and recovery states. We analyzed four variables that are involved in facilitating ES in generic, non-biological networks: (1) correlation between node degree and frequency, (2) disassortativity (i.e., the tendency of highly-connected nodes to link with less-connected nodes, or vice versa), (3) frequency difference of coupled nodes, and (4) an inequality relationship between local and global network properties, which is referred to as the suppressive rule. We observed that the four network conditions for ES were satisfied in the unconscious state. Conditions for ES in the human brain suggest a potential mechanism for rapid recovery from the lightly-anesthetized state. This study demonstrates for the first time that the network conditions for ES, formerly shown in generic networks only, are present in empirically-derived functional brain networks. Further investigations with deep anesthesia, sleep, and coma could provide insight into the underlying causes of variability in recovery profiles of these unconscious states.

Effective centrality and explosive synchronization in complex networks

Synchronization of networked oscillators is known to depend fundamentally on the interplay between the dynamics of the graph's units and the microscopic arrangement of the network's structure. We here propose an effective network whose topological properties reflect the interplay between the topology and dynamics of the original network. On that basis, we are able to introduce the effective centrality, a measure that quantifies the role and importance of each network's node in the synchronization process. In particular, in the context of explosive synchronization, we use such a measure to assess the propensity of a graph to sustain an irreversible transition to synchronization. We furthermore discuss a strategy to induce the explosive behavior in a generic network, by acting only upon a fraction of its nodes.

Effects of random rewiring on the degree correlation of scale-free networks

Random rewiring is used to generate null networks for the purpose of analyzing the topological properties of scale-free networks, yet the effects of random rewiring on the degree correlation are subject to contradicting interpretations in the literature. We comprehensively analyze the degree correlation of randomly rewired scale-free networks and show that random rewiring increases disassortativity by reducing the average degree of the nearest neighbors of high-degree nodes. The effect can be captured by the measures of the degree correlation that consider all links in the network, but not by analogous measures that consider only links between degree peers, hence the potential for contradicting interpretations. We furthermore find that random and directional rewiring affect the topology of a scale-free network differently, even if the degree correlation of the rewired networks is the same. Consequently, the network dynamics is changed, which is proven here by means of the biased random walk.

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.

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.

Experimental evidence of explosive synchronization in mercury beating-heart oscillators

We report experimental evidence of explosive synchronization in coupled chemo-mechanical systems, namely in mercury beating-heart (MBH) oscillators. Connecting four MBH oscillators in a star network configuration and setting natural frequencies of each oscillator in proportion to the number of its links, a gradual increase of the coupling strength results in an abrupt and irreversible (first-order-like) transition from the system's unordered to ordered phase. On its turn, such a transition indicates the emergence of a bistable regime wherein coexisting states can be experimentally revealed. Finally, we prove how such a regime allows an experimental implementation of magneticlike states of synchronization, by the use of an external signal.

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.

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.

Explosive synchronization in adaptive and multilayer networks

At this time, explosive synchronization (ES) of networked oscillators is thought of as being rooted in the setting of specific microscopic correlation features between the natural frequencies of the oscillators and their effective coupling strengths. We show that ES is, in fact, far more general and can occur in adaptive and multilayer networks in the absence of such correlation properties. We first report evidence of ES for single-layer networks where a fraction f of the nodes have links adaptively controlled by a local order parameter, and we then extend the study to a variety of two-layer networks with a fraction f of their nodes coupled with each other by means of dependency links. In the latter case, we give evidence of ES regardless of the differences in the frequency distribution, in the topology of connections between the layers, or both. Finally, we provide a rigorous, analytical treatment to properly ground all of the observed scenarios and to advance the understanding of the actual mechanisms at the basis of ES in real-world systems.

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.

Explosive synchronization as a process of explosive percolation in dynamical phase space

Explosive synchronization and explosive percolation are currently two independent phenomena occurring in complex networks, where the former takes place in dynamical phase space while the latter in configuration space. It has been revealed that the mechanism of EP can be explained by the Achlioptas process, where the formation of a giant component is controlled by a suppressive rule. We here introduce an equivalent suppressive rule for ES. Before the critical point of ES, the suppressive rule induces the presence of multiple, small sized, synchronized clusters, while inducing the abrupt formation of a giant cluster of synchronized oscillators at the critical coupling strength. We also show how the explosive character of ES degrades into a second-order phase transition when the suppressive rule is broken. These results suggest that our suppressive rule can be considered as a dynamical counterpart of the Achlioptas process, indicating that ES and EP can be unified into a same framework.

Disorder induces explosive synchronization

We study explosive synchronization, a phenomenon characterized by first-order phase transitions between incoherent and synchronized states in networks of coupled oscillators. While explosive synchronization has been the subject of many recent studies, in each case strong conditions on the heterogeneity of the network, its link weights, or its initial construction are imposed to engineer a first-order phase transition. This raises the question of how robust explosive synchronization is in view of more realistic structural and dynamical properties. Here we show that explosive synchronization can be induced in mildly heterogeneous networks by the addition of quenched disorder to the oscillators' frequencies, demonstrating that it is not only robust to, but moreover promoted by, this natural mechanism. We support these findings with numerical and analytical results, presenting simulations of a real neural network as well as a self-consistency theory used to study synthetic networks.

Basin of attraction determines hysteresis in explosive synchronization

Spontaneous explosive emergent behavior takes place in heterogeneous networks when the frequencies of the nodes are positively correlated to the node degree. A central feature of such explosive transitions is a hysteretic behavior at the transition to synchronization. We unravel the underlying mechanisms and show that the dynamical origin of the hysteresis is a change of basin of attraction of the synchronization state. Our findings hold for heterogeneous networks with star graph motifs such as scale-free networks, and hence, reveal how microscopic network parameters such as node degree and frequency affect the global network properties and can be used for network design and control.

Criterion for the emergence of explosive synchronization transitions in networks of phase oscillators

The emergence of explosive synchronization transitions in networks of phase oscillators recently has become one of the most interesting topics. It is widely believed that the large frequency mismatch of a pair of oscillators (also known as disassortativity in frequency) is a direct cause of an explosive synchronization. It is found that, besides the disassortativity in frequency, the disassortativity in node degree also shows up in connection with the first-order synchronization transition. In this paper, we simulate the Kuramoto model on top of a family of networks with different degree-degree and frequency-frequency correlation patterns. Results show that only when the degrees and natural frequencies of the network's nodes are both disassortative can an explosive synchronization occur.