Skip to main content
Scientific Reports logoLink to Scientific Reports
. 2013 Feb 15;3:1281. doi: 10.1038/srep01281

Explosive transitions to synchronization in networks of phase oscillators

I Leyva 1,2, A Navas 2, I Sendiña-Nadal 1,2,a, J A Almendral 2, J M Buldú 1,2, M Zanin 2,3,4, D Papo 2, S Boccaletti 2
PMCID: PMC3573336  PMID: 23412391

Abstract

The emergence of dynamical abrupt transitions in the macroscopic state of a system is currently a subject of the utmost interest. The occurrence of a first-order phase transition to synchronization of an ensemble of networked phase oscillators was reported, so far, for very particular network architectures. Here, we show how a sharp, discontinuous transition can occur, instead, as a generic feature of networks of phase oscillators. Precisely, we set conditions for the transition from unsynchronized to synchronized states to be first-order, and demonstrate how these conditions can be attained in a very wide spectrum of situations. We then show how the occurrence of such transitions is always accompanied by the spontaneous setting of frequency-degree correlation features. Third, we show that the conditions for abrupt transitions can be even softened in several cases. Finally, we discuss, as a possible application, the use of this phenomenon to express magnetic-like states of synchronization.


Many complex systems operate transitions between different regimes or phases under the action of a control parameter. These transitions can be monitored using a global order parameter, a physical quantity (e.g. scalar, vector, …) accounting for the symmetry of the phases. Phase transitions can be of first or second order according to whether the order parameter varies continuously or discontinuously at a critical value of the control parameter. In complex networks theory1, phase transitions have been observed in the way the graph collectively organizes its architecture (e.g. percolation2,3) and dynamical state (e.g. synchronization4,5,6).

Recently, Achlioptas et al.7 proposed a slight modification of the classical percolation models, leading to a first-order type transition, named ”explosive percolation”. In the model, the authors consider a preferential growth process for a network with fixed number of vertices. Edges are added one by one, by picking each time two randomly selected possible edges, and deciding to keep the one having a lower value of the product of the size of the two components that edge is joining, this way delaying as much as possible the formation of a giant connected component. At the percolation threshold, a macroscopic part of the system (a giant cluster of size O(N)) becomes connected. While the transition is found to be abrupt in Erdös-Renyi (ER) networks, heterogeneous structures as scale-free (SF) networks not always display explosive percolation, which crucially depends, instead, on the degree distribution8,9. The continuous or discontinuous nature of explosive percolation is still a matter of debate10.

The existence of abrupt transitions has also been investigated for the synchronization of networked phase oscillators5,6. In this context, a first-order transition has initially been described in the Kuramoto model11 with a particular realization of a uniform frequency distribution (evenly spaced frequencies) and an all-to-all network topology12. However, lately it has been shown that a positive correlation between the node degree and the corresponding oscillator's natural frequency may lead to a first-order transition also for SF networks13. These latter results have been further generalized to networked chaotic oscillators, and experimentally proved14.

In this report, we show that a sharp, discontinuous phase transition is not restricted to the above rather limited and apparently opposite cases, but it constitutes, instead, a generic feature of the synchronization of networked phase oscillators. Precisely, we initially give a condition for the transition from unsynchronized to synchronized states to be first-order, and demonstrate how such a condition is easy to attain in many circumstances, and for a wide class of oscillators' initial frequency distributions. We then show how such transitions are always accompanied by the spontaneous emergence of frequency-degree correlation features, and discuss how the considered condition can even be softened in several cases. Finally, we illustrate, as a possible application, the option of expressing magnetic-like states of synchronization with the use of such transitions.

Results

Our results are obtained with a network of N Kuramoto11 oscillators, which is described by:

graphic file with name srep01281-m1.jpg

where φi is the phase of the ith oscillator (i = 1, …, N), ωi is its associated natural frequency [drawn from a generic frequency distribution p(ω)], d is the coupling strength, and {aij} are the elements of the adjacency matrix that uniquely defines the nodes' interactions.

The classical order parameter for the system given by equation (1) is Inline graphic, and the level of phase synchronization can be monitored by looking at the value of S = 〈r(t)〉T, where 〈…〉T denotes a time average with Inline graphic. Furthermore, for each oscillator i, we denote by Inline graphic the set of oscillators linked to it.

As the coupling strength d increases, system (1) undergoes a phase transition from the unsynchronized Inline graphic to a synchronous (S ~ 1) state, where all oscillators ultimately acquire the same frequency. For this phase transition to display a first-order feature, we have to avoid that any oscillator behave as the core of a clustering process, where its neighbors begin to aggregate to the synchronous state smoothly and progressively, as in the classical routes described in Ref.15.

Along this paper, we realize such a condition by explicitly imposing certain constraints in the frequency differences between each node i and the whole set Inline graphic of oscillators belonging to its neighborhood. A first general constraint can be set as follows:

(C1) for each oscillator i, all nodes j belonging to Inline graphic satisfy | ωiωj |> γc.

Notice that condition (C1) is tantamount to imposing a minimal value for the frequency difference between linked nodes of the network.

We now fix N = 500, and illustrate the synchronization route for several frequency distributions p(ω), when condition (C1) is met in a frequency gap-conditioned (FGC) random network (see Methods section for the details on the network construction). In particular, we use the adjacency matrix obtained after the FGC network is constructed to simulate system (1), and monitor the state of the network by gradually increasing the value of d in steps δd = 10–4, from d = 0. At each step, a long transient is discarded before the data are acquired for further processing. Moreover, insofar as we are looking for a first-order phase transition (and thus for an expected associated hysteresis), simulations are also performed in the reverse way, i.e. starting from a given dmax (where the network is phase synchronized), and gradually decreasing the coupling by δd at each step. In what follows, the two sets of numerical trials are termed forward and backward, respectively.

Figure 1 reports the results obtained by setting p(ω) as a uniform frequency distribution in the interval [0, 1]. Panels a and b of Fig. 1 show S as a function of d. In particular, panel a (b) illustrates the case of a fixed mean degree 〈k〉 = 40 (of a fixed frequency gap γ = 0.4), and reports the results of the forward and backward simulations at different values of γ (〈k〉). A first important result is the evident first-order character acquired, in all cases, by the transition for sufficiently high values of γ.

Figure 1. Explosive transition to synchronization.

Figure 1

Phase synchronization level S (see text for definition) vs. the coupling strength d, for different values of the gap γ at 〈k〉 = 40 (panel a), and for different values of the average degree 〈k〉 at γ = 0.4 (panel b). In both panels, the continuous (dashed) lines refer to the forward (backward) simulations. See the Methods section for the construction procedure of the networks. In panels c and d, the degree ki that each node achieves after the network construction is completed (upper plots) and the average of the natural frequencies 〈ωj〉 for Inline graphic (bottom plots) are reported vs. the node's natural frequency ωi, for 〈k〉 = 100 and frequency gaps γ = 0.0 (panel c), and γ = 0.4 (panel d). The red solid line in panel d is a sketch of the theoretical prediction f(ω) (see text). Panels e and f show S (color coded according to the color bar) in the parameter space (d, γ) for (e) 〈k〉 = 20 and (f) 〈k〉 = 60. The horizontal dashed lines mark the separation between the region of the parameter space where a second-order transition occurs (below the line) and that in which the transition is instead of the first order type (above the line). The yellow striped area delimits the hysteresis region.

A second relevant result is the spontaneous emergence of degree-frequency correlation features associated to the passage from a second- to a first-order phase transition. While such a correlation was imposed ad hoc in Refs.13,14, here condition (C1) creates for each oscillator i a frequency barrier around ωi, where links are forbidden. The final degree ki is proportional to the total probability for that oscillator to receive connections from other oscillators in the network, and therefore to Inline graphic. This is shown in panels c and d of Fig. 1, where the degree ki that each node achieves after completing the FGC network construction is reported as a function of its natural frequency ωi, for 〈k〉 = 100. Precisely, the upper plot of panel c refers to the case γ = 0 in which no degree-frequency correlation is present. In the upper plot of panel d, instead, we report the case γ = 0.4 (a value for which a first-order phase transition occurs) and the (conveniently normalized) function Inline graphic, with p(ω) = 1 for ω ∈ [0, 1], and p(ω) = 0 elsewhere, which gives evidence of the emergence of a very pronounced V-shape relationship between the frequency and the degree of the network's nodes. By further inspecting the average frequency of each oscillator's neighborhood, comparison between the lower plots of panels c and d manifests that condition (C1) leads to the emergence of a bipartite-like network where low frequency oscillators are mainly coupled to high frequency oscillators.

Finally, panels e and f of Fig. 1 report S in the dγ space, for 〈k〉 = 20 and 〈k〉 = 60, and show that the rise of a first order phase transition is, indeed, a generic feature in the parameter space. The horizontal dashed lines in panels e and f mark the values of γc, separating the two regions where a second-order transition (below the line) and a first-order transition (above the line) occurs. Remarkably, while increasing 〈k〉 of the network facilitates the occurrence of the explosive transition, as both the values of γ and d for which the first-order phase transition takes place decrease, it also shrinks the width of the hysteresis (the yellow striped area in panels e and f). This is consistent with the results of Ref.12 for an all-to-all connected case, where in the limit N → ∞ a first-order phase transition has been predicted in the absence of hysteresis. It is worth noticing that similar scenarios of synchronization were observed in the past for chains of periodic16 and chaotic phase coherent17 oscillators. In this latter cases, synchronization was shown to appear or vanish in two ways: a soft transition (without cluster formation) for chains with very small frequency mismatches between the elements, and a hard transition for rather long chains with relatively large frequency mismatches.

We verified that fulfillment of condition (C1) leads to a first-order transition for a very wide class of distributions of the oscillators' natural frequencies. The panels a and b of Fig. 2 show, for instance, the case of the asymmetrical Rayleigh distribution, conveniently re-scaled to the interval [0, 1], given by Inline graphic (see the inset in panel a of Fig. 2), with σ = 0.25. The results highlight the presence of an abrupt transition (Fig. 2, panel a), together with the emergence of a clear frequency-degree correlation, well described, again, by the function f(ω) (Fig. 2, panel b).

Figure 2. Extension to different frequency distributions and network construction rules.

Figure 2

(Top row) S vs. d resulting from the forward (continuous lines) and backward (dashed lines) simulations of system (1) for different frequency distributions or network construction rules, and (bottom row) the corresponding distribution of the final node degree ki vs. the corresponding oscillator's natural frequency. (a)-(b) Rayleigh distribution for γ = 0.3. In panel b, the red solid line depicts the theoretical prediction f(ω); (c)-(d) uniform frequency distribution, but network constructed accordingly to the local mean field condition (see text) for γ = 0, and γ = 0.4. In panel d γ = 0.4; (e)-(f) Gaussian distribution with Z = 0.7. The insets in panels a and e report the corresponding distribution p(ω). See the text and the Methods section for the details on the specific construction procedure used in each case. In all cases, 〈k〉 = 60.

We now move to discussing several ways for even softening condition (C1), thus making it even more generally applicable, while still keeping the abrupt character of the transition. For the first extension, we consider again the case of a uniform frequency distribution in the interval [0, 1]. In this case, condition (C1) can be relaxed as follows:

(C2) for each oscillator i, all nodes j belonging to Inline graphic satisfy |ωi – 〈ωj〉| > γc, where 〈…〉 indicates the average value over the ensemble Inline graphic.

The new condition is tantamount to softening condition (C1) to the local mean field of frequency differences in the neighborhood of each network node. The panels c and d of Fig. 2 report the results of networks obtained with a modified construction procedure, in which pairs of randomly selected nodes are now linked if the value of |ωi – 〈ωj〉| (averaged over the set of nodes j already linked to node i, and the one candidate to be further linked), exceeds a gap γ. Again, an explosive transition occurs (Fig. 2, panel c) in correspondence with the emergence of frequency-degree correlations (Fig. 2, panel d).

Furthermore, it is worth noticing that a strict application of condition (C1) for non uniform frequency distributions implies that oscillators at different frequencies would in general have a different number of available neighbors in the network. For this reason, a natural extension of condition (C1) is to consider a frequency-dependent gap γ(ω) defined by Inline graphic. Panels e and f of Fig. 2 report the case of a Gaussian distribution limited to the frequency range [0, 1], centered at ω = 0.5, and given by Inline graphic, with σ = 0.13. The gap condition for the construction of the network is now to fix the value of Z, and accept the pairing of nodes when Inline graphic. The generic case is again that an explosive transition is obtained, with pronounced frequency-degree correlation features, as long as p(ω) is symmetrical.

An alternative way to explore the frontier between first- and second-order transitions is to start from the case of Ref.12 (an all to all network configuration with evenly spaced frequencies of the oscillators), and to study the robustness of explosiveness under removal of links by means of two pruning processes: a random pruning, where links are randomly chosen and removed, and a preferential pruning that removes links with a probability proportional to the inverse of the frequency difference of the connected nodes, so that close frequency oscillators are more probable to be disconnected (see the Methods section for details).

The results are shown in Fig. 3. For the random pruning process (panel a of Fig. 3), one can easily see that, as the connectivity of the network becomes low enough (i.e. as more and more links are pruned in a way that 〈k〉 progressively decreases) the transition gradually passes from first- to second-order. On the contrary, for the preferential pruning process (panel b of Fig. 3), the system maintains the explosive behavior due to the higher probability to uphold connections associated to large frequency mismatches, although the size of the discontinuity decreases. Once again, inspection of the nodes' degree and of the mean frequencies of the nodes' neighborhoods (panel c of Fig. 3) reveals the spontaneous emergence of correlation between topology and dynamics for the preferential pruning process.

Figure 3. Explosive synchronization in networks with evenly spaced natural frequencies.

Figure 3

S vs. d resulting from the forward (continuous lines) and backward (dashed lines) simulations of system (1), for different values of the average degree 〈k〉 (see legend in panel a) after performing a random pruning (a) and a preferential pruning (b) of links in an all-to-all connected network (see the text and the Methods section for the networks construction). Natural frequencies are evenly spaced in the interval [0, 1]. In panel c, the degree ki that each node achieves after the network pruning is completed (upper plots) and the average of the local natural frequency 〈ωj〉 for Inline graphic (bottom plots) are reported vs. the node's natural frequency ωi for the preferential pruning process with 〈k〉 = 10.

Discussion

In conclusion, we demonstrated the validity and generality of several conditions for the occurrence of abrupt phase transitions in networks of phase oscillators. Our study generalizes all previous results, and extends the possibility of encountering first-order phase transitions to a large variety of network topologies, as well as to a large variety of frequency distributions of the oscillators. This suggests practical methods for engineering networks able to display critical phenomena, and the emergence of dynamical abrupt transitions in their macroscopic states. Furthermore, the evidence for the emergence of frequency-degree correlations in connection with these abrupt transition, may shed light on the mechanisms underlying the relationship between topology and dynamics in many real-world systems.

One interesting application is expressing magnetic-like states of synchronization in such an ensemble of networked oscillators, provided that the coupling strength is set inside the hysteresis region of the first-order phase transition. For this purpose, one can again consider the case of an initial uniform distribution of the oscillators' frequencies, and modify system (1) as follows: Inline graphic, where Dp is the strength of a unidirectional connection to an external pacemaker (equal for all oscillators in the network), and φp is the phase of the pacemaker obeying Inline graphic. Initially, the system is let to freely evolve into the unsynchronized regime, with Dp = 0. The pacemaker is then switched on, and Dp is selected so that all oscillators are entrained to the pacemaker phase. At a subsequent time, the pacemaker is switched off again, and the state of the system is monitored. The results are illustrated in Fig. 4. Precisely, panel a shows that setting d outside (inside) the hysteresis region produces a final state that relaxes to the original unsynchronized behavior (that stays permanently in a synchronized configuration for sufficiently large Dp values). Panel b depicts the regions of the parameter space Dpωp for which these magnetic-like states can ultimately be produced.

Figure 4. Magnetic-like states of synchronization.

Figure 4

(a) Time evolution of the parameter r(t) under a pacemaker forcing for a value d = 0.004 outside the hysteresis region (upper plot), and a value d = 0.009 within the hysteresis region (bottom plot). ωp = 1.0 and Dp = 0.0005 (bottom red line), Dp = 0.005 (middle blue line), and Dp = 0.02 (top black line). The pacemaker is active from t = 50 to t = 350, as marked by the vertical dashed lines. (b) Colormap of S = 〈r(t)〉t>350 (coded as indicated in the color bar), showing the region of the parameter space Dp-ωp where the magnetic-like state of synchronization is maintained after removal of the pacemaker. The initial frequencies of the oscillators are taken from a uniform distribution in the interval [0, 1]. γ = 0.49 and 〈k〉 = 40.

Methods

Algorithms for the networks' construction

Frequency Gap-conditioned (FGC) random networks. The considered networks result from the following procedure: i) we assign natural frequencies {ωi}, drawn from a distribution p(ω), to the N oscillators; ii) we randomly pick a pair (i, j) of oscillators, and form a link between them only if the value of |ωiωj| exceeds a given gap γ; iii) we repeat point ii) until the desired number of links L in the graph is formed. After a final check on the connectedness of the resulting network, the procedure yields Erdös-Renyi-like topologies with an average degree Inline graphic.

Random and preferential pruning. The considered networks result from the following procedure: i) we start from an all-to-all network configuration, and assign to the N oscillators natural frequencies ωi evenly spaced spanning the interval [0, 1], i.e. Inline graphic; ii) from the total number of links Inline graphic, Inline graphic links are pruned to produce a network with desired mean degree 〈k〉. The links removal can either be performed randomly, or selectively by assigning to each link a pruning probability Inline graphic.

Author Contributions

IL, AN, ISN and SB devised the model and designed the study. IL and ISN carried out the numerical simulations. IL, ISN, JAA, JMB, DP, and MZ analyzed the data and prepared the figures. AN, MZ, DP, and SB wrote the main text of the manuscript.

Acknowledgments

Work supported by Ministerio de Educación y Ciencia, Spain, through grants FIS2009-07072 and from the BBVA-Foundation within the Isaac-Peral program of Chairs. Authors acknowledge also the R&D Program MODELICO-CM [S2009ESP-1691], and the usage of the resources, technical expertise and assistance provided by supercomputing facility CRESCO of ENEA in Portici (Italy).

References

  1. Dorogovtsev S. N., Goltsev A. V. & Mendes J. F. F. Critical phenomena in complex networks. Rev. Mod. Phys. 80, 1275 (2008). [Google Scholar]
  2. Cohen R., Ben-Avraham D. & Havlin S. Percolation critical exponents in scale-free networks. Phys. Rev. E 66, 036113 (2002). [DOI] [PubMed] [Google Scholar]
  3. Li G., Braunstein L. A., Buldyrev S. V., Havlin S. & Stanley H. E. Transport and percolation theory in weighted networks. Phys. Rev. E 75, 045103 (2007). [DOI] [PubMed] [Google Scholar]
  4. Boccaletti S., Kurths J., Osipov G., Valladares D. L. & Zhou C. S. The synchronization of chaotic systems. Phys. Rep. 366, 1 (2002). [Google Scholar]
  5. Boccaletti S., Latora V., Moreno Y., Chavez M. & Hwang D. U. Complex networks: Structure and dynamics. Phys. Rep. 424, 175 (2006). [Google Scholar]
  6. Arenas A., Díaz-Guilera A., Kurths J., Moreno Y. & Zhou C. S. Synchronization in complex networks. Phys. Rep. 469, 93 (2008). [Google Scholar]
  7. Achlioptas D., D'Souza R. M. & Spencer J. Explosive percolation in random networks. Science 323, 1453 (2009). [DOI] [PubMed] [Google Scholar]
  8. Cho Y. S., Kim J. S., Park J., Kahng B. & Kim D. Percolation transitions in scale-free networks under the Achlioptas process. Phys. Rev. Lett. 103, 135702 (2009). [DOI] [PubMed] [Google Scholar]
  9. Radicchi F. & Fortunato S. Explosive percolation in scale-free networks. Phys. Rev. Lett. 103, 168701 (2009). [DOI] [PubMed] [Google Scholar]
  10. Grassberger P., Christensen C., Bizhani G., Son S.-W. & Paczuski M. Explosive percolation is continuous, but with unusual finite size behavior. Phys. Rev. Lett. 106, 225701 (2011). [DOI] [PubMed] [Google Scholar]
  11. Kuramoto Y. Chemical oscillations, waves and turbulence (Springer, 1984).
  12. Pazó D. Thermodynamic limit of the first-order phase transition in the Kuramoto model. Phys. Rev. E 72, 046211 (2005). [DOI] [PubMed] [Google Scholar]
  13. Gómez-Gardeñes J., Gómez S., Arenas A. & Moreno Y. Explosive synchronization transitions in scale-free networks. Phys. Rev. Lett. 106, 128701 (2011). [DOI] [PubMed] [Google Scholar]
  14. Leyva I., Sevilla-Escoboza R., Buldú J. M., Sendiña-Nadal I., Gómez-Gardeñes J., Arenas A. et al. Explosive first-order transition to synchrony in networked chaotic oscillators. Phys. Rev. Lett. 108, 168702 (2012). [DOI] [PubMed] [Google Scholar]
  15. Gómez-Gardeñes J., Moreno Y. & Arenas A. Paths to synchronization on complex networks. Phys. Rev. Lett. 98, 034101 (2007). [DOI] [PubMed] [Google Scholar]
  16. Osipov G. V. & Sushchik M. M. Synchronized clusters and multistability in arrays of oscillators with different natural frequencies. Phys. Rev. E 58, 7198 (1998). [Google Scholar]
  17. Osipov G., Pikovsky A., Rosenblum M. & Kurths J. Phase synchronization effects in a lattice of nonidentical Rssler oscillators. Phys. Rev. E 55, 2353 (1997). [Google Scholar]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group

RESOURCES