Explosive synchronization in a general complex network

Higher-order interactions induce stepwise explosive phase transitions

In recent studies, it has been established that higher-order interactions in coupled oscillators can induce a process from continuous to explosive phase transition. In this study, we identify a phase transition, termed the stepwise explosive phase transition, characterized by the emergence of multiple critical phase plateaus in a globally frequency-weighted coupled pendulum model. This transition bridges the continuous and explosive phase transitions, arising from a delicate balance between attractive higher-order interactions and repulsive pairwise interactions. Specifically, the stepwise explosive phase transition occurs when the higher-order coupling is moderate, neither large nor small, while the pairwise interactions remain repulsive. Our analysis shows that stronger attractive higher-order interactions necessitate weaker repulsive pairwise interactions, leading to partial frequency locking among oscillators and triggering the stepwise transition. We construct an analytical framework using self-consistent equations to provide an approximation of the steady-state behavior. This study uncovers an alternative pathway to the desynchronization, and it provides additional insights into the phase transitions in coupled dynamical networks.

Synchronization transitions in adaptive simplicial complexes with cooperative and competitive dynamics

Adaptive network is a powerful presentation to describe different real-world phenomena. However, current models often neglect higher-order interactions (beyond pairwise interactions) and diverse adaptation types (cooperative and competitive) commonly observed in systems such as the human brain and social networks. This work addresses this gap by incorporating these factors into a model that explores their impact on collective properties such as synchronization. Through simplified network representations, we investigate how the simultaneous presence of cooperative and competitive adaptations influences phase transitions. Our findings reveal a transition from first-order to second-order synchronization as the strength of higher-order interactions increases under competitive adaptation. We also demonstrate the possibility of synchronization even without pairwise interactions, provided there is strong enough higher-order coupling. When only competitive adaptations are present, the system exhibits second-order-like phase transitions and clustering. Conversely, with a combination of cooperative and competitive adaptations, the system undergoes a first-order-like phase transition, characterized by a sharp transition to the synchronized state without reverting to an incoherent state during backward transitions. The specific nature of these second-order-like transitions varies depending on the coupling strengths and mean degrees. With our model, we can control not only when the system synchronizes but also the way the system goes to synchronization.

Synchronization transitions in phase oscillator populations with partial adaptive coupling

The adaptation underlying many realistic processes plays a pivotal role in shaping the collective dynamics of diverse systems. Here, we untangle the generic conditions for synchronization transitions in a system of coupled phase oscillators incorporating the adaptive scheme encoded by the feedback between the coupling and the order parameter via a power-law function with different weights. We mathematically argue that, in the subcritical and supercritical correlation scenarios, there exists no critical adaptive fraction for synchronization transitions converting from the first (second)-order to the second (first)-order. In contrast to the synchronization transitions previously deemed, the explosive and continuous phase transitions take place in the corresponding regions as long as the adaptive fraction is nonzero, respectively. Nevertheless, we uncover that, at the critical correlation, the routes toward synchronization depend crucially on the relative adaptive weights. In particular, we unveil that the emergence of a range of interrelated scaling behaviors of the order parameter near criticality, manifesting the subcritical and supercritical bifurcations, are responsible for various observed phase transitions. Our work, thus, provides profound insights for understanding the dynamical nature of phase transitions, and for better controlling and manipulating synchronization transitions in networked systems with adaptation.

Explosive transitions in coupled Lorenz oscillators

We study the transition to synchronization in an ensemble of chaotic oscillators that are interacting on a star network. These oscillators possess an invariant symmetry and we study emergent behavior by introducing the timescale variations in the dynamics of the nodes and the hub. If the coupling preserves the symmetry, the ensemble exhibits consecutive explosive transitions, each one associated with a hysteresis. The first transition is the explosive synchronization from a desynchronized state to a synchronized state which occurs discontinuously with the formation of intermediate clusters. These clusters appear because of the driving-induced multistability and the resulting attractors exhibit intermittent synchrony (antisynchrony). The second transition is the explosive death that occurs as a result of stabilization of the stable fixed points. However, if the symmetry is not preserved, the system again makes a first-order transition from an oscillatory state to death, namely, an explosive death. These transitions are studied with the help of the master stability functions, Lyapunov exponents, and the stability analysis.

Solvable dynamics of the three-dimensional Kuramoto model with frequency-weighted coupling

The presence of coupling heterogeneity is deemed to be a natural attribute in realistic systems comprised of many interacting agents. In this work, we study dynamics of the 3D Kuramoto model with heterogeneous couplings, where the strength of coupling for each agent is weighted by its intrinsic rotation frequency. The critical coupling strength for the instability of incoherence is rigorously derived to be zero by carrying out a linear stability analysis of an incoherent state. For positive values of the coupling strength, at which the incoherence turns out to be unstable, a self-consistency approach is developed to theoretically predict the degree of global coherence of the model. Our theoretical analyses match well with numerical simulations, which helps us to deepen the understanding of collective behaviors spontaneously emerged in heterogeneously coupled high-dimensional dynamical networks.

Discontinuous phase transition switching induced by a power-law function between dynamical parameters in coupled oscillators

In biological or physical systems, the intrinsic properties of oscillators are heterogeneous and correlated. These two characteristics have been empirically validated and have garnered attention in theoretical studies. In this paper, we propose a power-law function existed between the dynamical parameters of the coupled oscillators, which can control discontinuous phase transition switching. Unlike the special designs for the coupling terms, this generalized function within the dynamical term reveals another path for generating the first-order phase transitions. The power-law relationship between dynamic characteristics is reasonable, as observed in empirical studies, such as long-term tremor activity during volcanic eruptions and ion channel characteristics of the Xenopus expression system. Our work expands the conditions that used to be strict for the occurrence of the first-order phase transitions and deepens our understanding of the impact of correlation between intrinsic parameters on phase transitions. We explain the reason why the continuous phase transition switches to the discontinuous phase transition when the control parameter is at a critical value.

Order parameter dynamics in complex systems: From models to data

Collective ordering behaviors are typical macroscopic manifestations embedded in complex systems and can be ubiquitously observed across various physical backgrounds. Elements in complex systems may self-organize via mutual or external couplings to achieve diverse spatiotemporal coordinations. The order parameter, as a powerful quantity in describing the transition to collective states, may emerge spontaneously from large numbers of degrees of freedom through competitions. In this minireview, we extensively discussed the collective dynamics of complex systems from the viewpoint of order-parameter dynamics. A synergetic theory is adopted as the foundation of order-parameter dynamics, and it focuses on the self-organization and collective behaviors of complex systems. At the onset of macroscopic transitions, slow modes are distinguished from fast modes and act as order parameters, whose evolution can be established in terms of the slaving principle. We explore order-parameter dynamics in both model-based and data-based scenarios. For situations where microscopic dynamics modeling is available, as prototype examples, synchronization of coupled phase oscillators, chimera states, and neuron network dynamics are analytically studied, and the order-parameter dynamics is constructed in terms of reduction procedures such as the Ott-Antonsen ansatz, the Lorentz ansatz, and so on. For complicated systems highly challenging to be well modeled, we proposed the eigen-microstate approach (EMP) to reconstruct the macroscopic order-parameter dynamics, where the spatiotemporal evolution brought by big data can be well decomposed into eigenmodes, and the macroscopic collective behavior can be traced by Bose-Einstein condensation-like transitions and the emergence of dominant eigenmodes. The EMP is successfully applied to some typical examples, such as phase transitions in the Ising model, climate dynamics in earth systems, fluctuation patterns in stock markets, and collective motion in living systems.

Deterministic correlations enhance synchronization in oscillator populations with heterogeneous coupling

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

Synchronization transitions in Kuramoto networks with higher-mode interaction

Synchronization is an omnipresent collective phenomenon in nature and technology, whose understanding is still elusive for real-world systems in particular. We study the synchronization transition in a phase oscillator system with two nonvanishing Fourier-modes in the interaction function, hence going beyond the Kuramoto paradigm. We show that the transition scenarios crucially depend on the interplay of the two coupling modes. We describe the multistability induced by the presence of a second coupling mode. By extending the collective coordinate approach, we describe the emergence of various states observed in the transition from incoherence to coherence. Remarkably, our analysis suggests that, in essence, the two-mode coupling gives rise to states characterized by two independent but interacting groups of oscillators. We believe that these findings will stimulate future research on dynamical systems, including complex interaction functions beyond the Kuramoto-type.

Sensitive dynamics of brain cognitive networks and its resource constraints

It is well known that brain functions are closely related to the synchronization of brain networks, but the underlying mechanisms are still not completely understood. To study this problem, we here focus on the synchronization of cognitive networks, in contrast to that of a global brain network, as individual brain functions are in fact performed by different cognitive networks but not the global network. In detail, we consider four different levels of brain networks and two approaches, i.e., either with or without resource constraints. For the case of without resource constraints, we find that global brain networks have fundamentally different behaviors from that of the cognitive networks; i.e., the former has a continuous synchronization transition, while the latter shows a novel transition of oscillatory synchronization. This feature of oscillation comes from the sparse links among the communities of cognitive networks, resulting in coupling sensitive dynamics of brain cognitive networks. While for the case of resource constraints, we find that at the global level, the synchronization transition becomes explosive, in contrast to the continuous synchronization for the case of without resource constraints. At the level of cognitive networks, the transition also becomes explosive and the coupling sensitivity is significantly reduced, thus guaranteeing the robustness and fast switch of brain functions. Moreover, a brief theoretical analysis is provided.

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.

Generalized frustration in the multidimensional Kuramoto model

The Kuramoto model describes how coupled oscillators synchronize their phases as the intensity of the coupling increases past a threshold. The model was recently extended by reinterpreting the oscillators as particles moving on the surface of unit spheres in a D-dimensional space. Each particle is then represented by a D-dimensional unit vector; for D=2 the particles move on the unit circle and the vectors can be described by a single phase, recovering the original Kuramoto model. This multidimensional description can be further extended by promoting the coupling constant between the particles to a matrix K that acts on the unit vectors. As the coupling matrix changes the direction of the vectors, it can be interpreted as a generalized frustration that tends to hinder synchronization. In a recent paper we studied in detail the role of the coupling matrix for D=2. Here we extend this analysis to arbitrary dimensions. We show that, for identical particles, when the natural frequencies are set to zero, the system converges either to a stationary synchronized state, given by one of the real eigenvectors of K, or to an effective two-dimensional rotation, defined by one of the complex eigenvectors of K. The stability of these states depends on the set eigenvalues and eigenvectors of the coupling matrix, which controls the asymptotic behavior of the system, and therefore, can be used to manipulate these states. When the natural frequencies are not zero, synchronization depends on whether D is even or odd. In even dimensions the transition to synchronization is continuous and rotating states are replaced by active states, where the module of the order parameter oscillates while it rotates. If D is odd the phase transition is discontinuous and active states can be suppressed for some distributions of natural frequencies.

Heterogeneous Nucleation in Finite-Size Adaptive Dynamical Networks

Phase transitions in equilibrium and nonequilibrium systems play a major role in the natural sciences. In dynamical networks, phase transitions organize qualitative changes in the collective behavior of coupled dynamical units. Adaptive dynamical networks feature a connectivity structure that changes over time, coevolving with the nodes' dynamical state. In this Letter, we show the emergence of two distinct first-order nonequilibrium phase transitions in a finite-size adaptive network of heterogeneous phase oscillators. Depending on the nature of defects in the internal frequency distribution, we observe either an abrupt single-step transition to full synchronization or a more gradual multistep transition. This observation has a striking resemblance to heterogeneous nucleation. We develop a mean-field approach to study the interplay between adaptivity and nodal heterogeneity and describe the dynamics of multicluster states and their role in determining the character of the phase transition. Our work provides a theoretical framework for studying the interplay between adaptivity and nodal heterogeneity.

Relaxation dynamics of phase oscillators with generic heterogeneous coupling

The coupled phase oscillator model serves as a paradigm that has been successfully used to shed light on the collective dynamics occurring in large ensembles of interacting units. It was widely known that the system experiences a continuous (second-order) phase transition to synchronization by gradually increasing the homogeneous coupling among the oscillators. As the interest in exploring synchronized dynamics continues to grow, the heterogeneous patterns between phase oscillators have received ample attention during the past years. Here, we consider a variant of the Kuramoto model with quenched disorder in their natural frequencies and coupling. Correlating these two types of heterogeneity via a generic weighted function, we systematically investigate the impacts of the heterogeneous strategies, the correlation function, and the natural frequency distribution on the emergent dynamics. Importantly, we develop an analytical treatment for capturing the essential dynamical properties of the equilibrium states. In particular, we uncover that the critical threshold corresponding to the onset of synchronization is unaffected by the location of the inhomogeneity, which, however, does depend crucially on the value of the correlation function at its center. Furthermore, we reveal that the relaxation dynamics of the incoherent state featuring the responses to external perturbations is significantly shaped by all the considered effects, thereby leading to various decaying mechanisms of the order parameters in the subcritical region. Moreover, we untangle that synchronization is facilitated by the out-coupling strategy in the supercritical region. Our study is a step forward in highlighting the potential importance of the inhomogeneous patterns involved in the complex systems, and could thus provide theoretical insights for profoundly understanding the generic statistical mechanical properties of the steady states toward synchronization.

Route to synchronization in coupled phase oscillators with frequency-dependent coupling: Explosive or continuous?

Interconnected dynamical systems often transition between states of incoherence and synchronization due to changes in system parameters. These transitions could be continuous (gradual) or explosive (sudden) and may result in failures, which makes determining their nature important. In this study, we abstract dynamical networks as an ensemble of globally coupled Kuramoto-like phase oscillators with frequency-dependent coupling and investigate the mechanisms for transition between incoherent and synchronized dynamics. The characteristics that dictate a continuous or explosive route to synchronization are the distribution of the natural frequencies of the oscillators, quantified by a probability density function g(ω), and the relation between the coupling strength and natural frequency of an oscillator, defined by a frequency-coupling strength correlation function f(ω). Our main results are conditions on f(ω) and g(ω) that result in continuous or explosive routes to synchronization and explain the underlying physics. The analytical developments are validated through numerical examples.

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.

Explosive transitions to synchronization in networks of frequency dipoles

We reveal that an introduction of frequency-weighted inter-layer coupling term in networks of frequency dipoles can induce explosive synchronization transitions. The reason for explosive synchronization is that the oscillators with synchronization superiority are moderately suppressed. The findings show that a super-linear correlation induces explosive synchronization in networks of frequency dipoles, while a linear or sub-linear correlation excites chimera-like states. Clearly, the synchronization transition mode of networks of frequency dipoles is controlled by the power-law exponent. In addition, by means of the mean-field approximation, we obtain the critical values of the coupling strength within and between layers in two limit cases. The results of theoretical analysis are in good agreement with those of numerical simulation. Compared with the previous models, the model proposed in this paper retains the topological structure of network and the intrinsic properties of oscillators, so it is easy to realize pinning control.

Partial locking in phase-oscillator populations with heterogenous coupling

We consider a variant of the mean-field model of coupled phase oscillators with uniform distribution of natural frequencies. By establishing correlations between the quenched disorder of intrinsic frequencies and coupling strength with both in- and out-coupling heterogeneities, we reveal a generic criterion for the onset of partial locking that takes place in a domain with the coexistence of phase-locked oscillators and drifters. The critical points manifesting the instability of the stationary states are obtained analytically. In particular, the bifurcation mechanism of the equilibrium states is uncovered by the use of frequency-dependent version of the Ott-Antonsen reduction consistently with the analysis based on the self-consistent approach. We demonstrate that both the manner of coupling heterogeneity and correlation exponent have influence on the emergent patterns of partial locking. Our research could find applicability in better understanding the phase transitions and related collective phenomena involving synchronization control in networked systems.

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.

Ageing transitions in a network of Rulkov neurons

The phenomenon of ageing transitions (AT) in a Erdős–Rényi network of coupled Rulkov neurons is studied with respect to parameters modelling network connectivity, coupling strength and the fractional ratio of inactive neurons in the network. A general mean field coupling is proposed to model the neuronal interactions. A standard order parameter is defined for quantifying the network dynamics. Investigations are undertaken for both the noise free network as well as stochastic networks, where the interneuronal coupling strength is assumed to be superimposed with additive noise. The existence of both smooth and explosive AT are observed in the parameter space for both the noise free and the stochastic networks. The effects of noise on AT are investigated and are found to play a constructive role in mitigating the effects of inactive neurons and reducing the parameter regime in which explosive AT is observed.

Effects of topological characteristics on rhythmic states of the D-dimensional Kuramoto model in complex networks

Synchronization is a ubiquitous phenomenon in engineering and natural ecosystems. While the dynamics of synchronization modeled by the Kuramoto model are commonly studied in two dimensions and the state of dynamic units is characterized by a scalar angle variable, we studied the Kuramoto model generalized to D dimensions in the framework of a complex network and utilized the local synchronous order parameter between the agent and its neighbors as the controllable variable to adjust the coupling strength. Here, we reported that average connectivity of networks affects the time-dependent, rhythmic, cyclic state. Importantly, we found that the level of heterogeneity of networks governs the rhythmic state in the transition process. The analytical treatment for observed scenarios in a D-dimensional Kuramoto model at D=3 was provided. These results offered a platform for a better understanding of time-dependent swarming and flocking dynamics in nature.

Catalytic feed-forward explosive synchronization in multilayer networks

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

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.

Explosive synchronization in interlayer phase-shifted Kuramoto oscillators on multiplex networks

Different methods have been proposed in the past few years to incite explosive synchronization (ES) in Kuramoto phase oscillators. In this work, we show that the introduction of a phase shift α in interlayer coupling terms of a two-layer multiplex network of Kuramoto oscillators can also instigate ES in the layers. As α→π/2, ES emerges along with hysteresis. The width of hysteresis depends on the phase shift α, interlayer coupling strength, and natural frequency mismatch between mirror nodes. A mean-field analysis is performed to justify the numerical results. Similar to earlier works, the suppression of synchronization is accountable for the occurrence of ES. The robustness of ES against changes in network topology and natural frequency distribution is tested. Finally, taking a suggestion from the synchronized state of the multiplex networks, we extend the results to classical single networks where some specific links are assigned phase-shifted interactions.

Transition to synchronization in heterogeneous inhibitory neural networks with structured synapses

Inhibitory neurons form an extensive network involved in the development of different rhythms in the cerebral cortex. A transition from an incoherent state, where all inhibitory neurons fire unrelated to each other, to a synchronized or locked state, where all or most neurons define a tight firing pattern, is maybe the most salient process to analyze when considering neuronal rhythms. In this work, we analyzed whether different patterns of effective synaptic connectivity may support a first-order-like transition in this path to synchronization. Such an "explosive" phenomenon may be relevant in neural processes, as normal cognitive processing in different tasks and some neurological disorders manifest an increased power in many neuronal rhythms, supported by an extended concerted spiking activity and an abrupt change to this state. Furthermore, we built an adaptive mechanism that supports the generation of this kind of network, which rapidly creates the underlying structure based on the ongoing firing statistics.

Discontinuous phase transition in the Kuramoto model with asymmetric dynamic interaction

We investigate the critical behavior of the modified Kuramoto model with an asymmetric dynamic interaction which has been proposed to explain the difference between the synchronized frequency and the average intrinsic frequency. We find that the discontinuous phase transition arises when oscillators interact only with other oscillators whose phases are ahead. From the comparison with the conventional Kuramoto model in which the interaction possesses phase reflection symmetry, we conclude that the dynamical symmetry breaking and the dynamic change in interaction structure play important roles in changing the transition nature from continuous to discontinuous ones.

Universal scaling and phase transitions of coupled phase oscillator populations

The Kuramoto model, which serves as a paradigm for investigating synchronization phenomena of oscillatory systems, is known to exhibit second-order, i.e., continuous, phase transitions in the macroscopic order parameter. Here we generalize a number of classical results by presenting a general framework for analytically capturing the critical scaling properties of the order parameter at the onset of synchronization. Using a self-consistent approach and constructing a characteristic function, we identify various phase transitions toward synchrony and establish scaling relations describing the asymptotic dependence of the order parameter on the coupling strength near the critical point. We find that the geometric properties of the characteristic function, which depends on the natural frequency distribution, determine the scaling properties of order parameter above criticality.

Effective-potential approach to hybrid synchronization transitions

The Kuramoto model exhibits different types of synchronization transitions depending on the type of natural frequency distribution. To obtain these results, the Kuramoto self-consistency equation (SCE) approach has been used successfully. However, this approach affords only limited understanding of more detailed properties such as the stability. Here we extend the SCE approach by introducing an effective potential, that is, an integral version of the SCE. We examine the landscape of this effective potential for second-order, first-order, and hybrid synchronization transitions in the thermodynamic limit. In particular, for the hybrid transition, we find that the minimum of effective potential displays a plateau across the region in which the order parameter jumps. This result suggests that the effective potential can be used to determine a type of synchronization transition.

Asymmetric dynamic interaction shifts synchronized frequency of coupled oscillators

Interacting dynamic agents can often exhibit synchronization. It has been reported that the rhythm all agents agree on in the synchronized state could be different from the average of intrinsic rhythms of individual agents. Hinted by such a real-world behavior of the interaction-driven change of the average phase velocity, we propose a modified version of the Kuramoto model, in which the ith oscillator of the phase ϕ_i interacts with other oscillator j only when the phase difference \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\phi }}_{{j}}$$\end{document}ϕj − \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\phi }}_{i}$$\end{document}ϕi is in a limited range [−βπ, απ]. From extensive numerical investigations, we conclude that the asymmetric dynamic interaction characterized by β ≠ α leads to the shift of the synchronized frequency with respect to the original distribution of the intrinsic frequency. We also perform and report our computer-based synchronization experiment, which exhibits the expected shift of the synchronized frequency of human participants. In analogy to interacting runners, our result roughly suggests that agents tend to run faster if they are more influenced by runners ahead than behind. We verify the observation by using a simple model of interacting runners.

Cluster synchronization: From single-layer to multi-layer networks

Cluster synchronization is a very common phenomenon occurring in single-layer complex networks, and it can also be observed in many multilayer networks in real life. In this paper, we study cluster synchronization of an isolated network and then focus on that of the network when it is influenced by an external network. We mainly explore how the influence layer impacts the cluster synchronization of the interest layer in a multilayer network. Considering that the clusters are changeable, we introduce a term called "cluster synchronizability" to measure the ability of a network to reach cluster synchronization. Since cluster synchronizability is intimately associated with the structure of the coupled external layer, we consider community networks and networks with different densities as the coupled layer. Besides the topology structure, the connection between two layers may also have an influence on the cluster synchronization of the interest layer. We study three different patterns of connection, including typical positive correlation, negative correlation, and random correlation and find that they all have a certain influence. However, the general theoretical analysis of cluster synchronization on multilayer networks is still a challenging topic. In this paper, we mainly use numerical simulations to discuss cluster synchronization.

Inhibition-induced explosive synchronization in multiplex networks

To date, explosive synchronization (ES) in a network is shown to be originated from considering either degree-frequency correlation, frequency-coupling strength correlation, inertia, or adaptively controlled phase oscillators. Here we show that ES is a generic phenomenon and can occur in any network by multiplexing it with an appropriate layer without even considering any prerequisite for the emergence of ES. We devise a technique which leads to the occurrence of ES with hysteresis loop in a network upon its multiplexing with a negatively coupled (or inhibitory) layer. The impact of various structural properties of positively coupled (or excitatory) and inhibitory layers along with the strength of multiplexing in gaining control over the induced ES transition is discussed. Analytical prediction for the spread of phase distribution of each layer is provided, which is in good agreement with the numerical assessment. This investigation is a step forward in highlighting the importance of multiplex framework not only in bringing phenomena which are not possible in an isolated network but also in providing more structural control over the induced phenomena.

Explosive synchronization in frequency displaced multiplex networks

Motivated by the recent multiplex framework of complex networks, in this work, we investigate if explosive synchronization can be induced in the multiplex network of two layers. Using nonidentical Kuramoto oscillators, we show that a sufficient frequency mismatch between two layers of a multiplex network can lead to explosive inter- and intralayer synchronization due to mutual frustration in the completion of the synchronization processes of the layers, generating a hybrid transition without imposing any specific structure-dynamics correlation.

Critical neuromorphic computing based on explosive synchronization

Synchronous oscillations in neuronal ensembles have been proposed to provide a neural basis for the information processes in the brain. In this work, we present a neuromorphic computing algorithm based on oscillator synchronization in a critical regime. The algorithm uses the high-dimensional transient dynamics perturbed by an input and translates it into proper output stream. One of the benefits of adopting coupled phase oscillators as neuromorphic elements is that the synchrony among oscillators can be finely tuned at a critical state. Especially near a critical state, the marginally synchronized oscillators operate with high efficiency and maintain better computing performances. We also show that explosive synchronization that is induced from specific neuronal connectivity produces more improved and stable outputs. This work provides a systematic way to encode computing in a large size coupled oscillator, which may be useful in designing neuromorphic devices.

Synchronization transition in Sakaguchi-Kuramoto model on complex networks with partial degree-frequency correlation

We investigate transition to synchronization in the Sakaguchi-Kuramoto (SK) model on complex networks analytically as well as numerically. Natural frequencies of a percentage (f) of higher degree nodes of the network are assumed to be correlated with their degrees and that of the remaining nodes are drawn from some standard distribution, namely, Lorentz distribution. The effects of variation of f and phase frustration parameter α on transition to synchronization are investigated in detail. Self-consistent equations involving critical coupling strength (λ_c) and group angular velocity (Ω_c) at the onset of synchronization have been derived analytically in the thermodynamic limit. For the detailed investigation, we considered the SK model on scale-free (SF) as well as Erdős-Rényi (ER) networks. Interestingly, explosive synchronization (ES) has been observed in both networks for different ranges of values of α and f. For SF networks, as the value of f is set within 10%≤f≤70%, the range of the values of α for existence of the ES is greatly enhanced compared to the fully degree-frequency correlated case when scaling exponent γ<3. ES is also observed in SF networks with γ>3, which is never observed in fully degree-frequency correlated environment. On the other hand, for random networks, ES observed is in a narrow window of α when the value of f is taken within 30%≤f≤50%. In all the cases, critical coupling strengths for transition to synchronization computed from the analytically derived self-consistent equations show a very good agreement with the numerical results. Finally, we observe ES in the metabolic network of the roundworm Caenorhabditis elegans in partially degree-frequency correlated environment.

Exact explosive synchronization transitions in Kuramoto oscillators with time-delayed coupling

Synchronization commonly occurs in many natural and man-made systems, from neurons in the brain to cardiac cells to power grids to Josephson junction arrays. Transitions to or out of synchrony for coupled oscillators depend on several factors, such as individual frequencies, coupling, interaction time delays and network structure-function relation. Here, using a generalized Kuramoto model of time-delay coupled phase oscillators with frequency-weighted coupling, we study the stability of incoherent and coherent states and the transitions to or out of explosive (abrupt, first-order like) phase synchronization. We analytically derive the exact formulas for the critical coupling strengths at different time delays in both directions of increasing (forward) and decreasing (backward) coupling strengths. We find that time-delay does not affect the transition for the backward direction but can shift the transition for the forward direction of increasing coupling strength. These results provide valuable insights into our understanding of dynamical mechanisms for explosive synchronization in presence of often unavoidable time delays present in many physical and biological systems.

Mechanisms of hysteresis in human brain networks during transitions of consciousness and unconsciousness: Theoretical principles and empirical evidence

Hysteresis, the discrepancy in forward and reverse pathways of state transitions, is observed during changing levels of consciousness. Identifying the underlying mechanism of hysteresis phenomena in the brain will enhance the ability to understand, monitor, and control state transitions related to consciousness. We hypothesized that hysteresis in brain networks shares the same underlying mechanism of hysteresis as other biological and non-biological networks. In particular, we hypothesized that the principle of explosive synchronization, which can mediate abrupt state transitions, would be critical to explaining hysteresis in the brain during conscious state transitions. We analyzed high-density electroencephalogram (EEG) that was acquired in healthy human volunteers during conscious state transitions induced by the general anesthetics sevoflurane or ketamine. We developed a novel method to monitor the temporal evolution of EEG networks in a parameter space, which consists of the strength and topography of EEG-based networks. Furthermore, we studied conditions of explosive synchronization in anatomically informed human brain network models. We identified hysteresis in the trajectory of functional brain networks during state transitions. The model study and empirical data analysis explained various hysteresis phenomena during the loss and recovery of consciousness in a principled way: (1) more potent anesthetics induce a larger hysteresis; (2) a larger range of EEG frequencies facilitates transitions into unconsciousness and impedes the return of consciousness; (3) hysteresis of connectivity is larger than that of EEG power; and (4) the structure and strength of functional brain networks reconfigure differently during the loss vs. recovery of consciousness. We conclude that the hysteresis phenomena observed during the loss and recovery of consciousness are generic network features. Furthermore, the state transitions are grounded in the same principle as state transitions in complex non-biological networks, especially during perturbation. These findings suggest the possibility of predicting and modulating hysteresis of conscious state transitions in large-scale brain networks.

Hysteresis, characterized by distinct forward and reverse phase transitions, is ubiquitous in nature. For example, there are distinct temperatures for water freezing and ice melting. Similarly, it has been found that state transitions related to consciousness exhibit hysteresis. In particular, the concentration of general anesthetics required to achieve loss of consciousness is significantly higher than the concentration at which consciousness is regained. However, it is unknown whether this is trivially reducible to the pharmacology of these drugs or if it is something related to brain function itself. In this study, we took a novel, network-based approach and hypothesized that the hysteresis observed during anesthetic state transitions shares the same underlying mechanism as that observed in non-biological networks. Our computational modeling, analytic study, and high-density human EEG analysis suggest that various hysteresis phenomena during loss and recovery of consciousness can be explained in principled ways by generic network features. Identifying these network mechanisms of hysteresis in the brain also provides a unified framework for understanding the radically different conscious state transitions associated with sleep, anesthesia, and disorders of consciousness.

Rhythmic synchronization and hybrid collective states of globally coupled oscillators

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

Enhancing network synchronization by phase modulation

Due to time delays in signal transmission and processing, phase lags are inevitable in realistic complex oscillator networks. Conventional wisdom is that phase lags are detrimental to network synchronization. Here we show that judiciously chosen phase lag modulations can result in significantly enhanced network synchronization. We justify our strategy of phase modulation, demonstrate its power in facilitating and enhancing network synchronization with synthetic and empirical network models, and provide an analytic understanding of the underlying mechanism. Our work provides an alternative approach to synchronization optimization in complex networks, with insights into control of complex nonlinear networks.

Phase definition to assess synchronization quality of nonlinear oscillators

This paper proposes a phase definition, named the vector field phase, which can be defined for systems with arbitrary finite dimension and is a monotonically increasing function of time. The proposed definition can properly quantify the dynamics in the flow direction, often associated with the null Lyapunov exponent. Numerical examples that use benchmark periodic and chaotic oscillators are discussed to illustrate some of the main features of the definition, which are that (i) phase information can be obtained either from the vector field or from a time series, (ii) it permits not only detection of phase synchronization but also quantification of it, and (iii) it can be used in the phase synchronization of very different oscillators.

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.

How stimulation frequency and intensity impact on the long-lasting effects of coordinated reset stimulation

Several brain diseases are characterized by abnormally strong neuronal synchrony. Coordinated Reset (CR) stimulation was computationally designed to specifically counteract abnormal neuronal synchronization processes by desynchronization. In the presence of spike-timing-dependent plasticity (STDP) this may lead to a decrease of synaptic excitatory weights and ultimately to an anti-kindling, i.e. unlearning of abnormal synaptic connectivity and abnormal neuronal synchrony. The long-lasting desynchronizing impact of CR stimulation has been verified in pre-clinical and clinical proof of concept studies. However, as yet it is unclear how to optimally choose the CR stimulation frequency, i.e. the repetition rate at which the CR stimuli are delivered. This work presents the first computational study on the dependence of the acute and long-term outcome on the CR stimulation frequency in neuronal networks with STDP. For this purpose, CR stimulation was applied with Rapidly Varying Sequences (RVS) as well as with Slowly Varying Sequences (SVS) in a wide range of stimulation frequencies and intensities. Our findings demonstrate that acute desynchronization, achieved during stimulation, does not necessarily lead to long-term desynchronization after cessation of stimulation. By comparing the long-term effects of the two different CR protocols, the RVS CR stimulation turned out to be more robust against variations of the stimulation frequency. However, SVS CR stimulation can obtain stronger anti-kindling effects. We revealed specific parameter ranges that are favorable for long-term desynchronization. For instance, RVS CR stimulation at weak intensities and with stimulation frequencies in the range of the neuronal firing rates turned out to be effective and robust, in particular, if no closed loop adaptation of stimulation parameters is (technically) available. From a clinical standpoint, this may be relevant in the context of both invasive as well as non-invasive CR stimulation.

Bounding the first exit from the basin: Independence times and finite-time basin stability

We study the stability of deterministic systems, given sequences of large, jump-like perturbations. Our main result is the derivation of a lower bound for the probability of the system to remain in the basin, given that perturbations are rare enough. This bound is efficient to evaluate numerically. To quantify rare enough, we define the notion of the independence time of such a system. This is the time after which a perturbed state has probably returned close to the attractor, meaning that subsequent perturbations can be considered separately. The effect of jump-like perturbations that occur at least the independence time apart is thus well described by a fixed probability to exit the basin at each jump, allowing us to obtain the bound. To determine the independence time, we introduce the concept of finite-time basin stability, which corresponds to the probability that a perturbed trajectory returns to an attractor within a given time. The independence time can then be determined as the time scale at which the finite-time basin stability reaches its asymptotic value. Besides that, finite-time basin stability is a novel probabilistic stability measure on its own, with potential broad applications in complex systems.

The choroid plexus is an important circadian clock component

Mammalian circadian clocks have a hierarchical organization, governed by the suprachiasmatic nucleus (SCN) in the hypothalamus. The brain itself contains multiple loci that maintain autonomous circadian rhythmicity, but the contribution of the non-SCN clocks to this hierarchy remains unclear. We examine circadian oscillations of clock gene expression in various brain loci and discovered that in mouse, robust, higher amplitude, relatively faster oscillations occur in the choroid plexus (CP) compared to the SCN. Our computational analysis and modeling show that the CP achieves these properties by synchronization of “twist” circadian oscillators via gap-junctional connections. Using an in vitro tissue coculture model and in vivo targeted deletion of the Bmal1 gene to silence the CP circadian clock, we demonstrate that the CP clock adjusts the SCN clock likely via circulation of cerebrospinal fluid, thus finely tuning behavioral circadian rhythms.

The suprachiasmatic nucleus (SCN) has been thought of as the master circadian clock, but peripheral circadian clocks do exist. Here, the authors show that the choroid plexus displays oscillations more robust than the SCN and that can be described as a Poincaré oscillator with negative twist.

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.

Partial inertia induces additional phase transition in the majority vote model

Explosive (i.e., discontinuous) transitions have aroused great interest by manifesting in distinct systems, such as synchronization in coupled oscillators, percolation regime, absorbing phase transitions, and more recently, the majority-vote model with inertia. In the latter, the model rules are slightly modified by the inclusion of a term depending on the local spin (an inertial term). In such a case, Chen et al. [Phys Rev. E 95, 042304 (2017)2470-004510.1103/PhysRevE.95.042304] have found that relevant inertia changes the nature of the phase transition in complex networks, from continuous to discontinuous. Here we give a further step by embedding inertia only in vertices with degree larger than a threshold value 〈k〉k^{*}, 〈k〉 being the mean system degree and k^{*} the fraction restriction. Our results, from mean-field analysis and extensive numerical simulations, reveal that an explosive transition is presented in both homogeneous and heterogeneous structures for small and intermediate k^{*}'s. Otherwise, a large restriction can sustain a discontinuous transition only in the heterogeneous case. This shares some similarities with recent results for the Kuramoto model [Phys. Rev. E 91, 022818 (2015)PLEEE81539-375510.1103/PhysRevE.91.022818]. Surprisingly, intermediate restriction and large inertia are responsible for the emergence of an extra phase, in which the system is partially synchronized and the classification of phase transition depends on the inertia and the lattice topology. In this case, the system exhibits two phase transitions.

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.

A Mathematical Model for Storage and Recall of Images using Targeted Synchronization of Coupled Maps

We propose a mathematical model for storage and recall of images using coupled maps. We start by theoretically investigating targeted synchronization in coupled map systems wherein only a desired (partial) subset of the maps is made to synchronize. A simple method is introduced to specify coupling coefficients such that targeted synchronization is ensured. The principle of this method is extended to storage/recall of images using coupled Rulkov maps. The process of adjusting coupling coefficients between Rulkov maps (often used to model neurons) for the purpose of storing a desired image mimics the process of adjusting synaptic strengths between neurons to store memories. Our method uses both synchronisation and synaptic weight modification, as the human brain is thought to do. The stored image can be recalled by providing an initial random pattern to the dynamical system. The storage and recall of the standard image of Lena is explicitly demonstrated.

Dynamics of oscillators globally coupled via two mean fields

Many studies of synchronization properties of coupled oscillators, based on the classical Kuramoto approach, focus on ensembles coupled via a mean field. Here we introduce a setup of Kuramoto-type phase oscillators coupled via two mean fields. We derive stability properties of the incoherent state and find traveling wave solutions with different locking patterns; stability properties of these waves are found numerically. Mostly nontrivial states appear when the two fields compete, i.e. one tends to synchronize oscillators while the other one desynchronizes them. Here we identify normal branches which bifurcate from the incoherent state in a usual way, and anomalous branches, appearance of which cannot be described as a bifurcation. Furthermore, hybrid branches combining properties of both are described. In the situations where no stable traveling wave exists, modulated quasiperiodic in time dynamics is observed. Our results indicate that a competition between two coupling channels can lead to a complex system behavior, providing a potential generalized framework for understanding of complex phenomena in natural oscillatory systems.

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.