One node driving synchronisation.

Abrupt changes of behaviour in complex networks can be triggered by a single node. This work describes the dynamical fundamentals of how the behaviour of one node affects the whole network formed by coupled phase-oscillators with heterogeneous coupling strengths. The synchronisation of phase-oscillators is independent of the distribution of the natural frequencies, weakly depends on the network size, but highly depends on only one key oscillator whose ratio between its natural frequency in a rotating frame and its coupling strength is maximum. This result is based on a novel method to calculate the critical coupling strength with which the phase-oscillators emerge into frequency synchronisation. In addition, we put forward an analytical method to approximately calculate the phase-angles for the synchronous oscillators.

A remarkable phenomenon in phase-oscillator networks is the emergence of collective synchronous behaviour123456 such as phase synchronisation or phase-locking7891011. The Kuramoto model121314, a paradigmatic network to understand behaviour in complex networks, has drawn lots of attention of scientists1516171819. Many incipient works about Kuramoto model have assumed an infinite amount of oscillators coupled by a homogeneous strength. In 2000, Strogatz wrote20: “As of March 2000, there are no rigorous convergence results about the finite-N behavior of the Kuramoto model.” Since then, understanding the behaviour of networks composed by a finite number of oscillators2122232425262728 coupled by heterogeneously strengths2930 has been the goal of many recent works towards the creation of a more realistic paradigmatic model for the emergence of collective behaviour in complex networks.

However, most of the works about the finite-size Kuramoto model have relied on a mean field analysis, and consequently the emergence of synchronous behaviour has been associated with the collective action of all oscillators. Little is known about the contribution of an individual oscillator into the emergence of synchronous behaviour. But emergent behaviour in real complex networks can be tripped by only one node. Understanding the mechanism behind such a phenomenon in a paradigmatic, more realistic phase-oscillator network model is a fundamental step to develop strategies to control behaviour in complex systems. Besides, no analytical work has been proposed to solve the phase-angles of the synchronous oscillators. But a solution for the phase-angles is of great importance as, for example, they are key variables for monitoring generators in the power grids where a Kuramoto-like model is considered313233.

In this paper, we firstly provide a novel method to calculate the critical coupling strength that induces synchronisation in the finite-size Kuramoto model with heterogeneous coupling strengths. From our theory, we understand that the synchronisation of a finite number of oscillators is surprisingly independent of the distribution of their natural frequencies, weakly depends on the network size, but remarkably depends on only one key oscillator, the one maximising the ratio between its natural frequency in a rotating frame and its coupling strength. This lights a beacon for us that in order to predict, enhance or avoid synchronisation in a network of arbitrary size, all required is the knowledge of the state of only one node rather than the whole system. Under a practical point of view, if a pinning control3435 would be applied to enhance or slack synchrony in the studied network, the control function can be input into only one node. In addition, we put forward an analytical method to approximately calculate the phase-angles of synchronous oscillators, without imposing any restriction on the distribution of natural frequencies. This directly links the synchronous solution and the physical parameters in phase-oscillator networks.

Results

Software codes

All the software codes for this paper are available by searching at http://pure.abdn.ac.uk:8080/portal/

Critical coupling strength

We use to denote the N × 1 vector with all elements equal to one (zero), to indicate the index set . Given a vector with N elements, we use    to denote the mean value of the elements of . The finite-size Kuramoto model with heterogeneous coupling strengths for all-to-all networks is defined as,

where N > 0 is a finite integer number, K > 0 is the coupling strength, , , and , denote the vectors whose elements represent the oscillators’ natural frequencies, instantaneous phases, and coupling weights, respectively. Define the frequency synchronisation (FS), i.e., the phase-locking state, of the phase-oscillators described by Eq. (1) as,

Our goal is to find K, as the oscillators emerge into FS for a large enough K with as K > K.

Let , , indicate the instantaneous frequency of the oscillators when FS is reached. Divide by α on both sides of Eq. (1), then sum the equation from i = 1 to N, this results in . We rewrite Eq. (1) in a rotating frame, namely, let and , , such that as the oscillators emerge into FS, and we have,

Define the order parameter1213 by,

Multiplying e− on both sides of Eq. (4) and then equating its real part and imaginary part, respectively, we have

The mean field form of Eq. (3) is . Let , and , , such that, when FS is reached, i.e., , we have

Considering , where s(i) = ±1, we have, from Eqs. (5) and (7), that,

Define a function f as , where , and a set as representing the solution for the synchronisation manifold of Eq. (3). From Eqs. (6) and (7), we know, . Verwoerd and Mason26 proved that

This conclusion was obtained by a Kuramoto model with a mean field coupling strength, i.e., , , . However, the conclusion in (9) is still effective for the general case where α ≠ α. Because the proof for this conclusion was independent of , and the only restriction was 26, which is fulfilled when α ≠ α. The conclusion in (9) means that if Eq. (3) has at least one FS solution, then Eq. (8) holds with s(i) = 1, . This FS solution is obtained for , where K is the critical coupling strength for FS, which ensures that Eq. (8) holds with s(i) = 1, 26. Our following analysis is under the restriction that s(i) = 1, , which implies , i.e., , .

Define the key ratio by,

meaning that ζ is the one of ζ possessing the maximum absolute value. We call the m-th oscillator as the key oscillator. We assume ζ ≠ 0 by ignoring the particular case where ζ = 0 resulting in ω = 0 and ζ = 0, . Let x = sin ϕ, where x ≠ 0 and ϕ ≠ 0 obtained from ζ ≠ 0 and Eq. (7). Then we have, from Eq. (7), that . Substituting into Eq. (8), and considering s(i) = 1, , r can be calculated by

Because , and , , we have, from Eq. (7), that  , , implying , . Therefore, the m-th oscillator (the key oscillator) is the most “outside” one of all FS oscillators spreading on a unit circle, where the most inner oscillator possesses the smallest value of among all oscillators. As K is decreased from a larger value that enables FS in the network to a smaller one, (as well as ) increases correspondingly since from Eq. (7). For any i ≠ j, if , we have from Eq. (7), implying . This means that is determined only by the condition , and is independent of K. Thus, if we rank oscillators by their values of , this ranking is not altered as K is varied. This means that, regardless of the value of K, the key oscillator is always the most “outside” one. FS stops existing if no solution is found for , for any one oscillator. As K is decreased further, the first oscillator for which (and therefore, no solution is found for ) will be the key oscillator, because , , such that exceeds 1 at first. This means that K is the smallest K for which the key oscillator has a zero instantaneous frequency in the rotating frame, i.e., , resulting in Eq. (7) as i = m with restrictions  and ϕ ≠ 0. Therefore, K can be obtained by the following optimisation (OPT) problem in (12) to find the minimum K that implies with the restrictions that x ∈ [−1, 1] and x ≠ 0, where r is calculated by Eq. (11), namely,

where , if ω > 0, and , if ω < 0, where ε+ (ε−) indicates a positive (negative) infinitesimal. OPT in (12) can be numerically solved by selecting a small step for x, x, then increasing x from x to x by x, such that we get a series of values of f(x). The minimum f(x) is K.

Explosive synchronisation was studied in ref. 36 using a generalised Kuramoto model, which is a particular case of the model described in Eq. (3) by setting , . In this case, we have ζ = ±1, and from Eq. (7). Then OPT in (12) can be analytically solved, and the minimum of f(x) is 2, i.e., K = 2 when . This result remarkably coincides with the critical coupling strength proposed in ref. 36 for the backward process (namely, decrease K from a larger one to a smaller one) of the explosive behaviour. However, the critical coupling strength for the backward process is different from the one for the forward process (namely, increase K from a smaller one to a larger one) for the explosive synchronisation36. In this paper, we consider network configurations for which the critical coupling strength is the same for both the backward process and the forward process, i.e., no explosive synchronisation happens, then K obtained by OPT in (12) is also the critical coupling strength for the onset of FS in the forward process.

We further find, numerically, that OPT in (12) obtains its solution at . Consider , an approximate K can be analytically obtained by forcing , namely,

Let us now numerically demonstrate the exactness of the OPT in (12) to calculate K, and Eq. (13) to calculate K as the approximation of K, for different phase-oscillator networks. Let , , , where δ = 0 (δ > 0) indicates that all oscillators (not all oscillators) are in FS. The coupling weight α > 0, , is generated within110, without losing generality. Figure 1(a–c) show the results for three networks: Fig. 1(a), 10 oscillators with following an exponential distribution; Fig. 1(b), 50 oscillators with following a normal distribution; Fig. 1(c), 100 oscillators with following a uniform distribution. We calculate K by OPT in (12), and gradually decrease K from K = K + 0.2 to K  −  0.2. The results show that if K > K, δ = 0 with an acceptable error in numerical experiments for all cases, meaning that the oscillators are in FS. If K < K, δ > 0 implying that the oscillators lose FS for all cases. We note that the oscillators lose FS abruptly at K = K. This means that our method is effective to calculate K for all cases. Figure 1(d–f) demonstrate the effectiveness of Eq. (13) to analytically calculate an approximate K by forcing . Denote x as the value of x that provides K by OPT in (12). We define the relative error between 1 and as , and the relative error between K [Eq. (13)] and K [OPT in (12)] as . Figure 1(d–f) show the changes of η(x) and η(K) with respect to N(N = 3 to 200), with following exponential, normal and uniform distributions, respectively. The results indicate that K is near K, and is close to 1 for all cases. This means Eq. (13) works well to approximately calculate K for networks formed by arbitrary number of oscillators with any distributions.

Figure 1

(a–c) represent the results of δ (blue solid line) and K (red dash line) for networks formed by 10 oscillators with following an exponential distribution, 50 oscillators with following a normal distribution, and 100 oscillators with following an uniform distribution, respectively. (a), (b) and (c) are plotted based on average results of 5000 simulations with different initial phase-angles, but with the same and . (d–f) show the change of η(x) (blue line with circles) and η(K) (red line with triangles) for networks formed by N (N = 3 to 200) oscillators with following exponential, normal and uniform distributions, respectively. (d–f) are plotted based on average results of 100 simulations for each N, with different and .

One node driving synchronisation

Below, we show that K is independent of the distribution, weakly depends on the network size N, and mainly depends only on the key ratio of the key oscillator. For networks with different frequency distributions, diverse network sizes and various key ratios, we verify the dependence of K on the distribution, the network size N and the key ratio ζ. In order to present the results in a way such that they can be compared, we normalise ζ for these networks by making a parametrisation of α based on the value of ζ for each network. The surprising result is that, when we normalise ζ to be the same value for networks with different N and diverse distributions, K is roughly the same in these networks. Therefore, the key oscillator is the key factor for the behaviour of these networks. Next, we perform two sets of simulations to demonstrate this result. We use , and to denote the natural frequency vectors for networks constructed with a number of oscillators whose natural frequencies follow exponential, normal and uniform distributions, respectively, and correspondingly use , and to indicate the key ratios for these networks.

The first set of simulation includes 6 steps. (i), create all-to-all networks constructed by oscillators with natural frequencies , and , where N = 3 to 200. Thus, we have 3 * (200 − 2) = 594 networks in total, and each network has a key oscillator with a key ratio ζ. (ii), generate the coupling weights for all oscillators in the 594 networks by random numbers in [1, 10]. (iii), find the 594 key oscillators for the 594 networks, and create a set, , to contain all the 594 key ratios, i.e., , . (iv), find the maximum ζ in , mark it by ζ, and name this key oscillator as the “reference key oscillator” with label s. (v), change the values of α for all the key oscillators except for the reference key oscillator, such that all ζ are normalised as

where is a constant, and γ is a varying parameter which is set to be equal to 1 in the first set of simulation and will vary in the second set of simulation. Note that, this parametrisation process will enlarge all ζ except for ζ, such that all of these oscillators maintain their status of key oscillators in their own networks. (vi), calculate and record K for all the 594 networks.

In the second set of simulation, we further parametrise α as a function of γ for all the 594 key oscillators. We increase γ from its original value 1 to 20 by a small step, and simultaneously decrease each α by a proper ratio, such that Eq. (14) still holds. For each value of γ, we calculate and record K for all the 594 networks.

Figure 2(a–c) show the results for networks with frequency vectors given by , , and , respectively. The surfaces representing K are similar in all panels, which means that K is independent of the distribution. We note that K depends on N when N is small, but K is almost independent of N for most cases where . Thus, we say K weakly depends on N. However, if we keep N unchanged, we observe that K almost linearly increases with the growth of γ [i.e., the decrease of , and ] for all cases. In other words, K will increase if we decrease the coupling weight for only one key oscillator. The reason is that the key oscillator is the first one to lose FS when we decrease K, and a key oscillator with a smaller coupling weight is easier to lose FS, which in turn implies a larger K. As a conclusion, the behaviour of the key oscillator determines the FS of all oscillators, and the key ratio is the determinant physical parameter for the emergence of FS for all oscillators.

Figure 2

Exploring the determinant physical parameters for the emergence of the frequency synchronisation.

(a–c) show the results for networks formed by N(N = 3 to 200) oscillators with following exponential, normal, and uniform distributions, respectively. γ is the the parameter used to re-scale the key ratio. The surface represents the critical coupling, K, for different N and γ.

Master solution

When the oscillators emerge into FS, i.e., , the solution of Eq. (3) is

where is an arbitrary number, is the homogeneous solution of Eq. (3) by setting , and is a particular solution of the non-homogeneous Eq. (3). From Eq. (7), we have

where we exclude the unstable solutions for and for ζ < 0 (see Methods).

We name [Eq. (16)] as the master solution of Eq. (3), since it is an analytically expressible particular solution of Eq. (3), and it embodies all of other stable particular solutions, i.e., any stable particular solution can be expressed by . Note that r in Eq. (16) needs to be numerically calculated. Next, we propose an analytical method to approximately obtain the master solution.

Relabel the oscillators such that , and separate the oscillators into two groups: one group includes oscillators with labels from 1 to N′, where if N is even (or odd); the other group includes the remaining oscillators. Denote and for the first group and second group of oscillators, respectively. From Eqs. (6) and (7), we have . Thus, , implying . The non-negativity of μ1 comes from the fact that for any and any . Recall if N is even (or odd), we know if N is even (odd), implying if N is even (odd). For simplicity, we denote for both cases. When the oscillators emerge into FS with a given K′ (K′ > K), our model treats the whole system as two frequency-synchronous oscillators coupled by a common coupling strength K′, with natural frequencies μ1 and μ2, respectively. We assume that the two-oscillator system also follows the model described by Eq. (3) with coupling weights α1 = α2 = 1 which results in ζ1 = μ1, and ζ2 = μ2. Thus, from Eq. (7), we have

where r′ is the order parameter of the two-oscillator system. From Eq. (7), we have , where we exclude the case where (see Methods). Thus, we have from Eq. (5). Since , we have whose solution is

where and indicate a locally stable branch and a locally unstable branch of the FS solution for the two-oscillator model, respectively (see Methods). We only consider the stable branch (). Furthermore, we use the order parameter of the two-oscillator system to be an approximation of the order parameter [Eq. (5)] of the N-oscillator system, i.e., . Thus, the analytical approximation for the master solution in Eq. (16) is,

The corresponding approximate FS solution [Eq. (15)], is

Figure 3(a) shows the numerical results of the order parameter for a network with 50 oscillators where follows a normal distribution and α, , is a random number within110. K is indicated by the magenta dash-dot line. When , the approximate order parameter, λ1 [Eq. (18)] is close to the numerical one, r [Eq. (5)]. This means λ1 can effectively approximate r. Define an N × 1 vector, , with elements , representing the absolute error between [Eq. (19)] and [Eq. (16)]. Define as the standard deviation of . Figure 3(b,c) show the results of the average absolute error and σ, respectively, at K = K + 0.1 which ensures the emergence of FS. Networks are formed by N(N = 3 to 200) oscillators, with following exponential, normal and uniform distributions. and σ are small for all cases, which means that the error between and is small in all cases. Moreover, the larger K is, the smaller the error between λ1 and r is [Fig. 3(a)], which will further imply a smaller error between and ϕ, . This means our method is effective to solve the phase-angles for oscillators as they emerge into FS, for networks formed by an arbitrary number of oscillators with any distribution.

Figure 3

(a) The order parameter and its approximation for 50 oscillators. r (red line with triangles) is numerically calculated by Eq. (5) as s(i) = 1. . λ1 (blue line with circles) and λ2 (green line with squares) are calculated by Eq. (18) as . The value of K and 2μ1 are represented by magenta dash-dot line and black dash line respectively. (b,c) show, for different networks, the change of the average absolute error between in Eq. (19) and in Eq. (16) and the standard deviation (σ) of  , as a function of K, respectively. Networks with N (from 3 to 200) oscillators with following exponential (dash red line), normal (green solid line) and uniform (black line with “+”) distributions, respectively.

Discussion

In this paper, we presented our studies on the synchronisation for a finite-size Kuramoto model with heterogeneous coupling strengths. We provided a novel method to accurately calculate [OPT in (12)] or analytically approximate [Eq. (13)] the critical coupling strength for the onset of synchronisation among oscillators. With this method, we find that the synchronisation of phase-oscillators is independent of the natural frequency distribution of the oscillators, weakly depends on the network size, but highly depends on only one node which has the maximum proportion of its natural frequency to its coupling strength. This helps us to understand the mechanism of “the one affects the whole” in complex networks.

In addition, we put forward a method to approximately calculate the phase-angles for the oscillators when they emerge into synchronisation. With our method, one can easily obtain the solution of phase-angles for frequency-synchronous oscillators, without numerically solving the differential equation.

Methods

Excluding the unstable solutions

The FS solution of Eq. (3), i.e., the solution of Eq. (7) is

A rigorous analysis for the stability of the FS solutions was given by ref. 28 for a mean filed coupled Kuramoto model, i.e., α = α = 1, , . From the conclusion of ref. 28, we know that the FS solution of Eq. (3) is locally unstable if at least one s(i) = −1 in Eq. (8). In other words, if the FS solution is stable, then s(i) = 1, , implying that , i.e., , . Therefore, we exclude the case that for the solution of the two-oscillator system in the paper.

However, the stability analysis of the FS solution for the general case where α ≠ α is difficult and is still an open problem. In our numerical experiments, the stable solution we obtained is only the one that , . Thus, we exclude the solutions that for , and that for ζ < 0.

The stability analysis for the two-oscillator system

The two-oscillator system also follows the Kuramoto model with α1 = α2 = 1, namely,

Let be a FS solution of Eq. (22). Let , where is a small perturbation on . Linearise Eq. (22) around , we have,

where the Jacobian matrix J is

The two eigenvalues of J are e1 = 0 and . If the FS solution is stable, we have e2 < 0 implying , i.e., .

We have, from Eq. (18), that . Substituting this condition into Eq. (17), we get and . If , we have that and . Because from Eq. (18), we approximately have and , implying . If λ1 grows larger as K increases from 2μ1, will become larger. However, implies instability of the FS solution of the two oscillators. This means that r′ ≈ λ2 describes an unstable FS solution. On the other hand, r′ ≈ λ1 ensures the stability of the FS solution.

Additional Information

How to cite this article: Wang, C. et al. One node driving synchronisation. Sci. Rep. 5, 18091; doi: 10.1038/srep18091 (2015).